QUARTERLY OF APPLIED MATHEMATICS 





Vol. XVI JANUARY, 1959 No. 4 





SOME APPLICATIONS OF DIFFERENTIAL GEOMETRY TO THE 
INTERPRETATION OF PHYSICAL PHENOMENA* 


BY 
P.O. BELL 


Lockheed Aircraft Corporation, Missle Systems Division 


1. Introduction. ‘The motivation of the present investigation was the conviction 
that a differential geometric analysis of such physical phenomena as heat conduction 
in a homogeneous medium and electrical charge distribution on closed conductors 
should yield interesting quantitative geometric interpretations. Although the study 
has been somewhat rudimentary, the results have been, in the author's opinion, of 
sufficient interest to warrant extending them to their mathematical analogues in n-dimen- 
sional Euclidean space. For the sake of elegance of presentation, the three-dimensional 
results will appear as specializations of the n-dimensional theory. 

Let ¢, ¥ denote single-valued point functions of class C* in a domain D of a Euclidean 
space FR, of n-dimensions. Let V,-, denote the hypersurface on which @ is constant 
which passes through a given point P of D. Let the family of such hypersurfaces which 
pass through points of D be denoted by Ff’. Let o denote are length of the orthogonal 
trajectory 7 of F at P, the positive and negative senses of ¢ along 7 being taken to be 
those in which @ decreases and increases, respectively. The positive function Q defined 
by the relation 


Q = —dd/da (1.1) 


will be called the “gauge function” of F at P. A locus of points on V,_, at which Q has 
a constant value is an nm — 2 dimensional variety V,_, which will be called an “isogauge”’ 
of V,_, . An isogauge of V,_, with respect to F is, in general, an n — 2 dimensional 
variety. It can be n — | dimensional only if it coincides with V,_, . In such a case, F 
is a family of concentric hyperspheres, each V,_, being defined by a constant value of 
a differentiable function of the radius. 

Let M denote, the mean curvature of V,_, at P, and let p and Ky denote the unit 
principal normal and curvature, respectively, of 7 at P. Let Ay, A, denote the Laplacians 
of ¥ with respect to R, and V,_, respectively. The following theorems will be proved. 


Theorem 1. The relation which exists between the two Laplacians Ay, A, is given by 


iu. - ov + Ay — Kp-Vy, (1.2) 


fered Oo 


in which p- Vy is the projection of the gradient of y on the principal normal of T at P. 
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Corouuary (1.1). 7f @ is harmonic 


0 
M — (In Q), (1.3) 
Oo 
in which Q = —0¢/dc, and M is the mean curvature of V,_, at P (a result due to H. Wey! 


[5, p. 181)}). 
T heore m A Let r denote the orthogonal trae clory on V,, i al P of the family of rsogauge Ss 
of V,_, , and let s denote its arc length. The principal normal p of T at P is tangent to r, 


and the curvature of T at P is given by the relation 


. ra) 
Kr = —(in Q), (1.4) 
os 
An application to steady-state heat flow may be described as follows. Let F be a 
family of isothermal surfaces with respect to heat flow in a homogeneous isotropic 
ce. The fundamental hypothesis in steady-state heat flow states 


medium in ordinary spa 
that heat flows orthogonally across the isothermal surfaces at a rate proportional to 
—0¢/dc, in which ¢ denotes temperature. It follows that Corollary (1.1) and Theorem 2, 
in combination, yield the following interesting result. 

Theorem 3. In the case of a steady-state heat flow in a homoge neous tsotroprc medium 


the ratio of the rates of flow at any two points P, , P, of the medium is given by the relation 
4 


2./ = cap { - | M do + | Kk. ds), (1.5) 


\ va ‘ 


+ 


an which the first ante gral 1s along the arc of T from P,, to the point P of intersection of T 


with the isothermal surface which passes through P, and the second is along the arc of 7 
from P to the intersection Q, of + with the “isogauge’’ which passes through P, 


It is known that on a closed charged conductor S, of class C~ the density 6 of charge 
f no other charges than that on the conductor are present, is pro- 
Irom corollary (1.1) it follows that 6 is given by the relation 


at a point P, of S 
portional to (—0¢/dc) 


6 = kQ, exp | — M do, (1.6) 
in which © ts the value of the gauge function at the point P, to which o, corresponds, and the 


f 


integration ts from P, to P, along T. 


It is known [4, p. 191] that electrical charge placed on an ellipsoidal conductor Ly 
becomes distributed so that its density at a point P, of E, is proportional to the “sup- 
porting function” of FE, at P, , that is, the distance h, from the center of 2, to the tangent 
plane to E, at P, . The equipotential surfaces induced by the charged ellipsoid are 
confocal ellipsoids. Let HE denote an equipotential ellipsoid whose potential is given 


by ¢ = c. At points P, , P of E, , E, respectively, the gauge functions are given by 

Qo = koho , 0, = kh, (1.7) 
in which ky , k, are constants, and Q is proportional to the density of charge of FE, at P. 
On putting £, , E in the roles of S, , S for interpretation of (1.6) relative to confocal 
ellipsoids and making use of (1.7), the following relation is obtained. 
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mh,/h = exp (/ — M dc), (1.8) 


in which m ky /k. is a constant which corresponds to the selection Hy , E. 

The simple yet fundamental role of the “supporting function” in characterizing 
physical phenomena associated with confocal ellipsoids suggests its usefulness in a 
general study of isothermal surfaces. Such a study is presented in the present paper in 
which use is made of the supporting function together with an associated ‘natural’ 
coordinate system. The representation employed may be described as follows. Let a 
one-parameter family of surfaces of class C° be defined by a functional relation, 


h = h(6,w, >), @ = constant, (1.9) 


in which relative to a point P(6, w, @) the numbers h, 6, w, @ have the following meanings. 
One surface S on which ¢ is constant passes through P. The numbers h, 6, w are spherical 
coordinates, with respect to some fixed rectangular cartesian coordinate frame, of the 
foot of the perpendicular line from the origin to the tangent plane to S at P; that is 
to say, h is the supporting function of S at P. A curve described by a point P as ¢ varies, 
but @ and w are held fixed, is called a ¢-curve. A curve is a ¢-curve relative to F if, and 
only if, at its points the tangent planes to F are parallel. 

Let a system F of surfaces, not necessarily isothermal, be represented by an equation 
of form (1.9). At a point P of a surface S of the system, let y denote the acute angle 
between the normal to S at P and the tangent to the ¢-curve at P, and let K denote 
the curvature at P of the curve of section of S by that normal plane of S which contains 
the tangent to the ¢-curve at P. The following relations will be established 


h, dhs = hy, + (K tan’ y)hi , (1.10) 
oles 
‘ dh, 
h; Ao + Mh, + — = 0, (1.11) 
where 
h, = 0h/00, h, = dh/dw, hz = dh/dd. 


If @ is a harmonic function, relations (1.10), (1.11) reduce to 


Oh, 
Oo 


hy, + hi(M + K tan’ y) = 0, + Mh, = O. (1.12) 

The coefficients of Eqs. (1.12) are geometric quantities whose values are independent of 
the choice of coordinates. The quantity —1/h, will be shown to be equal to the gauge function 
Q of F at P. It follows that the first equation of (1.12) is equivalent to the following 


first-order equation in @ and Q 


AQ , ' 

— + M + K tan’ y = 0. (1.13) 
0d 

This equation leads to the following integral equation which defines the difference 


between the gauge functions at two points P, , P, of a ¢-curve 


— a [ "(M + K tan’ y) dé. (1.14) 


“@e 
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An interesting physical interpretation of (1.14) may be stated as follows 


Theorem 4. In steady-state heat flow through a homogeneous isotropic medium, the 

di fie rence between the rates of flow Po ; p; al points a A spectively of a flow line is 
given by a relation of the form 

a-p=k | (M + K tan? y) dé, (1.15) 

in which k is a positive constant and ¢, , ¢, are temperatures at P,, , P, respectively, and the 


integration is along the flow line. 


The results stated above may be specialized to obtain the analogous geometric 
interpretations of physical phenomena applicable to the systems of cylinders which 
constitute the integral surfaces of the Laplace equation in two variables. 

2. The Laplacians A and A, of a function 1 with respect to the surfaces @ = constant. 
In a domain D of a Euclidean space RP, let ¢, ¥ denote single-valued point functions of 
class C", let the family F of hypersurfaces V,_, (introduced in See. 1) be coordinate 
hypersurfaces x” = ¢ = c whose orthogonal trajectories 7' are the x"-curves of PR, 
Let V,_; be referred to an orthogonal system of coordinate curves. The fundamental 


metric then assumes the form’ 
ds’ = g,,(dx')’ = g,,(dx")” + gaa(dx")’. (2.1) 


Theorem 1 will be proved and some specializations will be noted in this section. 
The Laplacian of a function y with respect to PR, at P is given by the well-known 


formula 





; Oy a» OW rw 
? Or OX OX 
in which the functions g'’ are defined by relations of the form 
E eo -& 1 
gg 6 2.3) 
l, k#1 
and the functions I’; are the Christoffel functions of the second kind defined by the 
relations 
o — a O”d,,; Od O”d;; ’ 
J g (27, k] f ( ee ee err 3 (2.4) 
: . Ox OX Ox 
In view of the conditions of orthogonality 
g;; = 0, 1A~Y (2.5) 
and Eqs. (2.3) and (2.4), Eq. (2.2) assumes the form 
P oy ae ov = ow a OY i a 
Ay qd ~~ | Gale ] wna & J a : + Aiy, 2.6) 
ar or Or Ox’ 
1Greek indices will have the range 1,2,...,n-1. Latin indices (except n) will have the range 1,2,...,7 


Repeated indices in a term denote the usual summation convention over the respective ranges. 
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in which A,y is defined by 
i , @O , 
Ay g"(- “ v Tt | ee 4 (2.7) 
Ox OX 


ox” 


and the Christoffel functions are given by the formulas 


I” “ ) — 7 9G» 
on ( win, We <6. 9 
9 2 ax 
= a Q”” O9as 6 
ier gq’ lap, | alae terres (2.8) 
, 2 ox 
a va Gg” OGun 
ca q Inn, w| ie : 
1" ( o . 
= Ov 


The Laplacian of Y may, therefore, be written in the form 


ay , |g cs OUae oy sal Oy 


Ay q”" { g’"T!, —: + Ay (2.9) 
(dx) 2 Ol On , OX 
Phe following relations will be established 
ie . O”d., aw av 
Py Bac OY _ yy OV (2.10) 
2 Or Ox da 
ay - ay ay , ay 
—— Ky = ( ~ (2.11) 
Oo OX Ax" OX 


in which p' denotes the contravariant components of the unit principal normal of T at P. 
Let ¢ denote the unit normal to V,_, at P, and let ¢' denote its contravariant com- 
ponents. The following equations are well known [1, p. 32] 


= tes wst’s”, (2.12) 
in which g denotes the determinant whose element in row ¢ and column j is g,; 


Equations (2.12) will be used to establish (2.10). Since ¢ is the unit tangent to 7 
| £ 


at P, its contravariant components are defined by 


Oa Or Ox E 1/2 ‘ « 
c=; / (, ae ) 7s teas (2.13) 
On ‘ Or OV 


On substituting these values in (2.12) we find 


a) q l Odaa 
M g'7 2 (£) is” 2.14) 
ox Wu 20 aah Gun) Ox 
\loreover, since do Jun(d.e’)°, we have 
OW Oy dx” 3 oy : (9.15 
° ae Qrn a (2.10) 
do Ox" do Ox 
since g L/gi (i 1,2, +--+, m), Eq. (2.10) is an obvious consequence of (2.14) and 


2.15 
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The components Ap’ of the first curvature vector of 7’ at P are defined by the 
intrinsic derivative of ¢' along 7. Thus, we have 


K yp’ _ oF > a 5 sg (2.16) 


Evaluation of the intrinsic derivative yields 


- : st 1/2 at wh yt 2 
Krp’ = 6 a / Jun = (< s+f¢ ri.) / Jon 
a Ox _ 
2.17) 


6, Og, Re 
9,7 4.7 a : 
29nn OF eS 


In view of (2.8) it follows from (2.17) that Kyp' dp/dxr' is given by 


; ; 0 ‘i , oy aa ; 
K rp ¥ = Kyp’ - Y =¢ TT. re (2.18) 
Ox Ox’ Ox” 


Again making use of (2.8), simple calculations yield the following equations 


ay , ay 1 oq” Ow 
Breas < fe a nee 


a q na _m\2 i oe 
Oa (ox ) 2 ox Ox 
(2.19) 
nn OY (g"")’ 8gnn OW ~( ay r 4 
( = - has | ~4n\2 °+4x£onn a nye 
9 (Ox )~ 2 Ox Ox g (Ox") Ox 


By subtracting the members of (2.18) from the corresponding members of (2.19), Eq. 
(2.11) is obtained. Finally, in view of Eqs. (2.10) and (2.11), which are now established, 


Eq. (2.9) may be written in the form of the relation which constitutes Theorem 1 in 


which p- Vy denotes p’ dy/dz'. 


In the special case in which y = ¢, since p" = 0 and ex” is independent of 
v,a2,-::,2”', the scalars Aid and p' 0¢/dx' vanish. Therefore we have 
0 o 0 , 
ap = <$-M—. (2.20) 
ole Oo 


Corollary 1.1 is an obvious consequence of (2.20). 
In case the system of hypersurfaces ¢ = x” = constant are parallel, x” may be taken 
to be the arc o measured from x” = 0 along a normal. It follows that 


Gm = 1 =", r., = 0, (¢= 1,2, --- ,n). 


In this case the relation between the Laplacians of ¢ becomes 


on fa] 
a (2.21) 
Oo Oo 


Ag a Ad 


In particular, for a V, in a surface, relative to a system of curves geodesically parallel 


[1, p. 57] to V, , (2.21) reduces to 


Od , Ob Od 
Ad = 3 K 5 (2.22) 


feleg Oa Os” 


in which K is the curvature of V, at P and s is are length of V, [3]. 
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3. Proof of Theorem 2. On substituting from (2.8) into (2.17), the components of 
the first curvature vector of 7’ at P are found to be given by 


— ew Aln Q) ' 
Kyp" = 9 —. p =0, (3.1) 
Ox~ 
in which (gan) \’” Q = —d/dc. On forming the inner product of the members of 


(3.1) with gus drg/ds in which dxg/ds denotes the components pz the following relation 


is obtained 


Ky a(In Q)/ds, (3.2) 
in which s denotes the are length of the orthogonal trajectory on V,_, at P of the family 
of isogauges of V . This completes the proof of Theorem 2. 


4. The Laplacian Ad in terms of the supporting function. Let h denote the support- 
ing function of a one-parameter family of surfaces of class C* defined by a functional 
relation of form (1.9). The quantities @, w, @ defined in association with (1.9) will serve 
here as curvilinear coordinates of the point P of any surface. To express the Laplacian 
A@ in terms of the partial derivatives of h with respect to 0, w, ¢, the components of 
the fundamental metric tensors g,; and g'’ must first be calculated. For this purpose we 
make use of the definitions of h, 0, w, @ in terms of a fixed auxiliary rectangular cartesian 
coordinate system (2, ¥, Since (@, w, @) are curvilinear coordinates of the contact 
point P of a tangent plane toa@ = constant surface S, the rectangular cartesian coordi- 
nates (2, y, 2) of P satisfy the following three equations. The first of these equations is 


the equation of the tangent plane at P, and the next two are obtained from it by partial 


differentiation with respect to @ and w, respectively 
rsinw cos 6+ ysinwsin 6+ 2 cosw = h, 
—xrsinwsin 6+ ysinw cos6=h, , (4.1) 
xr cosw cos 6 + y coswsin 6 — zsinw Ra 


in which h. Oh/dw, h, = dh/oe. 
The following set of equations, which has the same solution as that of (4.1) and 
whose coefficients of x, y, z are functions solely of 6, is obtained by combining Eqs. 


(4.1) in the ways indicated by the forms of the right members. 


x cos 6+ ysin 0 hsinw + h. cosa, 


—~xsin 6+ y cos 6 h, esc w, (4.2) 
z = h cosw — hz sinw. 
The solution of (4.2) is given by 
x (hsinw + h, cosw) cos 6 — h, esc wsin 8, 
y = (hsinw + h, cosw) sin 6 + h, ese w cos 8, (4.3) 
z= h cosw — hzsinw. 


On calculating the g’s by making use of the relations 


ax” dx" _ , 
9,= 8 *,j=1,2,3 (4.4) 


Ak o o , 
du’ du’ 








ww 
#8) 
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in which x 1a /, z,uU 6,u° = w, u’ = 4, the following values are obtained 
g ese’ w + f fio = f(g +e ese w), 
q P+ f? ese’ w, fis = —(fhos + ehig ese’ w), . 
1.5) 
Gus ql + fh,, ese w 133 hs + hs. + hh; csc’ a, 
G ] hoH : ese Ww, H’ (eq — fy 
in which e, f, g are defined DV the equations 
h sin’ w h, coswsinw + h,,), 
} h, eot w } g —(h + h, 
the subscripts 1, 2, 3 of A denoting partial derivatives with respect to 6, w, ¢, respectively. 
The functions e, f, g are actually the coefficients of the second fundamental form of S, 
as may be see! DV direct calculations of these coefficients. Let x, xX, X,3 denote the 
vectors whose 7th components are Ox /Ou', Ox /Ou Ou respectively. rom the 
definition of the second fundamental form of S together with some orthogonality rela- 
tions [6, pp. 92-94] it follows that the coefficients of the form are given by the set of 
determinant relations 
XopXiX i] a, B |, 2 1.6 
in which the g’s are defined in (4.5) 
The roliowing expressions for the components of the vectors x x x xX Tue 
are obtained by differentiating Eqs. (4.3 
Tr é S SIN GO - 1 COS w COS GO, 
reese @sin 6 - 7 COB wW COS G, 
» cos 6 f cosw sin 86, 
gecos@sin 6 1 CSC w COS @, 
f sin 
J sin 
e CS i, COS w) Cos 6 + (é sC @ i COS @w sin 4, 
e€&, CSC w f cosw) cos 6 + (e est w f, cos w) sin @, ba 
fF, Sin 
(f ese J, COS w) Cos 6 + (f, esew + g COs w) sin 6, 
a dq Cos f, es@w) cos 8 4 i ¢sCw gd COS w) sin 6, 
x g, Sin 
Z, g sin w Jo COS w) COS 8 + (fz CSC w f ese w cot w) sin 8, 
x f ese w cot w f, ese w) cos 0+ (g sin w J. COS w) sin 8, 
Lo2 = gJ2SINw tr g COS. 
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The three-rowed determinants which appear in (4.6) may now be written in the form 


rig eesewsin @— fcoswcos 6  f csewsin 6 — g COsw Cos 6 
x. ».X,.x ’ cC CSC @ COs 6 — f coswsin 6 —g COS sin @ — f ese w cos 6 |, (4.8) 
Las f sin a) gq sin @ 


in which x4, (2 |, 2, 3) are given by (4.7). 


a 


After transforming this determinant by elementary transformations on the rows 
(the transformations being indicated by the transformed elements of the first column), 


the following form is assumed 
X4gX1X (4.9) 
Vaatrag COL W COS 6+ (xo6 +23 cot w sin 6) tan 6 0 0 
+2, cot wsin 0 —e esc w cos 6 —f est w cos 0 
: : ‘ 
Loe fsnw g sin @ 


(4.7), and (4.9), simple calculations yield the values e, f, g 


In view of Eqs. (4.5), 
d.. of the second fundamental form, respectively. 


for the coefficients d,, , die, 
The components of the contravariant metric tensor, calculated according to the 


definition 
g’ = (cofactor of g;,;)/G, 


are found to have the following values 


(f + q° sin w + (gh, _— fliz,)” hs] , 


g (Cho, — fhys)(ghis — fhos)/h3 — fle + g sin’ w))/H’, (4.10) 
q (ghy, — fh Hh; , gq” = lc ; + f sin’ w + (fh,; _ eho)” h;] rr. 
g (eh, — fh, \/Hh; . gq” = J hs ‘ 


Since the absolute value of d¢/d¢ is equal to the length of the gradient of ¢, where 
the components of the gradient of ¢,are 6; , we have the relations 


(2#) _ (Vo)* =— g 55. — gq: = l he ba (4.11) 


Oo 
For a system of closed surfaces ¢ = constant, according to the convention adopted for 
the sign of o, the direction of increasing o is outward. Therefore, if the origin is an interior 
point of the system, h increases as o increases, and it follows that 
0¢/d0 = 1/hz = —Q. (4.12) 
The components 7’ of the unit normal to the surface S(¢@ = c) at P are given by 


n* = hag”. (4.13) 


From (4.12) and (4.13) it follows that 
a°¢ Ah h, ‘ 
ad 1dh;_ 1 (* ) a) (4.14) 


ac he do ~~ hy Ou’ 


3 
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in which u' = 6, u = w, u’ = ¢. On substituting the values of the g’s as given by (4.10) 


into (4.14) and combining the result with (4.12) to evaluate the expression for A@ given 
by (2.20) one readily obtains the relation 
—h; Ap = Mhs + has + (eh2s — 2fhishes + ghis)/H. (4.15) 
In view of (4.14), Eq. (2.20) may now be written in the form 


Oh. 


do 


h; Ad + Mh, + 0 (4.16) 
which appears in Sec. 1 as (1.11). 

The derivation of (1.10) proceeds as follows. Consider a point P’ on the orthogonal 
trajectory T of S at P whose coordinates are (@ + dé, w + dw, @ + dd). As P’ varies 
along 7’ tending toward P the following limiting relations are approached 

dw/d6 = g/g", do/dw = g°/g", (4.17) 
since the contravariant components of the gradient of ¢ are given by g'*. Let Q denote 
the point of intersection of S with the ¢-curve passing through P’. The coordinates of 
Q are (0 + dé, w + dw, ¢). It is now clear that the limit of the direction of QP as P’ 
and Q tend to P is given by the first equation of (4.17). It is clear from the definitions 
of P’ and Q relative to P that the limit of the direction of QP is the direction at P of 
the line of intersection of the tangent plane to S at P with the normal plane to S which 
contains the tangent to the ¢-curve at P. The normal curvature of S at P which cor- 
responds to this direction will be denoted by K. It follows that K is given by the relation 

K = (e(g’*)’ + 2fg°g” + gg) )/gasg" gy’ - (4.18) 
On substituting from (4.4) into (4.12) we find 
K eho, — 2fhishos + ghis)/H|hos + his ese’ w]. (4.19) 

Let y denote the acute angle between the normal to S at P and the tangent to the 

¢-curve at P. The covariant components ¢, of the unit normal are given by 

t; = 8:/(g")'”. (4.20) 
The contravariant components é° of the unit tangent to the ¢-curve at P are defined by 

&° = 33/(gas)’’- (4.21) 
It follows that cos y is given by the relation 
COS ¥ ¢é° = 1/(9°gss)'”” —h;/(hz + hig + his ese’ w)'”’. (4.22) 
Therefore tan’ y is given by the relation 

tan’ y = (hoz + his ese’ w)/h; . (4.23) 

In view of this and (4.19), relation (4.15) may be written in the form 


—h; Ad = Mh, + Riss +- hiK tan“y (4.24) 


which, in view of (4.16), is equivalent to (1.10). 


If @ is a harmonic function, Eqs. (4.16) and (4.24) assume the forms 
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dhs + Mh, = 0, (4.25) 
00 
haz + (M + K tan’ y)hi = 0, (4.26) 


respectively. 

In view of (4.12), Eq. (4.26) can be written in the form 
dQ 
ag 

Equation (1.14) and Theorem 4 (stated in Sec. 1) now follow immediately. 

5. The equation of Laplace in two variables. Consider a family F of ¢ = constant 
cylinders whose generators are parallel to a fixed axis, say the z-axis. The supporting 
function of a contact point P of a tangent plane to the @¢ = constant cylinder is inde- 
pendent of the position of P on the generator. The only admissible value of w is, therefore, 
7/2 and the number triple (@, 7/2, @) determines a generator rather than a point of 
the generator. In order to apply results of Sec. 4 to this case it will be convenient to 
determine a point P of a generator by specifying its z coordinate. Thus we augment 


+ M+ Ktan’y = 0. (4.27) 


— 


the original curvilinear coordinates 6, @ by the coordinate vu’ = z. 
The auxiliary rectangular coordinates of P satisfy the set of equations 
xecosé+ ysin@=h 
' ' (5.1) 
—xsin 0+ y cos 6 = h, , 
in which h h(0,o), h, = dh/00, h, = dh/dz, h, = dh/dg. The solution (x, y) of (5.1) 
is given by 
x =h cos @ — h, sin @, " 
(5.2) 
y = hsin 6+ h, cos @. 
The g’s of the metric tensor in terms of the coordinates u' = 0, u° = z, u® = ¢, may 
be easily calculated. For this purpose, the following values of the partial derivatives 


are needed. 


Ox , oY Oz 

—— -(h + h,,) sin @, — (h + h,,) cos 86, — 

ou Ou Ou 

Ox Ou Oz = 
==, 0, —; 0, —, = I, (5.3) 
ou ou Ou 

Ox ; . OY , Oz 

—, h, cos @ — hy, sin 6, —; = h,sin 0+ h,; cos 8, —;z = 0 

Ou Ou Ou 


Since the g’s are defined in terms of these partial derivatives by a set of equations of 
form (4.2) in which 2’ = x, x = y, 2° = z, the following values are readily obtained 
Jis (h + hi )his ’ J23 = 0, 933 = hs + his ’ (5.4) 

gi = (h+h,,)’, fiz = O, Jo2 = 1. 
The quantity (h + h,,)° will be shown to be the square of the radius of curvature 


of the curve C of normal section of a ¢ = constant cylinder. For this purpose, let a 








332 P. O. BELL [Vol. XVI, No. 4 


denote the angle of inclination of the tangent to C at P with respect to the x-axis. Then 
a = 6+ 77/2. From the definition of curvature K of C, we have 


Kk = da ds = dé ds, 


where differentiation is with respect to are length of C at P. On differentiating (5.2) 
with respect to s we find 
10 dy _ dé 


(h + h,,) cos 0 . (5.5) 


dx ; ( 
—(h + h,,) sin 6 
ds 


ds ds ’ ds 
Squaring (5.5), and adding yields 


dé 


) = K°*(h + ar: 
ds 


1=(h+ hay 
This is what was to be shown, since r° is-by definition 1/K~. 
It will be shown that the necessary and sufficient condition that @ be a solution of 
the equation of Laplace is that the function h be a solution of the equation 


yp - ° 
= « Dp +d, (5.6) 


0g 
in which r is the radius of curvature of C at P, and p and q are defined by the relations 
oh _ Op 


Pe ae’ 1 00 


If the origin of rectangular coordinates is an interior point of the system of ¢@ = 
constant cylinders, Eq. (4.12) holds. Recalling that the positive sense of o is that in 
which @ decreases, and that the gauge function diminishes as P moves along a flow 


line (o increasing) we find 


dQ ao 0 , 


ae “ (hs) <0, @>O0. (5.7) 
Oo Oo Oo 
A relation of form (4.14) holds for a system of ¢ = constant surfaces, independently of 
the choice of the coordinates u', u’, and consequently, in the present coordinates u' = 86, 
u” = z. Thus 
29 _ dhs 13 3 - 
—h, 3 = = hs(hsig”” + hssg”). (9.3) 
Oo Oo ‘ : 
But in view of (4.16), dh;/dc0 may be replaced by —Mh, if and only if ¢ is harmonic. 
On making this replacement and dividing by A; , Eq. (5.8) assumes the form 
re " 7 
—hy SS = —M = hag” + hag”. (5.9) 
Oo 


From (5.7), M is found to be negative. The mean curvature M which is the sum of the 
principal curvatures of the cylinder @ = constant at P, is given by the relation 


M = -lI/,, (5.10) 


since one of the principal curvatures is the curvature of C at P and the other is zero, 
the curvature of the generator at P. 
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1959] 
The contravariant components g’’ are defined by the relations 

g'' = (cofactor of g;;)/| gi; |, t,j = 1, 2,3. 
In view of (5.4) the determinant of the g’s is given by the equation 
r 0 this 
Ji; = | 0 1 0 
lrhig O hi + hi; 





It follows that g'*, g** are given by the relations 


13 


g° = —hy rhs , g =! ne . (5.11) 
On substituting from (5.10) and (5.11) into (5.9) and clearing, the resulting equation is 
rhs, = (hi + hi,), (5.12) 


which is (5.6), since h; = 0h/0% = p, his = Op/00, hs, = Ap/dd. 
The theorems stated in the introduction can be easily specialized for the two dimen- 
sional case. The statements of these two dimensional theorems will be omitted here. 
Additional results which are not specializations of the theorems hereto proved, but 
whose proofs involve only routine calculations will be presented here (without proof). 
Theorem 5.1. The ratio of the gauge functions at points P, , P of a d-curve is given by 
the integral 


2,/2 = exp [ (K? + K2)"? ds (5.13) 
in which the integration is with respect to arc length of the ¢-curve from P, to P, and K, Ky 
are the curvatures at P of the @ = constant curve (isothermal curve) at P and the orthogonal 
trajectory of the @-curves at P, respectively. 

Theorem 5.2. If the density of charge at each point P of the curve C, of a clossd conductor 
is proportional to a constant power of the supporting function h of P, then the exponent of 
the power is unity and the curve Cy is an ellipse, and the equipotential curves induced by 
the distribution of charge on C,, are ellipses confocal with C, . 
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BOOK REVIEWS 


Ndherungsmethoden der héheren Analysis. By L. W. Kantorowitsch and W. I. Krylow. 
VEB Deutscher Verlag der Wissenschaften, Berlin, 1956. xi + 611 pp. $11.25. 


This book is a German translation from the third Russian edition (1950), the latter being only a 
minor revision of the second Russian edition (1941). It presents a wealth of methods for the approxi- 
mate solution of partial differential equations, the main interest being centered about elliptic equations, 
especially Laplace’s and Poisson’s equations, and the biharmonic equation, with the usual boundary 
conditions. In Chapter I methods based on the representation of the solution as an infinite series are 
treated with emphasis on Fourier series and orthogonal functions. A section of this chapter is devoted 
to infinite systems of linear equations. Chapter II contains methods for the approximate solution of the 
Fredholm integral equation together with an application to the problems of Dirichlet and Neumann. 
In Chapter III the differential equation is approximated by a corresponding difference equation and 
methods for solving the latter are studied. Chapter IV is one of the most interesting of the book and is 
devoted to variational methods. There is given a clear account of the methods of Ritz and Galerkin; 
in contrast to many other presentations of this subject the question of convergence is carefully examined 
and estimations of error are given for some simple but important cases. The remainder of the book, 
Chapters V-VII, deals with various questions concerning conformal mapping, application of conformal 
mapping to boundary problems for harmonic functions, and the alternating method of Schwarz. 

The book presents an excellent survey of different analytical tools of importance in numerical 
analysis, the level being somewhat more advanced than is ordinarily the case in text-books on this 
subject. Theorems are proved in a stringent way and much emphasis is laid upon estimations of error. 
Several theorems can, however, be shown to be true under more general conditions than are stated here, 
but the generality of the assumptions is entirely sufficient for the applied mathematician. Every new 
method is illustrated by one or two well-chosen problems, which are worked out in detail. There are 
now exercises for the reader. 

Since different methods mentioned in the book ean often be applied to the same problem, it would 
have been valuable if the advantages and disadvantages of the respective methods had been discussed. 
As is usual in Russian mathematical text-books references are made almost exclusively to Russian 
papers. Since the book except for minor changes is the same as the 1941 edition, the development of 
numerical analysis during the last fifteen years has not influenced the contents, for instance the unified 
treatment of many methods and theorems which can be made by means of functional analysis is not 
used and the possibility of applying a computer is not discussed. These minor objections are, however, 
in no sense important. As a whole, the book is a very interesting and useful contribution to the liter- 


ature on numerical analysis 


Experimental designs. By William G. Cochran and Gertrude M. Cox. John Wiley & Sons, 
Inc., New York, and Chapman & Hall, Ltd., London, 1957. xiv + 611 pp. $10.25. 


This book is ‘‘intended to serve as a handbook which is consulted when a new experiment is under 
consideration.” The first edition was highly successful in this aim and the second edition is an augmenta- 
tion of the first edition to take care of advances in the design of comparative experiments in the last 
few years. 

This is the best book in existence on the down-to-earth practice of experimental design. There are 
aspects of the second edition which are unsatisfactory and these arise because the material of the first 
edition was transferred to the second edition in toto. Some of the material of the first edition could have 
been improved and the continuity of old and new material improved. 

The book will be of very high value to experimenters and experimental statisticians and one can 
also surmise that workers in mathematical statistics will obtain some benefit in knowledge and outlook 
by perusing the book, since statistics is a subject which has in the past gained so much from under- 
standing the needs of research workers. 

O. KEMPTHORNE 


(Continued on p. 396) 
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MULTIPLE FOURIER ANALYSIS IN RECTIFIER PROBLEMS* 
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Laboratory for Electronics, Inc., Boston, Massachusetts 


ROBERT L 


Abstract. The non-linear problem of the multiple Fourier analysis of the output 
from a cut-off power law rectifier responding to a several-frequency input is reviewed 
for the one- and two-frequency problems and is briefly investigated for the three-frequency 
problem. The solutions for the modulation product amplitudes or multiple Fourier 
coefficients are obtained in exact although transcendental form. An account of the mathe- 
matical properties of these multiple Fourier coefficients or Bennett functions, including 
hypergeometric representations and power series expansions for them as well as recur- 
rence relations satisfied by them, is given in the paper together with line graphs of the 
first ten basic functions for the one-frequency problem and of the first fifteen basic 
functions for the two-frequency problem. Further applications of the theory are also 
given to the computation of average output power with the aid of the multiple Fourier 
coefficients or Bennett functions studied in the paper, and the work is concluded with 
some brief remarks concerning the interpretation of the results in terms of the theory 
of almost periodic functions and the generalized Fourier series of Bohr under appropriate 
conditions. Numerical tables of the functions graphed have been prepared and are 
ivailable separately in the United States and Great Britain for applications requiring 
great accuracy. Finally, the entire theory is based on the original method of the expansion 
of the rectifier output in multiple Fourier series introduced by Bennett in 1933 and 1947. 

1. Introduction. A non-linear problem of continuing interest in theoretical elec- 
ronics, and one much studied in the one- and two-frequency instances, is the multiple 
lourier analysis of the output from a cut-off power law rectifier responding to a several- 
frequency input under various assumptions about the cut-off bias, about the exponent 
of the power law, and about the input amplitude ratios. Among others Bennett [1, 2], 
Kaufman [3, 4], Lampard [5], Salzberg [6], and Sternberg, Kaufman, Shipman, and 
Thurston [7, 8, 9, 10, 11] have obtained exact although transcendental solutions for the 
modulation product amplitudes or multiple Fourier coefficients in the first two cases 
mentioned and have tabulated some of these quantities under the name of Bennett 
functions, while Feuerstein [12], following a method of Bennett [1] and Rice [13], has 
recently given a quite extensive treatment of the n-frequency problem in terms of gen- 
eralized Weber-Schafheitlin integrals [14] suitable for numerical evaluation. 

In this paper, after reviewing some of the more interesting results noted above for 
the one- and two-frequency problems associated with a cut-off power law rectifier, 
with one or two new results added, we briefly investigate the three-frequency problem 
by a transcendental method of extension of previous results using, particularly, certain 
integral recurrence relations which hold between the multiple Fourier coefficients or 
Bennett functions of different multiplicities. Since the frequencies and phase angles are 
readily found, the real problem with which we are concerned in each instance consists 
of evaluating the single, double, or triple integrals defining the Fourier coefficients in 
the multiple Fourier series expansion of the output from the rectifier. The mathematical 
properties of these multiple Fourier coefficients or Bennett functions are investigated 
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systematically, and hypergeometric representations and power series expansions for 
them and recurrence relations satisfied by them are given in the paper, together with 
line graphs of the first ten basic functions for the one-frequency problem and of the first 
fifteen basic functions for the two-frequency problem. In addition to providing for the 
multiple Fourier series expansion of the output from the rectifier, it is shown that the 
average output power may also be computed in terms of the same multiple Fourier 
coefficients or Bennett functions as those studied, provided, in the two- and three- 
frequency cases, that the input is non-periodic. No attempt is made to discuss modulation 
product problems more general than those associated with rectifiers, the reader in such 
instances being referred to Sternberg and Kaufman |7, 8, 9] for a general approximation 
process applicable to any continuous modulator, the solution again being formulated 
in terms of Bennett functions as here. The paper is concluded with some brief remarks 
concerning the interpretation of the results in terms of the theory of almost periodic 
functions and the generalized Fourier series of Bohr, again under the proviso, in the 
two- and three-frequency cases, that the input be non-periodic. Numerical tables giving 
units at values of the arguments 


values of the functions graphed good to 1 X 10 
spaced one-tenth of a unit apart have been prepared and are available separately from 
the sources noted below in the United States and Great Britain for applications requiring 
great accuracy. Finally, we may note in passing that the entire theory is based on the 
original method of the expansion of the rectifier output in multiple Fourier series intro- 
duced by Bennett [1, 2] in 1933 and 1947. 

The general rectifier problem is formulated precisely and the solution in simple or 
multiple Fourier series is outlined in Sec. 2. The review of results for the one- and two- 
frequency problems is presented in Secs. 3 and 4 while the three-frequency problem is 
treated in Sec. 5. The results of the theory which apply to the computation of average 
output power for a single frequency periodic input or a multiple frequency non-periodic 
input are given in Sec. 6, and the concluding remarks bearing on the connections with 
the theory of almost periodic functions and the generalized lourier series of Bohr are 
given in Sec. 7. Copies of the numerical tables of the Bennett functions graphed in the 
paper and available to the public have been deposited with Mr. D. H. Lehmer of the 
Unpublished Mathematical Tables Vile in the United States and with Mr. Andrew 
Young of the Grace Library of the University of Liverpool in Great Britain and with 
the authors. 

2. Formulation of the problem. Consider a cut-off power law rectifier having an 


output versus input characteristic Y = Y’(X; X,) of the form 
. ' ; cee Ae) ; X Res 
y(X;X,) =4 ? . > Me vy > 0. (2.1) 
lo t<x. 


Let the input to the rectifier be a one-, two-, or three-frequency function of time, 


z(t), of one of the forms 


x(t) = P cos (pt + 86,), Pe > ©. (2.2) 
a(t) = P cos (pt + 6,) + Q cos (qt + 8,), Pa G> 0. (2.3) 


or 


a(t) = P cos (pt + 6,) + Q cos (qt + 6,) + R cos (ri + 86,), P>Q2>R>0. (2.4) 
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The output from the rectifier, y(t) = Y"[x(t); Xo] has then a single, double, or triple 
Fourier series expansion of the corresponding form 


o 


y(t) = 4P"AL?(hA) + P” YS AL (h) cos ont + bn); (2.5) 
y(t) = 3P"A(h, BE +P” >O* ALLA, b) 008 emnl + Gamn)s (2.6) 
Ol 


y(t) 2P” Aoool(h, ki , ke) + P” >* AL mnt(h, ky , hq) COB (Woaemnrl + Gaomnt); (2.7) 


000 
mn,l=0 


where hh = X,/P,k = Q/P > 0, and k, = Q/P > k, = R/P > 0, and where the modu- 
lation product angular frequencies, phase angles, and amplitudes are given by the 


formulas 


WO, = mp, o. = m0, , (2.8) 
Wemn = mp + nq, damn = m0, + nO, , (2.9) 
Ween = mp tng + lr, Gaamnt = MO, + nO, + 16, , (2.10) 


and the formulas 
Aa'(h) = : (cos u — h)” cos mu du, 
T J, (2.11) 
R, :cosu > h, O<u<e, 


A he. EN 2 I (cos u + k cosv — h)” cos mu du cos nv dv, 
mee (2.12) 


Ka 


t, :ecosu+keosv > h, 0< u,v <7, 


p> « 


9 . 
A celh, e 5 Ke) =; II (cosu + k, cosv + ks cos w — h)’ 
w Je 


™ (2.13) 


-cos mu du cos nv dv cos lw dw, 
®R; :cosu + k, cosv + k, cos w > h, 0O<u,v,w se, 


where in all of these formulas the indices m, n, l take all integral values m, n, 1 > 0 
and, finally, where the asterisks on the summation signs in (2.6) and (2.7) indicate that 
in these multiple Fourier series we sum only on all distinct arrangements of plus and 
minus signs, equivalent arrangements being taken only with the plus signs and with 
the zero order terms, particularly, having been removed from the sums. 

We will assume the genesis of these Fourier expansions to be well-known and will 
refer the reader to Bennett [1] or Sternberg and Kaufman [7] for additional details of 
such matters. It may be worth noting, however, that the existence of these Fourier 
series follows at once from the continuity of the kernel functions in the integrands of 
the quantities A‘? (h), AL&?.(h, k), and A{”,,,,(h, ki , k2) within the corresponding regions 


=mn 


®, , R, , and ®, while the convergence of these Fourier series to the functions y(t) follows 
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at least at all points of continuity of the y(t), in the case of (2.5) from the elementary 
theory, in case of (2.6) from the theory of the double Fourier series given by Hobson [15] 
or Tonelli [16] and, finally, in the case of (2.7) may presumably be surmised from Hobson’s 
remarks concerning comparative properties of the double and multiple Fourier series in 
general. 

As noted in the introduction the whole problem of multiple harmonic analysis with 
which we are concerned for rectifiers consists of the evaluation of the multiple Fourier 
coefficients or the integrals A‘”(h), AL”.(h, k), and AL”,.,,(h, hk: , ke) in (2.11), (2.12), 
and (2.13) as functions of their several parameters. Aside from certain relatively simple 
closed form solutions obtainable in the one-frequency problem none of these integrals 
can be evaluated in closed form in terms of elementary functions, except for specialized 
values of the parameters h, k, k, , and k, , but rather in general the best that can be 
obtained are solutions in a transcendental form of one type or another with the final 
computations to be done by numerical methods. Although one might wish for more 
simple results in the two- and three-frequency problems, the power of the methods to 
be presented should not be underestimated for either theoretical or practical purposes. 
Finally, in recognition of the original work of Bennett [1, 2] and following previous 
usage, we term the functions A‘ (h), AL? (h, k), and AL”?,.,, (h, k, , k2) Bennett functions 
of the vth kind and of multiplicities one, two, and three respectively. 

3. Review of results for the one-frequency problem. In the one-frequency problem 
for the rectifier (2.1) we have for the input the signal (2.2) and for the output the Fourier 
series (2.5) and have to determine particularly the Fourier coefficients 


An? (h) = 2/ (cos u — h)” cos mu du, (3.1) 
GR; o. 


R, :cosu > h, 0O<u<z_, 


as functions of h = X,/P for fixed v > 0. The variable h takes all real values and, following 
Kaufman [3], we note three cases of the functions A‘’’(h) according to the scheme: 
(0) h > 1, (a) |h| < 1, and (7) h < — 1. Clearly case (0) may be disposed of at once for 
in this case the rectifier is biased so strongly that we have A{”’(h) = 0 for all m. 

Consider now cases (a) and () for all »y > 0. To begin with, the functions A{” (h) 
satisfy a number of recurrence relations. Thus, for determining the (v + 1)th kind 
functions in terms of those of the vth kind we have the recurrence formulas 


4A,’*? = 4A,” — $hA,” 
atin = ere (3.2) 
2mAn” = + DAL, — @+ DAN, 
while for purposes of evaluating the (m + 1)th order functions in terms of those of the 


mth and lower orders we have the recurrence formula 

(m+v+ IAS), = 2mhAL? — (mM —v — 1)AL?,, (3.3) 
where in each of these relations m > 1. For proofs we refer to Kaufman [3]. For actual 
evaluation of the functions A‘ (h) two general solutions are available. In particular, 
in case (a) Bennett [2] has given the formula 


(7/2) — 2\'" Tw +1) ms y+i/2p71 ee , -. 2 a ‘ 
A‘(h) = (2) re tay WF + ma mie + HR - WL, BA) 
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while in case (©) Lampard [5] has given the result 


eS ee 2 = 1 

An’(h) = > h Fi[3(m — v), 3m —v+1);m4+1;h"J, (3.5) 
m! Tv — m + 1) 

where in each of these formulas F(a, b; c; x) denotes the Gaussian hypergeometric function 

and where m > 0. Various alternative forms of these solutions have also been noted by 

the same authors, Bennett, in particular, indicating certain special results for integral 

and half integral values of v; see also Salzberg [6]. 

Consider next cases (a) and () for integral values of v. In view of the recurrence 
relations (3.2) for determining the functions A‘’*’(h) in terms of the functions A{&’(h), 
it is clear that for such values of v the solutions may be completely expressed in terms of 
those for a basic set such as A‘ (h) or A,,(h) = A{(h). Actually we choose the latter 
as being more important in their own right. In case (a) following Bennett [2] we may 
simplify (3.4) or integrate directly with Salzberg [6] to obtain the special results 


4 A,(h) _ : (e a. h?)'/? es h cos”: hj, 


A,(h) = : [cos ' h — h(1 — h’)'”], (3.6) 
T 


II 


A,,(} 
ACh) m— 1 m+ 1 


mr 


| {ei lim — 0s ‘h) _ sin [(m + 1) cos Alt 


where m > 2, while in case (~) we have 


SA (h) = | hi], A,(h) = 1, A,,(h) = 0, (3.7) 
where again m > 2. These formulas represent an essentially complete solution for the 
functions A,(h) = A?(h). For many purposes, however, the corresponding power 


series expansions are more illuminating. Expanding the elementary functions in (3.6) 
and combining terms we have, for example, in case (a) for the first few functions A,,(h) = 


A‘ (h) the series 


1 aa l ee... l pe L 4 i 6 + 
5 A,(} ~ l . h+ . he + D4 h* + 30 h+ . (3.8) 
1,(h) zm tig dre bare ---| (3.9) 
ee Ta” ae 56 , i 
SIS. ws bees 3 iw ss... 
A,(h) 1/2 he + ra + 35 hi + | (3.10) 
= aoe | 
A,(h) = 713 h +7 h + 94 + ; (3.11) 
8 2 we. Sea Pw :.. 
A,(h) = - 15 +h 4 h* + 4 hi + | (3.12) 
Mm @ck— te 8a | 
Adi} = =] W524 =H = re ate 
1 (h) =t | 5 tah ra gh , (3.13) 





each of which converges uniformly for | h | < 1 and, indeed quite rapidly, the terms 
given above, for example, being sufficient to evaluate the functions }A,(h) to A;(h) to 
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better than 1 X 10° units for | h| < 34. Finally, in addition to the above formulas and 
series for the functions A,,(h) = A‘'’(h) in cases (a) and (@) Kaufman [3] has given 


the reflection relations 
2A(—h) = 4A(hk) +h, Ai(—h) = 1 — A,(h), A,(—h) = (-—1)"A,(h), = (3.14) 


for determining the functions of negative bias in terms of those of positive bias where 
again m > 2. 

For elementary applications and general interest line graphs of the first ten functions 
bAo(h) = $A,’ (h) and A,,(h) = A®(h) are presented in Fig. 1 while for applications 
of greater accuracy the reader may refer to the tables of these functions mentioned in 


the introduction. 





06 




















1 


Fic. 1. The functions 4 Ao(h) and An(h, k) 


4. Review of results for the two-frequency problem. In the two-frequency problem 


for the rectifier (2.1) we have for the input the signal (2.3) and for the output the double 
I £ I 


Fourier series (2.6) and have to determine particularly the Fourier coefficients 


9 rp 
tk, jes [| (cosu + k cosv — h)” cos mu du cos nv av, 
wT Je 
Rs (4.1) 
R, :cosu+ kcosv>h, 0<u,v <7, 


as functions of h X,/P and k = Q/P for fixed v > 0. The variable h takes all real 
values while the variable k ranges over the interval 09 < k < 1 and, as before, we note 


three cases of the functions A”, (h, k) according as: (0)h > 1+ hk, (a)i\h| <1 +k, 
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or (2) h < —1 — k. Clearly again case (0) may be disposed of at once, the rectifier 
being biased so strongly that we have A{’?,(h, k) = 0 for all m, n while in all cases the 


funetions A ‘”. (h, k) are independent of the double sign so that we may drop it hereafter 


and write simply A‘%(h, k). 

Consider now cases (a) and (@) for all vy > 0. As before the functions A‘? (h, k) 
satisfy a number of recurrence relations. Thus, first of all comparing the definitions 
(3.1) and (4.1) of the functions A{’’(h) and A‘,’(h, k) one readily observes the funda- 


mental formulas 
ee 
t’(h, &) | A, (h — k cos v) cos nv dv, 
T et) 
(4.2) 


n 


ge” 7 ' 
A’ Ah, k) | A, (h’ — k’ cos u) cos mu du, 
0 


where h’ h/k and k’ = 1/k, which express the A‘?)(h, k) in terms of the A{”’(A) and, 
hence, in a sense constitute a complete solution of the two-frequency problem in terms 
of the solutions of the one-frequency problem. Next, for determining the (vy + 1)th 
kind of functions in terms of those of the vth kind we have much as before the recurrence 


formulas 


(a) 2Ann) = 3Ayo + 3KAni — $hRAoo 

(b QnA." (y+ 1)Av,, —w+t DAL, (4.3) 
(c nA’? = (vw + IKALY’, — w+ DEAL: 

where in (b) m > 1, n > O, and in (c) m > 0, n > 1, while for purposes of evaluating 


the (m + n + 1)th order functions in terms of those of the (m + n)th and lower orders 


we have the recurrence formulas 


a (m+tn+rv4+ )1DA 2mhA \, 
9 ’ (») 
— 2mkA,,-) — (mM —n—v — DAg ia, 
(b m—nt+v+3)Anion-) 2m + 1IhAS?, a- 


— 20m+ 1I)kA,, — (mtn—v— DAL, 14) 
(4. 


(ec n+m+vt 1Av iss = Q2nh’A,, 
— 2Qnk’A’,, — (n— m—v— IA, 4- 
s— 3 m > 
(d (n mtv+t+ 3)An i nse = An + IN’An i ner 
, ¥ (yr) 
— An+ 1I)k’AS)., —mtm—v— DAL a», 
where h’ h/k and k’ 1 /k and where in (a) and (ec) m > 1, n > 1, in (b) m > 0, 
m 1, and in (d) m > 1, > 0. The proofs of these relations, though computationally 


lengthy, are entirely elementary and may be effected by combined use of the recurrence 
relations (3.2) and (3.3) for the functions A{”’(h) and the formulas (4.2) for expressing 


the funetions A‘?(h, k) in terms of the functions A&?(h); see also Kaufman [4]. For 
actual evaluation of the functions A ‘?(h, k) only one general solution seems to be avail- 


able. Thus, in case (@) Lampard [5] has given the formula 
ae rv + 1) 
Aa, 0 ee ee eee A 
min! Tw—- m—n+ 1) (4.5) 


-F [km +n —v, (mtn —v4+1);m+1,n4+1;h?, Ph), 


"4l2 
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where F’,(a, b; c, d; x, y) denotes Appell’s fourth type of hypergeometric function of 
two variables [17] and where m > 0, n > 0. Other forms of this solution have also been 
noted by the same writer. 

Consider next cases (a) and (~) for integral values of v. Similarly as before the 
recurrence relations (4.3) for determining the functions A$%’*")(h, k) in terms of the 
functions A‘”(h, k) make it clear that for such values of v the solutions may be com- 
pletely expressed in terms of those for a basic set such as A‘,(h, k) or A,,,(h, k) = 
A‘ (h, k) and again we make the latter choice. Even for these functions it is not possible 
to give a complete solution free from integral signs or numerical procedures in case (a). 
However, in the subcase of case (a) in which | h | + k < 1, using a method amounting 
to substitution of the series (3.8) to (3.11) and so on for the functions A,,(h) = AS (h) 
into the formulas (4.2) for the functions A‘”?(h, k), Sternberg and Kaufman [7] have 
derived double power series expansions of the forms 








| l 2 * ‘ a. 6 sir 
; + 5h t+ 54 / + 9 oo 
| I 2 l 27.2 3 47.2 
pty hk +ghk +a5hk +--- 
. Aw (h, k) = as h+ a , (4.6) 
2 2 “lates 9 bh a. 
| 64 ~ 128 
| ak I ke+. 
| 256 | 
2+ Cee ee hv + - | 
3 20 56 
| | | 
‘ee +ih k? + s. h'k? + | 
1 oh ° | 
A,(h,k) = =+-4 (4.7) 
2 T ) 3 4 15 27.4 ( 
| +39" + 64 | 
} = | 
128 | t | 
1 3 sa e: > se = 
at io” ~ 112" | 
oe «2a ~ Say — 
k hk 8 16 64 
Ao (h, k) = — : 38 
( ) 5 4 . 4 3 ” 25 sie \ (4.8) 
4° #13 
| 25 76 
1024 
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! co. 2 
i= 20 = 20 = 2 = | 
> - 9 
2p — 2 yy — Sig | 
re | 8 16 64 
A, (h, k) = —? 7 ” (4.13) 
A ee 
64 128 | 
_ 35 
1024 ~ 
a ae SES 
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ftw. £ Se 
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A.Ah, kb) = —3 - _ (4.14) 
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rT | 15 is! 1445 12h! 
| 1024 12288 * 
245 
~ Q4576 7 


and so forth, each of which converges uniformly for | h | + k < 1 and, in fact, quite 
rapidly, the terms given above, for example, being sufficient to evaluate the functions 
1 Aoo(h, k) to Aos(h, k) to better than 1 X 10° units for | h | + k < }. In case (~) 


as before we have 


5 Agolh, k) =|h x Ajo(h, k) = 1, 


; (4.16) 

Ao(h, k) = k, Aah, k) = 0, 
where m + n > 2. For half-integral values of » in the sub-case of case (a) in which 
h | + k < 1 similar double power series expansions for the functions A‘(h, k) also 


may be obtained by combined application of the formulas (3.4) for the functions A {”’ (h) 
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v. 


and (4.2) for the functions A‘(h, k), the half-integral property of »v making possible 
the explicit computation of the coefficients in the series, while in case (~) the formula 
(4.5) for the functions A‘Z)(h, k) also applies with some slight simplifications of the 
gamma functions. Finally, in addition to the above formulas and series for the functions 
Ann(h, k) = A“? (h, k) in cases (a) and (~) Bennett [2] and Sternberg and Kaufman [7] 


have given the reflection relations 
4 Aoo(—h, k) FAoolh, k) +h, A,(—h, k) = 1 — A, (h, k), £17) 
\ eid 
Aoi(—h, k) k — Ao,(h, k), Amni(—h, k) = (-1)"""A,a(h, k), 


for determining the functions of negative bias in terms of those of positive bias where 
the aid of previous results for the functions A,,(h) = AS (h) and the formulas (4.2) 
for the functions A‘”)(h, k). 

For elementary applications and general usefulness line graphs of the first fifteen 
functions 4A (h, k) = 4A$o(h, k) and A,,,(h, k) = A‘2(h, k) are presented in Figs. 
2 to 16 while for applications requiring greater accuracy the reader as before may refer 
to the tables of these functions described in the introduction. 

5. Some new results for the three-frequency problem. In the three-frequency 
problem for the rectifier (2.1) we have for the input the signal (2.4) and for the output 
the triple Fourier series (2.7) and have now to determine the Fourier coefficients 


again m + n > 2. These formulas also may be derived in an elementary manner with 


*<mnl 


9 e 
Ave nt(h, ky , ke) = 3 iff (cosu + k, cosv + k, cos w — h)’ 
WT JJ 


ss (5.1) 


-cos mu du cos nv dv cos lw dw, 


R®; :cosu + k, cosv + k, cosw > h, 0<4u,9,w <7, 
as functions of h = X,/P, k, = Q/P, and k, = R/P for fixed v > 0. The variable h 
again takes all real values while the variables k, and k, range over the intervals 0 < k, < 
k, < 1 and, as before, there are three cases of the functions A 2”2,,,,,(h, k, , k2) according as: 


OAPI +k,+k,,(a)|h| Sl +h +hkh.,or(e)h < —1 —k, — k, . Similarly 
as before in case (0) we have A2”,,,,(h, k, , k2) = 0 for all m, n, l while in all cases the 


functions A” ,.,(h, k,; , k.) are independent of the double signs so that we may write 


simply A?) (h, k, , k 

Consider now cases (a) and () for all vy > 0. As before the functions A”) (h, k, , k2) 
satisfy a number of recurrence relations. Thus, to begin with comparing the definitions 
(4.1) and (5.1) of the functions A‘? (h, k) and A,?)(h, k, , k.), one readily observes in a 


manner similar to that previously used the fundamental formulas 


as : 
Avth, k, ,k) ==" | Ay(h — kz cos w, k,) cos lw dw, 
T Jo 
' = ' ss 
Ag” th, k; , ke) Ayi(h — k, cosv, kz) cos nv dv, (52) 
a di 
ky wg P)/Le / / 7 
i(k, hy 5 he) = Ant (hi — ki cos u, kik.) cos mu du, 
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where hi = h/k, and kf = 1/k, , which express the A,{y) (A, k, , k,) in terms of the A  (h, k) 
and, hence, in a sense constitute a complete solution of the three-frequency problem in 
terms of the solutions of the two-frequency problem. Next, for determining the higher 
kind and higher order functions in terms of those of the vth kind and (m + n + lI)th 


order we have respectively the recurrence formulas 


(a OY heed = 4 A too + 1k, Aoto + the. loos = bh Asoo ’ 
(b) 2mAva? = V+ NAR iat — OF DAwsiint » (5.3) 
(ec QnAsnt (vv + Dk Anni — & + DkiAnnsi.t » 
(d) y Pala (vy + 1)k,A 1 — + DAR as 
where in (b) m > 1,n > 0,1 > 0, in (c) m > 0,n > 1,1 > O, and in (d) m > 0, n 2 0, 


1 > 1, and the recurrence formulas 


(a mt+nt+Il+p 4 DAS a act QmhA,,”; 
— 2mk, At 1.1. — 2mMk.AP wt 
—(m—n—l—v— 1lAx iat, 
(b m n+t+l+v+t+ 3)Anien-1 = 2m + IhAL neve 
— Hm + DEAS. «.t — Am + Yk Aansaie~1.t-1 


—(m+n—l—v— lAna-i.t; (5.4) 


(¢c mtn — ltot 3Agtent-1) = Am + YhAgir.s.t-1 


— Am + 1k. Avi n.e — Wm + Dk Anes n-t.t-1 


—(m—-n+l—v — DALn 0K; 
(d m r l+yr4 prs 1 = Am + WhAgion-1.1-1 
— Am + 2)ko n-t.t — 2m + Wk Ansent-1 
—(m+nt l —»v — DAgir.2-1.0-1 > 
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(c n+ m+p-e 3)A,, 5 sae = Wn + DAL AS. 
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= 1/k, , and ki = 1/k, and where in formulas (5.4)-(a), 


n > 1,1 > 1; in formulas (5.4)-(b) and (c) m > 0, 
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> 1,1 > 1; in formulas (5.5)-(b) and (c) m > 1, n > 0,1 > 1; in formulas (5.6)-(b) 
and (c) m > 1,” > 1,1 > O while in formulas (5.4)-(d), (5.5)-(d), and (5.6)-(d) we 
have respectively m > —l,n >1,l1>1;m>1,n > —-1,l > ljandm>1I1,n > 1, 
| > —1. As before the proof of these relations, though computationally lengthy, is 
entirely elementary and may be carried out by combined application of the recurrence 
relations (4.3) and (4.4) for the functions A ‘%’(h, k) and the formulas (5.2) for expressing 
the functions A ,{’)(h, k, , kz) in terms of the functions A‘ (h, k). For actual evaluation ot 
the functions A,.")(h, k, , k2) no general solutions free from integral signs or numerical 
procedures appear to be known, the formulas (5.2), the Sternberg and Kaufman approxi- 
mation process [7, 8, 9] and Feuerstein’s generalized Weber-Schafheitlin integral repre- 
sentations [12] being the most general and useful points of departure for numerical work. 

Consider next cases (a) and () for integral values of v and for similar reasons as 
previously, consider particularly, the basic set of functions Aya:(h, ki , ke) = 
A!) (h, ky , ke). Clearly again it is not possible even here to give a complete solution free 
from integral signs or numerical procedures in case (a). However, in the subcase of case 
a) in which | h| + k, + k, < 1, using methods similar to those used before of substitution 
of the series (4.6) to (4.11) and so on for the functions A,,,(h, k) = A‘%2(h, k) into the 
formulas (5.2) for the functions A ,“”) (h, k, , k2), we may again obtain a solution in multiple 
in particular, for the case of half-wave rectification in which the bias 


n= ft, 


power series; 


} 0 we have double’ power series expansions of the forms 
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and so forth, each of which converges uniformly for k, + k, < 1, the k, and k, being 
positive and, in fact as before, the convergence is quite rapid, the terms given above, 
for example, being sufficient to evaluate the functions 3A o90(0, k, , kz) to Ayoo(0, ky, , ke) 
to better than 1 X 10°° units for k, + kz < 3. In case () as before we have 


 Aooolh, k, cs = | h I, Ajyoolh, k, ’ kz) — 1, Aoiolh, k, , k2) - k, ’ 


(5.15) 
Aoolh, ky , ke) = ke , Amnulh, k, ,k.) = 0, 
where m + n + / > 2. As before for half-integral values of v in the sub-case of case (a) 


in which | hk > + k, + k, < 1 comparable double or triple power series expansions for 
the functions A") (0, ky , ko) or ALPi(h, k, , k2) also may be obtained by combined use 
of the formulas (3.4) for the functions A{”’ (h), formulas (4.2) for the functions A‘? (h, k), 
and formulas (5.2) for the functions A") (h, k; , k2), the half-integral character of »v 
again making possible the explicit computation of the coefficients in the series but in 
case (~) no special results exist. Finally, in addition to the above formulas and series 
for the functions A,,,.(h, ky , k2) = AS? (h, ky , ke) in cases (2) and (~) we may note 


the reflection relations 
§ Aooo(—h, ky , ke) = 4Aooolh, ky , ke) +h, 
Ajool—h, ky , ke) = 1 — Ayjool(h, ky, , ke), 


Aool—h, ky , ke) =k, — Aorwlh, k, , ke), (5.16) 

Aon(—h, ky , ke) = ke — Aoolh, k, , k:), 

Awat(—h, ky , ke) = (—1)"*"** A. th, ke, , he), 
for determining the functions of negative bias in terms of those of positive bias where 
again m + n + 1 > 2. These formulas likewise may be established in an elementary 
manner by the use of previous results for the functions A,,,,(h, k) = A<)’(h, k) and the 


formulas (5.2) for the functions A,‘’)(h, k, , k). 
6. Applications of the theory to the computation of power. ‘The theory of the 


preceding sections concerning the functions A‘ (h), AS (h, k), and A,?)(h, ky , k2) and 
the graphs of the functions A,,(h) = A‘?(h) and A,,,(h, k) = A‘2(h, k) given in this 
paper not only may be applied to the computation of modulation product amplitudes 
but also to the computation of the average output power of the rectifier (2.1) when 
responding to a several-frequency input of the form (2.2), (2.3), or (2.4), provided, in 
the two- and three-frequency problems, that the input is non-periodic. Thus, applying 
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the Parseval theorem [15, 16] or Bessel equality for simple or multiple Fourier series 
and assuming for the moment that the average output power @, in each case may be 
taken to be one-half of the sum of the squares of the corresponding Fourier coefficients, 


we have in the one-, two-, and three-frequency problems respectively the formulas 


> = P*?’Ay”*(h) + 5pe p> Ay”’*(h) = = PY AS (h), (6.1) 
P = 1 pwymry 5 , | p» 2 ASU L pe 424. fy) (6.2 
} = ‘. 2 T 95 f sant > Te) 5) Aon (h, &), (0.2) 
and 

¥)2 
P r pe a te ie. xe) 

(6.3) 

- =P > Hy Ae Aes Tent = PASS (h, ky = es). 


provided, for the validity of (6.2), that the input frequency ratio p/q in (2.3) is irrational 
and, for the validity of (6.3), that the input frequency ratios p/q, q/r, and r/p in (2.4) 
are all irrational or, in short, that the inputs in question are non-periodic. The proofs of 
these formulas follow at once from the theorems cited while the justification of the 
assumptions made concerning the definition of average output power will be given in 
the last section below. We note in passing that by considering separately the leading 
terms and the summations in the left hand members of (6.1), (6.2), and (6.3) we may 
separate out the alternating and direct current output power from the average output 
Finally, it may be remarked that by applying the Sternberg 


power without difficulty. 
to the squared characteristics of more 


and Kaufman approximation process |7, 8, 9 
general continuous modulators than rectifiers, output power computations for these 
more general problems may also be carried out in terms of the same Bennett functions 
as here. 

7. Some concluding remarks. In concluding this paper on multiple harmonic 
analysis in rectifier problems it seems worthwhile to note that the simple and multiple 
Fourier series expansions of the rectifier output (2.5), (2.6), and (2.7) are also generalized 
Fourier series in the sense of Bohr [18, 19, 20] provided, in the two- and three-frequency 
problems, that the inputs are non-periodic and that v > 0. The significance of this remark 
ave for the 


and also its justification lies in the fact that under the conditions noted we | 


functions A?’ (h), AYO(h, i), and AS’) (h, ky , k.) the further formulas 
A, (h) = 2 lim z [ Y’[cos (pt + 6,);h] cos (w,t + ¢,,) db, 7.1) 
ae a . ‘ 
Aoi a 2 lim : [ } COS (pl tT O om 
T= 1 P (7.2) 


+ k cos qt + Cals h| COS (Ww t+ ¢.,,,) dat, 


i m 
A ee h. } . k.) 2 908 1’ | } |COS (pl - 6, ) (7.3) 


+ k, cos (qt + @,) + k, cos (rt + 6,); h| COS (Waamatt + Paamnt) dt, 
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where for the validity of (7.2) and (7.3) we thus assume that the frequency ratio p/q 
in (7.2) and the frequency ratios p/q, g/r, and r/p in (7.3) are all irrational and that 
v > 0. To establish these formulas we begin by noting with the aid of Hobson [15] and 
Tonelli [16] that each of the corresponding Fourier series expansions (2.5), (2.6), and 
(2.7) converges uniformly for »y > 0. To complete the proof of (7.2), for example, we 
now multiply (2.6) by cos (@amnt + @amn’) and note that the resulting series also 
converges uniformly for vy > 0 so that for each e« > 0 there exist integers M > m’ and 
N > n’ such that the remainder of the new series after the terms 

P’ AL ywlh, k) cos @aunt + day) COS Wamenrl + amin’); (7.4) 
never exceeds }¢«P" in absolute value regardless of the value of ¢. Consequently, inte- 
grating termwise up to and including the terms (7.4), multiplying through the resulting 
relation by 2/7P’, and taking limits we may thence interchange the order of limiting 
and summation operations in this finite part of the series to conclude that 


~ 7 


] F 
Azin(h, k) — 2 lim T | Y"[cos (pt + @,) 


+ k cos (qt + 0.); h] Cos (@ainrn'l + damn’) at (7.5) 
ee x 
< lim =, [ teP” dt = «, 
P’ T+ 1 J0 


all of which operations are permissible and valid provided the ratio p/q is irrational 
and vy > 0 so that (7.2) is thereby proven. Similar proofs may be given of course for 
formulas (7.1) and (7.3) also. 

Consider now the implications of the formulas (7.1), (7.2), and (7.3). First of all, 
these formulas applied to the constant terms in the various Fourier series expansions 
considered, assert that under the stated conditions time averages are equal to phase 
averages. Secondly, since squaring Y"(X; Xo) is equivalent to replacing v by 2v, it is 
now clear that the formulas (6.1), (6.2), and (6.3) previously given under the same con- 
ditions for the average output power @, were indeed correct, at least for v > 0. In closing 
we may note, finally, that in a sense we have here a simplified example of an ergodic 
situation while the input function and output series are always almost periodic functions 


in the sense of Bohr 
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ON CURVED SHOCK WAVES IN THREE-DIMENSIONAL GAS FLOWS* 


BY 
R. P. KANWAL** 


Indiana University 


Introduction. In this paper we discuss the problem of shock waves in three- 
dimensional steady rotational flows of an ideal gas with viscosity and heat conduction 
neglected and subject to no extraneous force. If a shock wave arises, e.g., as the result 
of the presence of an obstacle in the field of the flow, it will be assumed that this divides 
the flow under consideration into a region 1 and a region 2. We suppose region 1 to be 
traversed by the gas before contact with the shock surface and region 2 to be that into 
which gas enters after passing through the shock surface. It is assumed that this surface 
is given by a continuous and differentiable function of coordinates and that it has a 
continuous and differentiable unit normal which we suppose directed from region 1 to 
region 2. Furthermore, we assume the required differentiability conditions to be satisfied 
by the velocity components u; , the density p, pressure p and entropy S in each region 
1 and 2 so that it is possible to express the surface covariant derivatives of these quantities 
along either side of the shock surface. 

If the flow in front of the shock is known then the flow behind the shock is determined 
by the well-known Rankine-Hugoniot relations. The main object of this paper is to 
obtain formulas for the determination of the gradients of velocity components, pressure, 
density and entropy behind the shock surface when the flow in front is known. These 
formulas have been obtained for plane flows by Thomas [1]. In his analysis Thomas 
assumed that y, the ratio of specific heats, is constant throughout the flow. Truesdell [2] 
has extended these results to fluids obeying an arbitrary equation of state. The same can 
be done for the three-dimensional gas flows also, although we have carried out the calcu- 
lations for the case when y is constant. Truesdell has further put Thomas’s results in 
more compact form by introducing certain dimensionless variables which have been 
found useful in the following discussion also. 

The required derivatives are obtained by differentiating the Rankine-Hugoniot 
relations along the shock surface. Gauss-Weingarten formulas and various other results 
of the geometry of surfaces concerning the principal normal curvatures are found useful 
in the analysis. The complete analysis has been given only for the case of uniform flow 
upstream of the shock. As anticipated, the flow behind the shock wave is found in general 
to be rotational. The explicit determination of the vorticity components is carried out. 
This leads to the formulation of a theorem regarding the characterization of surfaces 
behind which the flow remains irrotational. It is found that the plane, the right circular 
cone, the cylinder and the developable helicoid are the only such surfaces. Furthermore, 
it is found that the component of vorticity along the normal to the shock surface vanishes 
at every point of the surface—a property shared by the contact discontinuities which 
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bo 


are surfaces across which the normal component of the velocity is continuous while the 
tangential components are discontinuous. 
Formal methods of tensor analysis have been used throughout the following discussion. 
1. Equations of motion. We consider the differential relations governing the steady 
flow of a perfect gas devoid of viscosity and heat conduction and subject to no extraneous 
force, namely [3, 4], 
Di t+ puju;,; = Y, (equations of motion), (1) 
p..u; + pu;,; = Y, (equation of continuity), (2) 
in which p, p and u; denote the pressure, density and velocity components respectively. 
We assume the motion referred to a system of rectangular coordinates x‘; then the comma 
in the above equations and in the following, represents partial differentiation. It is to 
be understood in the above and in the following discussion, unless the contrary is stated, 
that an index which occurs twice in a term is to be summed over the admissable values 
of the index. Since there is no distinction between covariant and contravariant indices 
within a rectangular system, we may write any index as a subscript or a superscript 
without modifying the value of the term in which the index occurs. 


In addition there is the equation of state, viz. [3, 4], 
p = exp (S/Je,)p’, (3) 


where S is the entropy, J the mechanical equivalent of heat and y is the ratio of two 


specific heats c, and c, assumed constant. The entropy remains constant on a streamline 


9 


but varies from streamline to streamline i.e., S;u; = 0. By differentiating (3) and using 
the Eqs. (1) and (2) we obtain 

U; , ;UiU; aan CUL,. -_ Q, (4) 
where c = yp/p, is the adiabatic speed of sound. 

2. The Rankine-Hugoniot relations. Denote by u,; , p, and p, the velocity com- 
ponents, density and pressure on the side of the wave bordering region 1 and correspond- 
ingly by wz; , p2 and p, the values of these quantities on the side of the wave bordering 
region 2. Put 

Uin - Uri:E , U2, = Ut’, 
so that u,, and u,, are the normal components of the velocity on sides 1 and 2 of the 
wave surface, respectively. Then the Rankine-Hugoniot relations for the stationary 


wave under consideration can be written as 


9 3 — se 
fu] = ets — ae 5 
—(y + I)pitin 
9; 2 
“\ PiU, — YP1) (ce 
[p] = met (6) 
’ we 


[p] 2p1(Pitin a Pt a 
2yP: + (yy — I)pitin 


Here the bracket [ | denotes the difference of the values on the two sides of the shock 


surface of the quantity enclosed, e.g., [u;| = uz; — UW; 
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3. Coordinate system on the shock surface. In the following analysis it has been 
found convenient and advantageous to have the lines of curvature as the Gaussian 
coordinate curves on the shock surface. Denote these coordinates by y’, y*. Then the 
lines of curvature are y' = constant and y* = constant and the parametric equations 
of the surface are given by 

zy = x'(y’, y’), i= - 2, ol (8) 

As is customary we shall use Latin letters for the indices referring to the space variables 
and Greek letters for the indices referring to the surface variables. Thus the Latin 
indices will assume values 1, 2 and 3, and Greek indices the values 1 and 2. 

The unit surface tangent vectors to the coordinate curves are [5a] 6%g,, and 6%q.. 
respectively where 

\1/2 
Jap = 1/(Gap)’ , 


and a,g are the components of the first fundamental form of the surface 


ax’ ax’ 
aap = a..@ 9p ? 
ay* ay 





(9) 


while 6 is the Kronecker delta. The corresponding space components of the unit tangent 
vectors are x}g,, and x3g2. respectively, where we have put 

3 

I, = =. (10) 

ay" 

Then 2; are the components of a contravariant space vector and a covariant surface 

vector. We shall call these quantities the components of projection tensor [6]. We have 

already assumed that the normal vector to the shock surface is directed from region 1 

to region 2. We further choose the orientation of the surface normal and Gaussian 

coordinates in such a way that xi , x and —' have the same orientation as the z'-, 2’- 


and 2x°-axis. 


Now consider the velocity field u; and let ¢; denote the part of the velocity field 
which is normal to the surface and let v; denote that part of u; which lies in the tangent 
plane of the surface. Evidently we have 


u;=Yy+ 6;. 
By multiplying both sides of this relation by x} we obtain [6] 


Lots = Ve; (11) 


a‘? 


‘v; = v, . We say that the tensor x} projects u,; onto the tangent 
plane to the surface, v.gaa (a not summed) being the tangential component of the 
velocity along the y* curve. 

Following Truesdell [2] we introduce the dimensionless variables: 


where we have put x 


5 —2eitin— YP) _ Peace 
= 9 ; ae ) s = as u ’ 
YP: + (y 1) pits, in (12) 
a Uin Uon 
te ™ of ’ M,, = — ’ M,, ea 
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where a is not summed, 6 is the shock strength and 7, are the components of the obliquity 
of the shock. With these abbreviations the shock relations can be written as 


— Unk; 
U;] = ae (13) 
[u;] ee 
[p] . ; 14 
p| = = Dillia ; (14) 
yi rx Pil, 
[ = -£ (15 
lp] bp; . (15) 
Multiplying both sides of (13) by £; and x; we easily deduce 
Uin Uin fan 
Un, = a (16) 
1+ 6 P2 
and 
V20 = Via (17) 


4. Differentiation of the shock relations. Differentiating the relations (5) to (7) with 
yj’ and y” we get relations along the shock surface of the form 


respect to 7 
U; fa = Uy,;%a + Ac. = AP ; (18) 
De =n1,;7%2 +B, = B, (19) 
p26 =p 7i+C, =Ct, (20) 


where for simplicity we have omitted the subscript 2 on the quantities appearing in 
the left members of these equations. The explicit formulas for the quantities A;, , B, 
and C, are readily obtained by differentiating the right members of the equations (5) 
to (7), respectively, with the help of the Gauss-Weingarten formulas of differential 
geometry (see Sec. 6). We see that these quantities are expressible in terms of the known 
parameters of the flow in the region 1 and their partial derivatives along the shock 
surface, the projection tensor x; , the surface normal é and the two principal curvatures. 
Following Thomas [1], any quantity which can be so expressed will be said, in the follow- 
ing, to have an allowable determination. 

5. Derivation of the partial derivatives behind the shock surface. Consider the 
Eqs. (1), (2), (4), (18), (19), and (20). From this set of six relations we can determine 
the partial derivatives of the velocity components u, , 
immediately behind the shock surface. For this purpose we first eliminate the derivatives 
to obtain the following equations, which we represent as two sets of equations for 


Pi 
convenience of reference, namely 
pst, = C2, (21) 
pu; + pu. = 0, (22) 
and 
ue. = Al, (23) 
pu; ;2qu; = —Bt, (24) 
pu; UU; — YP, = O. (25) 


the density p and the pressure p‘ 
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Now from (23) to (25) we can obtain u;,; and then from (21) and (22) p,; can be found; 
after this the Eqs. (1) yield the values of the derivatives p,; . To effect the above deter- 
mination we define a matrix || C;; || such that 
Cia=Na; Ci = Ui/U,, (26) 
where \,, denote the space components of the unit tangent vector to the curve y* and 
are given, as remarked above, in terms of the projection tensor by 
Nse = Aa = BaGae no summation on a. 

The determinant of the matrix || C;; || is easily seen to be unity. Hence we can define the 
quantities D;, by 

Dili; = 6;; ; D, Ci = 63; . (27) 
The second relation (27) follows from the first relation and conversely. Since we shall 
have to use the explicit form of the matrix || D;; ||, it can be easily seen to be given by 


| (Astls — Ages) (Altus — Afus) , 























Zz u u = 
D:, Da Da || , , | 
| 3 = 1 = 
Djs Dos Des | af (Asu, Ass) (AjyUs Aju) t, | (28) 
i ! Un Un 
13 Dos Dsl} |i, —— 
: Dis Ds ! (Atte — AZu,) (AZu, — Alea) P | 
| Un U, oe 
Now define the quantities B,;,; by 
B;; _ Uy mC eC mj ° (29) 
Then from (23), (24), (25) and (29) we have 
Bt gas 
Bas — LaABIJaaGss ’ By; _ a ’ 
_ A tatiJaa _ Uk .k 
Bsa — Un ] By, — M , 


where a, 6 are not summed. 
Furthermore, by multiplying both sides of (29) by D;,D;, and using the relation 


(27) we readily obtain 
Ut ,m _ B,D Dim ° (31) 
Now the first eight elements of the matrix || B;; || have an allowable determination 


as can be readily seen. However, B,, is given in terms of the u,, ; it therefore remains to 
determine the quantity B,, . But this can easily be done by contracting the indices 


and m in (31), viz., 
Use = B;;,DuDir ’ 
which, by the help of the matrix (28), gives the value 


Ux. = B,(1 + x1) + B,(1 + X2) + 2B 2)XiX2 = 2B as)X1 uF 2B .23)X2 + Bg, ’ 
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where x’s are defined by (12) and for brevity we have introduced 


Buin = UB, + B,.. 


Thus we arrive at the value of B,; 
B B,, + x1) + B,.(1 + x>) + 2Buwxixe — 2Buax: — 2Bi23)x2 (29) 
~~ = — => menaecimcacanien 4 

M,; -1 
which furnishes an allowable determination of B,, . Hence the relations (31), in which 
B,; is given by (32), provide an allowable determination of the gradients of velocity 


components behind the shock wave. 
To effect the determination of the gradient of density behind the shock surface we 
observe that the set of equations (22) and (23) can be written in the form 
pCi, = d;, 
where 
d, = giuiC% ; dy = Jo2C3 
d; = —pu,.,./U, = —pB,,M;/u, . 


n 


Hence we have 


WwW 
wn 


p= & By. (33) 


as the allowable determination of the quantities p,; . Finally from (1) the derivatives 
p.; behind the shock surface can be determined by the equations 


Dp. = —puu;,; = —pu,B;,D,, (34 

6. Calculation of the invariants d; and B;;. It is possible to give an explicit formu- 
lation of the invariants d; and B;; under the assumption that the flow in front of the 
shock is uniform, i.e., the velocity, pressure, and density are constant in region 1. To 
effect this formulation we need some results of the differential geometry which are briefly 
stated below [5a]. 

We have already noted the expressions for the first fundamental form of the surface, 
the projection tensor and the space components of the unit tangent vectors to the 


coordinate curves. They are 
Ox 
Gap = LqXp ; rt, = Aa 
OY 


= Jaala ; a not summed, 


respectively. In addition we need the expression for the normal curvatures in the direc- 

tions of the coordinate curves (which are the lines of curvatures and hence these curva- 
tures are the principal normal curvatures). They are 

b= Cafe, 3 he. @ BJ acs, (35) 

where b,, are the components of the second fundamental form of the surface, i.e.; 

6 4s 97 

bas = 3€ CijLa pl Ls ° (36) 

. , . . . . 

In the above relation e”” are the components of the surface permutation tensor with 

the properties 
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Mea Pa 1/o; & = -1/@", 


where a is the det | das ||, while e;;, are the components of the permutation tensor 
of the space with the properties 


= (0 if any two of the indices 7, j and k are equal, 





(i C554 
(11) €123 = C231 = Cri2 = +1, 
(111) €i30 = C321 = Cris = —1, 


‘ .s denotes the surface covariant differentiation of x} . As we have taken the lines 


and 2° 
of curvature as the coordinate curves on the shock surface we have 


Qe = 6, = @. (37) 
| Moreover if a*” denote the components of the tensor conjugate to a,, , namely, 
aga" = 83, 
we obtain 
a'' = a,,/a, a* =a’ =0, a®™® = a,,/a. (38) 


In this notation the space components of the unit normal vector to the surface are 


given by the relation 


8 k 
£° = 6° CipeLatg . (39) 
The well-known Weingarten formulas give 
7 78 
Ba = —O bed ; (40) 


where ; denotes, as before, the surface covariant differentiation. Keeping these results 
and the assumption of uniform flow in front of the shock in mind, we get 


: 
Uinsa = Unikie = —U,,07"b,y2, = 08""b,, , (41) 


where we have used the relations (11) and (17). Expanding this result and making use 


of the expressions in (35) and (38) we obtain 


Dicia —v,k, ; Uinsa = — Voke. (42) 
Similarly, the relation (40) yields 
é. = —&k,2; ; Ein = —kyx,. (43) 


Now for the uniform flow 


Hence we can readily find [7]. 


lin = (rate _ aes) 
4iia —(¥vy +- 1) pitin “5a (44) 


_ — {Bat ‘ Ip (—kaXatin + go.k.)}, 


Pillin PiUin 
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in which there is no summation on a. 


From (14), (30), and (44) we obtain 
B,, = 40% 


?| @2; . 

palma Sie = lp] ky-t,°—— ™ bu,k, ’ 
3) Pitin Qa, 

Bi» = By; = 0, 


(45) 
Baz = Oke « 
9) 


Also from (6) and (42) we get 





Ap UinWaka , 
B, = — ee no summation on a, 
~ + I 
yielding 
Big $9, ,v,k, 
B, = — oo - (46) 
pu, 77 % 
and 
. B 22 Agovok ~ 
B., = — = (47) 
pu, y+ 1 
where the relation (16) viz., p,t,, = pote, has been used. 
Furthermore, from (44) we obtain 
u ] 6 Uin 
Aiea = =A Bo + ——- | —2, £.° + vk. }?, 
Un Pitiiy 1+ 6 Un 
i 
in which there is no summation on a. From this we get 
A sjsU411 4 0 
Bs, = ag Vig ky a a <)? (48 
Ur yt! 1+ 6 
and 
sath ee 4 5° 
B., = ee = Jooks [artis ak-eereed F (49) 
un, _ + ] ] —- @ 
Lastly from the formula (32) we get the value of B,; as 
/ / \ 5 
j 8 hh ae 
Bz, = 4 6u,(1 + x3) — {J — - }y Ve. 
: ( +1 1+ 6/ ) © E 
(90 
l J ¢ (1 < 2 | ‘e) I 5° Nh 
+ 4 6u,(1 + x2) — + Xs} fhe - 
yt+1°14+e/J 
The matrix || B, is thereby completely determined. 
Correspondingly the d; are given by the expressions 
d Ci9n = 2g vik 


(y + l)u 


¥~ ip & + Dal 
292k 
2922 = Gt+1 Us, iy — Ip — vy + Ips 


§» 


pB;,M? 
d, = ——>—". 


(51) 
Un 
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From the relations (45) to (50) and the relations (51) the invariant character of 


B,, and d; can be easily seen. 
7. Variation of velocity, density and pressure along the shock surface. The rate of 
change with respect to the are length along the lines of curvature of the velocity com- 


ponents u,; is given by the following formula 
Us: LaJaa = BimDiiDniLaGaa 5 a not summed. 
When expanded this becomes 
“u; A; = ByD,; , ua; As = B,D; « (52) 


Similarly the rate of change of density and pressure is given by 


and 
D.irte = —pu,B;,D, ir. , 
which, when expanded, becomes 
pA; = —pu,B,; , D.iA2 = —pu,B,; . (53) 
From (52) we readily find 
gg. = uu; Ai = By , (54) 
99.2 = uu; A> = Bs , 


where g = u,u,; . Thus when the flow in the region 1 is uniform the expressions for these 
quantities are known from the results of Sec. 6. These results moreover give the physical 
meaning of the invariants B,;; and d; . For example, B;, and B,, are the rates of change 
of half the square of the velocity magnitude along the lines of curvature as shown in (54). 


Furthermore from (53) we get 


p.; = p'u2B;sBisD,.Di: 5 
Df p 53D) F k (55) 
= pun {Bis i B:, i (—Bisxi - BosX2 + B, )*}. 
Thus p,, 0 only if B B,, = Bs, = 0. But from (46), (47), and (12) we get 
47,u,,k 4M wks 56 
a. a 47M Uink , B. = —— “2 (56) 


7. * ¥+1 

Thus we find that Truesdell’s [2] theorem in plane flows is carried over to spatial gas 

flows, viz., when a uniform flow of a fluid crosses a shock, the pressure gradient cannot 

vanish at any point on the rear side of the shock if the shock is curved and oblique. 
8. Determination of the derivatives of entropy behind the shock. The value of 

the entropy S as given by the relation (3) is 


S = Je, log (2). 


p 
Thus 
[8] = Je, log {+ Sip + Sous) ait 
pill + 4) 
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where we have used the shock relations (6) and (7). 
In the case of uniform flow in front of the shock we derive, by differentiating both 
sides of (57), the result 


\ 


( 5) 5.u7. \\ 
S We (ve. log ( Re me *) ( E. ’ (58) 
fit a9 Ji 


where £, can readily be evaluated with the help of the relations (15), (20) and (42). 
We also know that 


which states that entropy is constant along the streamlines. This result and the relations 
(58) can be put in a compact form with the help of the matrix (26), as 


a f, (59) 
where 
fii =E£,, f.=E£., f, = 0. (60) 
Inverting the relation (59) with the help of the relation (27) we obtain 
S,=f.De. (61) 


9. Curvature of the streamlines behind the shock. The curvature Kk of a streamline 
in the flow under consideration is given by [8], 


2 


~ I: ig 
Ku; = ~(2s ~= teed), 
\pg qd / 
where A; and yu; are the components of the unit tangent vector and principal normal 


vector to the streamlines. Hence 


q\ p q 


When we make use of the results (32) and (55) we get the value of K’ as 


9 { 9 
~_ a ae . > B:. ) 
K* = -34Biz + Bis + (—Biax: — Bosx2 + B33)? — ———*—. 
7 | , l+xi +x 
In this work we have confined ourselves to the evaluation of the first derivatives of 
the parameters of the flow behind the shock, while to determine the torsion of the stream- 
lines we require as well the second derivatives [8] which have not yet been calculated. 
10. The expression for the vorticity vector behind the shock. The components of 
vorticity behind the wave are given by 
i zh ijk / i 1 5° ¢ 
w= 6"; = é' BB; DyD., (—kyygisA2 + kev2go2\3) -—— (62) 


1+ 6 
Now for w' to be zero it is necessary as well as sufficient that either 6°/(1 + 4) or 
(—kyigirrz + kvego2d2) 


vanish. In the latter case coefficients of both \j and A} must vanish because \j and i} 
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1959] 
are perpendicular to each other. Thus in order that the flow be irrotational behind the 
shock surface at least one of the following conditions must be satisfied 

Q 6 0, no shock wave; 

(ii k, and k, are zero, i.e., the shock surface is a plane; 


v, and v, are zero, i.e., the obliquity is zero which again leads to a plane; 


As (iv) and (v) lead to the same result, we will discuss, for the sake of definiteness, the 
case (Vv), viz., k, = 0 = v, . k, = 0 implies that the shock surface is a developable surface 
and as such its generators and their orthogonal trajectories form its two congruences 
f lines of curvature [5a]; while v, = 0 implies that the orthogonal trajectories to its 
‘nerators are plane curves lying in the planes which are normal to the direction of the 


given uniform flow in front of the shock wave. 

Now with the exception of cylinders and cones every developable surface is the 
tangent surface of some curve [5b] and the orthogonal trajectories of the tangent surface 
of a curve are the involute of the curve [5c]. But the necessary and sufficient condition 
that the involutes of a twisted curve be plane curves is that the curve be a cylindrical 
helix [5d]. Moreover, the planes of the involutes of a cylindrical helix are normal to the 
generators of the cylinder on which the helix lies [5e]. Therefore, if we take the helix 
to lie on the cylinder whose generators are in the direction of the flow then the tangent 
surface of this helix, developable helicoid, satisfies the condition (v). As far as a cone 
and a cylinder are concerned it can be easily seen that the right circular cone with its 
axis parallel to the direction of the uniform flow and any cylinder with generators parallel 
to the direction of the uniform flow satisfy the conditions (v). We thus get the theorem: 
the only shock surfaces behind which the flow remains irrotational are, a plane, a right 
circular cone, a cylinder and a developable helicoid. Among these shocks the conical shocks 
have been extensively studied [9]. In the case of cylindrical shock surface whose generators 
are parallel to the direction of the uniform flow it is obvious that there is no discontinuity 
in the flow because flow is merely tangential to the surface. 

Another result follows if we multiply both sides of the Eq. (62) by £; getting wf; = 0. 
Hence, when a uniform flow breaks across a shock, the component of vorticity along 
the normal to the shock wave vanishes at every point of the surface. It is interesting 
to note that this property is also shared by contact discontinuities, i.e., vortex sheets, 
which are singular surfaces across which the tangential component of the velocity 
changes abruptly while the normal component of velocity is continuous and the expres- 


J 


sion for vorticity components on the sheets is given by [10] 
1 ijk 
w =e''E ful, 
u,| being the jump in velocity across the vortex-sheet. The component of vorticity 
normal to the three-dimensional pseudo-stationary and unsteady shock fronts has also 


been found to vanish [11, 12]. 
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Note added in proof: The author has recently extended these results to unsteady shocks 
in fluids obeying an arbitrary equation of state [12]. 
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SOME INTEGRATED PROPERTIES OF SOLUTIONS OF THE WAVE 
EQUATION WITH NON-PLANAR BOUNDARIES* 


BY 
LU TING 


Department of Aeronautical Engineering and Applied Mechanics, Polytechnic Institute of Brooklyn 


Abstract. Integrated properties of solutions of the wave equation with non-planar 
boundaries are found and applied to three dimensional supersonic flow problems and 
two dimensional diffraction problems. 

For the problem of supersonic flow outside a cylindrical surface with generators 
parallel to the flow direction, a theorem is proved concerning the integrated properties 
of the linearized pressure distribution and the prescribed normal velocity on the surface. 
The theorem is a generalization of the integral relationships obtained previously and is 
useful in the evaluation of total lift and drag of wing-body combinations when the linear 
dimensions of the cross section of the body are not small as compared to the chord length. 

For the diffraction of a pulse or a weak shock over a rectangular notch, a pressure 
integral theorem is obtained. Its usefulness is demonstrated in reducing the labor of 
obtaining the mean pressure distribution along any depth inside the notch at different 
instants for various width-height ratios of the notch. 

Introduction. Theorems concerning certain integrated properties of the linearized 
pressure field due to planar source distributions in a supersonic stream have been pre- 
sented by Lagerstrom and Van Dyke [1] and by Bleviss [2]. Extension of the theorems 
and their applications to a class of three dimensional problems involving biplanes or 
cruciform wing arrangements were presented by Ferri [3] and by Ferri and Clarke [4]. 

Ferri, Clarke and Ting [5] obtained a theorem regarding the pressure integral along 
the line of intersection of a forward Mach plane with a specified non-planar surface 
which represents a prismatic body of rectangular cross section mounted on a planar 
wing with supersonic edges. The pressure integral is related to the integral of the source 
distribution which in turn is related to the integral of the prescribed normal velocity 
on the top (or bottom) surface of the body and that on the wing surface. A similar rela- 
tionship can be established if the normal velocity is prescribed on the side walls of the 
body. The wing surface together with the surface of the body above (or below) the 
wing represents a cylindrical surface in the form of a “single’”’ step with generator parallel 
to the direction of the undisturbed supersonic flow. 

By a further extension of the theorem in [5], it can be shown that a similar integral 
relationship is valid if the cylindrical surface is of the form of a “stairway” with a finite 
number of steps. As the number of steps becomes infinite and the size of each step 
infinitesimal, then a sequence of “stairways” is obtained whose limit can approximate 
the shape of any given cylindrical surface. Hence, the generalized integral relationship 
to be presented in this paper is conjectured. However, the argument which leads to 
the conjecture cannot be accepted as a proof of the validity of the relationship unless 
it can be shown that the limit of the sequence of the corresponding disturbance pressures 
or potentials exists and equals the corresponding value for the given cylindrical surface. 


*Received November 12, 1957. 
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Such a convergence proof is by no means simple. An attempt to verify the relationship 
directly by the same procedure used in the preceding investigations has not been suc- 
cessful owing to the difficulty of finding the proper source distribution which corresponds 
to the prescribed boundary condition on the cylindrical surface. The fact that the 
velocity potential is a solution of the wave equation is expressed implicitly in the preceding 
investigations through the relationship between the integral of the pressure distribution 
and that of the source distribution. In the present paper the integral relationship between 
the pressure distribution and the prescribed normal velocity is verified by observing the 
fact that the velocity potential obeys the wave equation. Thus, the necessity of finding 
the proper source distribution is avoided. 

Generalization of integral relationship. As shown in Il'ig. 1, a cylindrical surface, 
y = F(z), is placed in a supersonic stream with its generator parallel to the z-axis, 


Cylindrical surface, y= F(z), abcdnma 


Mach plane, x+By= M8 ehlie 

Pione x0, begfb 

Plane z= K,(8) ckgc 

Plane z=Ke(f) bjfb 

Domain of surface integral, If, shaded area 
Path of line integral, r, ijkl 

















Ml stm, 


Fic. 1. The generalized integral relationship 


which is the direction of the undisturbed supersonic flow with velocity U and Mach 
number M. q,[z, y F(z), 2] represents the prescribed small normal velocity on the 
cylindrical surface with 

g,/U «1 (1) 


and 


q,(z, F(z),z] = 0 for zx <0. (2) 
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It will be assumed that q,[x, F(z), z] is piece-wise continuous and that for any given 
Mach plane, x + By = M8, at a finite distance 8 from the origin, there exist two finite 
values K,(8) and K,(8) with K, > K, , such that on the Mach plane, the disturbance 
potential g, due to q, , is confined inside the region K, > z > Kz, Le., 


g(x < Mp — By, y,2z) = 0 


for 
o >z>K,(6) and K,.(g) >z> —o. (3) 


The disturbance pressure, p(x, y, z), is related to the prescribed normal velocity on 
the cylindrical surface by the following integral relationship. 

Theorem I. The integral of the y-component of the disturbance pressure force acting 
on the portion of the cylindrical surface, y = F(z), which lies ahead of any given Mach 
plane parallel to the z-axis, x + By = MB is equal to pU/B times the integral of the 


prescribed normal velocity on the cylindrical surface, over the same domain of inte- 


gration, 1.€. 


a ure 
I p(x, y,zn-jdS - ‘3 I q(x, y, 2) dS, (4) 
. , . 


I r 


. 


where n is the unit vector normal to the cylindrical surface, j is the unit vector parallel 
to y-axis, and T’ represents the domain of integration on the cylindrical surface for 


0<2< MB — By. 
Proof. The disturbance potential g(x, y, z) obeys the wave equation 


x a “2 a 
oo a oe ery (5) 


dx aye dz 
If a vector Q(z, y, z) is defined as 
Q = (Be, + ¢,)(Bi — j) — ok (6) 
then Eq. (5) yields 
div Q = 0. 
If the volume which is confined by the cylindrical surface y = F(z), the Mach plane 


x + By = Mg, the planes z = K,(8), z = K.(8) and x = 0, is designated by V and 


its surface by &, then Gauss’ theorem [6] states 


—[[a-qus = [ff aiv qav =o, (7) 
; ! 


Ss 


where n denotes the unit vector normal to the surface S and points inward. Equation (7) 
is valid if Q is continuous inside JV, i.e., if ¢ and its first derivatives are continuous. 
This is certainly true when q,[z, F(z), 2] is continuous. It will be shown later that Eq. 
(7) is valid even if q, is piecewise continuous. 

Since ¢ and its first derivatives vanish on planes x = 0,2 = A, and z = K,, Q van- 
ishes thereon. Since the unit normal vector n; to the Mach plane « + By = Mg is 
n, = (i+ Bj)/M, it is clear that on the Mach plane —Q-n = Q-n, = 0. Consequently 
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on the surface S, n-Q vanishes except the portion on the cylindrical surface, and Eq. 
(7) becomes 


The domain of integration on the cylindrical surface has been extended beyond the 


region of K, < z < K, because ¢ vanishes for z > K, orz > K,. 
For the cylindrical surface y = F(z), n-i 0, and Eq. (7) becomes 
|} (Be. + ¢,)in-j) + en-k] dS = 0. 9) 
With disturbance pressure p = —ply, and the disturbance velocity, 


q=ei+ed tock, 


Kq. 9) leads to 


B Se 
I 

iy [fonds = ffnas. 
I 3 


The remainder of the proof of the theorem is to show that Eq. (7) is valid even when 
qn is only piecewise continuous. 

If ¢,{z, F(z), z] has a jump discontinuity across a curve on the cylindrical surface I, 
and the first derivatives of the disturbance potential are discontinuous across a surface 
S. inside the volume V, the discontinuity surface S, will be a characteristic surface—a 
Mach cone or an envelope of Mach cones. 

A unit vector, n, , normal to the discontinuity surface, S, , is also normal to a Mach 


cone and can be written in general as 
n, = (i — B cos \j — Bsin \k)/M (11) 
where \ is a parameter. Due to the discontinuity surface, a term 


[/ 6(Q-n,.) dS 


e 


should be added to the right side of Eq. (7). Here 6 express the jump of the quantity 
inside the parenthesis across the surface S, 
From Eqs. (6) and (11), it follows: 


Q-n. = [(Be, + ¢,)(1 + cos A) + ¢, sin AJB/M 


= grad ¢-t, 
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where 
tr. = [((Bi + j)\(1 + cosa) + sin AKIB/M 
and ¢, is the derivative of ¢ in the direction of <, . 
With «.-n, = 0, or +, tangential to the surface S, , g, represents a derivative of ¢ 


along the surface S, . Since ¢ is continuous across the surface S, , ¢, is also continuous, 
i.e., 5(¢,) = 0. Hence, the additional term to Eq. (7) vanishes, 


[| 6(Q-n,) dS = I |e. | 6(y,) dS = 0 


and Eas. (7) to (10) are valid when gq, is piecewise continuous.* Thus concludes the proof 
] \ \ f| 


of the theorem. 
If T, and [, are the areas on the cylindrical surface y = F(z) ahead of the Mach 
- B,M and x + By = £6.M, respectively, the theorem or Eq. (4) yields 


[| pn-jdS = v I q, aS. (12) 
"4-Ts r.-T. 


r, 


planes x + By 


When tle two Mach planes are very close to each other, i.e., 1/8, — MB, = dx, the 
following result is obtained. 

Corollary. The integral of the disturbance pressure along the curve of intersection, 
T, of the cylindrical surface y = F(z) with the Mach plane z + By = MG is related to 
the integral of the prescribed normal velocity q,[z, F(z), z] as follows 


[ psin on-jdT = - [ q, sin o dT, (13) 
J1 B Jr 


T 


where o represents the angle between the z-axis and the tangent to the curve T. 


Equation (13) is obtained readily from Eq. (12) when dS is replaced by dx dT sin o 


and the area T, — I, becomes a strip along the curve T with thickness equal to dz sin o. 
On the basis of linearized theory, the prescribed normal velocity g, on the cylindrical 
surface y = F(z) can be created by deforming the cylindrical surface slightly into a 
surface y = F(z) eE (x, y, 2), if 
qn(z, F(z), 2] = eUE,|x, F(z), z)(1 + F2)'” (14) 
where e < l. 


With this equation the integrals of disturbance pressure on an almost cylindrical 
surface in a supersonic stream are obtained from Eqs. (4) and (13). 

Application to wing-body interference. In studying wing-body interference, the 
body is usually represented by an almost cylindrical surface with generator parallel 
to the direction of the supersonic stream (z-axis). The wing, which is planar with super- 
sonic edges, coincides with the z-z plane while the z-y plane is the plane of symmetry. 
The surface of the body above (or below) the wing together with the upper (or lower) 
surface of the wing form the cylindrical surface y = F(z) (Fig. 2). 


*This is a sufficient condition but not a necessary condition for the validity of the theorem. 
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SHADED AREA: DOMAIN OF INFLUENCE OF BODY ON WING 
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Fic. 2. Wing-body interference 


If 7, and 7’ y represent the line of intersection of the Mach plane, zr + By 8M, 
with the body and the wing respectively, the corollary, Eq. (13) becomes: 


j ' — i ty = : 
| psin on-j dT + | p dz 7" | g,sino dT + [ [ d dz, 15) 
J TE JT > JT; JT 

where #(z, z) is the inclination of the wing surface with respect to X-AXIS. 


This relationship can serve as a check to the numerical results of the pressure distri- 
bution obtained by analytic methods [7, 8, 9, 10}. 

If the wing has a straight unswept trailing edge, the total lift acting on the region I 
of the upper (or lower) surface of the wing-body combination ahead of the Mach plane 
passing through the trailing edge is expressed by the term on the left side of Eq. (4). 
This lifting force which consists of the total lift acting on the wing and the lift acting 
on the part of the body is equal to the integral of prescribed normal velocity represented 
by the right side of the same equation. Therefore, if the pressure distribution on the 
body is obtained by analytical methods [7, 8, 9, 10], the total lift and the center of 
pressure can be obtained by virtue of the theorem and the corollary, i.e., Eq. (15). The 
labor of calculating the pressure distribution on the wing is saved. 

In the special case where the inclination of the wing surface inside the region of 
influence of the body is constant in the spanwise direction, the total drag can be evaluated 
with the help of the theorems without calculating the pressure distribution on the wing 
due to the interference of the body [5, 11 

For wings with subsonic edges at an angle of attack, the theorem, in general, will 
not be he Ipful in the evaluation of lift or drag. However, for the special case of a rectan- 
gular wing, the theorem can be applied to obtain expressions for the lift, the drag, and 
the center of pressure of the wing without calculating the detailed pressure distribution 
[11]. 

Integral relationship in diffraction problems. Since the linearized two-dimensional 


unsteady flow obeys the wave equation, 
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7” Vo 1 0 : 
<¥ 4+", =0, (16) 
Ox oy C~ at 


it is expected that an integral relationship similar to that for the linearized supersonic 
three-dimensional flow can be found. In order to demonstrate its usefulness, an integral 
relationship is established for a special diffraction problem—the diffraction of a weak 


shock or pulse over a rectangular notch. 
As shown in Fig. 3, a plane weak shock (2 = Ct), which travels parallel to the ground 
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Fic. 3. Diffraction of a weak shock over a rectangular notch. 


(y 0) is passing over a two-dimensional rectangular notch of width w and depth A. 
The pressure, density and speed of sound of the air at rest in front of the shock are P, p and 
C respectively. Deviations from the states at rest are defined as disturbance quantities. 
For example, if the notch is absent, the disturbance pressure and potential are 
ePH(Ct — x) and «C(1/y). (2 — COH(Ct — x) respectively, where « is the strength of 
the shock and H(Ct — x) is the unit function, i.e., H = 1 for Ct > xz and H = 0 for 
Ct < x. The disturbance due to the notch will be confined inside the forward cone 
Ct (2° + y°)'*. Across the cone, the disturbance potential and its first derivatives 


are continuous. 
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Experimental investigations of this problem in shock tubes were reported by Smith 
[13] and Coulter [14]. Theoretical solution to this problem can be obtained by the method 
outlined in [9] and [12]. Nevertheless, an integral relationship will be established and 
save much labor in getting numerical results. 

Theorem II. Twice the value of the line integral of disturbance pressure across the 
notch at a depth d(0Q > y = —d > —A) at an instant ¢ equals the difference of the 
integral of disturbance pressure on the ground from the corresponding integral if the 
notch is absent at the “retarded” instant t = tj — d/C plus the same at t = t, — 
(2h — d)/C, i.e., 


nw Ct 0 ct 
2 | p(x, —d, to) dx = |e [ H(Ct — x) dx — | pdx — / p az |. 
d Ct w t=ta-—d/C (17) 


70 Ct 


aCt +0 aCt 
+ | | H(Ct — x) dx — | pdx — | p az |, ; 
J-Ct J-Ct Jw tet o—(2h-d)/C 
Proof.* With characteristic coordinates a = Ct — yand B = Ct + y, the wave equation 
becomes 


4Gasp = Gre: (18) 


The a- and §-axes together with the z-axis form an orthogonal system. The plane 


8 = Ct, — d = const. 8 passes through the depth y = —d at the instant ¢ to . On 
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Fic. 4. Integral relationship for diffraction problems. 


this plane, the cone is represented by a parabola x* = a@, and the shock front is repre- 
sented by a straight line 2x + a = £8, which is tangential to the parabola at point x = 8, 
(Fig. 4). 

If S designates the area confined by the parabola, the ground a = £, , the side walls 


*For the convenience of those who are interested in diffraction problems, the proof is given here 


without relying on the previous theorem. 
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of the notch z = 0 and xz = w, and the line a = 8) + 2d at the depth y = —d, then 
the divergence theorem in the plane 8 = £, states (p. 88, Ref. [6]): 


I 


[ (4e5 dx + e. da) 


| [4gg cos (A, 2) — ¢, cos (A, a)] dr 


= | (4¢as — Per) dS -_ 0, 
8 
where A designates the contour of the area S. If A, , denotes the part of the contour on the 


ground, A, on the side walls, A; along the depth y = —d, and A, and A; denote the 
portions of the parabola above and below the ground level respectively (Fig. 5) then 


A = A, + A, + As + Agt+ As. 





p=0 p z=«P P=0 
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fh = Shock front 





Fig. 5. B-planes, Ct + y = Cty — d. 


On the side walls, ¢. = 0 and cos (A, x) = O, therefore, the line integral of Eq. (19) 
along A, vanishes. Across the parabola (Mach cone), ¢ and its first derivatives are 
continuous, therefore, along A; , gs , ¢, and the corresponding line integral vanish.* 
Along A, , ¢. = €C/y and gz; = —eC/(2y) and the corresponding line integral becomes: 


. 


| [4g3 dx + vy, dx] = 4eC(Ct, — d)/y. 
JA 


*When Ct > d + w, As disappears from A, however, it does not affect the result, because the line 
integral along A; vanishes anyway. 
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Along A; , ¢, = 0 and gs = (vy, + ¢./C)/2 = —p/(2pC), and Eq. (19) becomes 
oy nCto-h 
| p(x, 0, tp — d/C) dx + | p(x, 0, te — d/C) dx 
e (Cte-d J u 
— | ePH(Cl, — d) dx (20) 


= — | D(x, —d, to) dx + pC | g(x, —d, ty) dx, 
the terms on the right side corresponding to the line integral along A; . lor t) < (d~ + 
w’)'’*/C, As has been extended to the full width of the notch, because the integral in 
the additional segment vanishes. 
By the method of images, the bottom of the notch plane y = —A, can be removed 
and considered as a plane of symmetry, 1.e., 


g(a, y, 2) = glx, —(2h + y), 2]. (21) 
Instead of the plane 6 = Ct, — d, the other characteristic plane a Clo + dis 
considered and the result corresponding to Eq. (20) is obtained, 
> nCt »Ct 
| p(x, —2h, t) dx + | wr, —2h, i) — | ePH(Ct — x) dx 
« Ct uw « Ct b=le rh i 4 (99 


.u aw 


= = | p(x, —d, to) dx — pC | y,(x, —d, ty) dx. 


/0 


With the condition of symmetry Eq. (21), the sum of Eq. (20) and Eq. (22) is equiva- 
lent to Eq. (17). Hence, the theorem is proved. 
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The mean value of disturbance pressure across the width of the notch at depth, 
/ —d, is defined as 


p*(y = —d,) = : | p(x, —d, t) dx 


JO 


and Theorem II, Eq. (17), becomes 


nCt »0 Ct 
Pry —d,t = r. | ep | H(Ct — x) dx — | Pp dx — / Pp ir |. 
aW J-Ct J-Ct w t=te-d/C 
Ct 0 ct (23) 
+ eP | H(Ct — z) dz — / pdx — / Pp ar |, : 
2w J-ct -Ct ” t=te—(2h—-d)/C 


With this relationship, the mean value of disturbance pressure along any depth 
inside the notch can be obtained from the pressure integral on the ground. The latter is 
obtained by the method presented in [9] and [12]. Detailed steps are given in [11]. The 
variations of mean pressure across the width of the notch at y/h = 0, —}, —}, —? and 

-1 with respect to time are shown in I’igs. 6 and 7 for notches with width-height ratios 
equal to 1 and 3 respectively. 
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Concluding remarks. lor the problem of supersonic flow outside a cylindrical surface 

with generators parallel to the flow direction, a theorem is proved concerning the inte- 
grated properties of the linearized pressure distribution and the prescribed normal 
elocity on the cylindrical surface. This theorem can be extended readily to the case 
where there is more than one cylindrical surface present in the flow field and will be 
useful in the evaluation of total lift and drag of wing-body combinations and additional 
cylindrical bodies representing engines and (or) pylons. 

For the diffraction of a pulse or a weak shock over a rectangular notch, a pressure 
integral theorem is obtained. By the same procedure, the theorem can be generalized 
for a notch of any shape, although the analytical solution of the diffraction problem is 
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not available. By virtue of the theorem, the time history of the mean pressure along 
any depth inside the notch can be obtained from the measured pressure distribution 


on the boundaries. 
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VIBRATIONS OF TWISTED BEAMS II* 


BY 
WILLIAM BOYCE AND GEORGE HANDELMAN 


Department of Mathematics, Rensselaer Polytechnic Institute 


1. Introduction. The basic equations governing the transverse vibrations of a 
straight, twisted beam rotating about an axis, passing through one end and perpendicular 
to the undeflected central axis of the bar, have been discussed in an earlier paper [1].** 
It is the purpose of this note to explore some of the implications of these equations and 
to answer some unresolved questions raised in the previous work. We shall not be con- 
cerned with the problems of numerical calculations, for these have been studied exten- 
sively by others (see, for example, [2] and [3]). The results presented here are of a some- 
what more qualitative nature. 

The twisted beam is described in terms of a straight center line which is the locus of 
the centroids of the cross-sectional planes taken normal to the line. A cross-section is 
specified by means of the are length s measured along the center line from a fixed origin 
0 on the axis of rotation. Two triads of orthogonal unit vectors, as shown in the figure, 











are used in the analysis. The first is a moving triad i, j, k, in which i is directed along 
the undeflected center line at a generic point Q, positive in the direction of increasing s; 
whereas, the vectors j and k have the directions of the principal axes of inertia of the 
cross section at Q. In general, the triad rotates about the i axis as the center line is 
traversed with j(s) at Q making an angle ¢(s) with j(0). The second triad of unit vectors 
is a rotating, untwisted frame consisting of I, coincident with i; J, lying along the axis of 
rotation; and K, perpendicular to I and J forming a right-handed system. The vector 
j(0) forms the angle @, with J. 

If the displacement u(s, ¢) is assumed to have harmonic time dependence, u(s, t) = 

v(s) exp (7At), and w(s) = i X v, the equation of motion [see [1], Eq. (4.8)] is given by 
da’ f d-w\ .af dw | 2 ae 
5{\B: ? r <L(s) ——? — m(s)A(s)Q°@-w — m(s)A(s)\X w = 0. (1.1) 

ds | ds } 

*Received November 29, 1957. This work has been partially supported by the Air Force Office of 
Scientific Research of the Air Research and Development Command under Contract No. AF 18(600)- 
1586. 

**Numbers in square brackets refer to the list of References given at the end of the paper. 
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The mass per unit volume and cross-sectional area are denoted by m(s) and A(s) respec- 
tively, while Q is the constant angular speed of rotation. When referred to the i, j, k 
triad, the dyadic B™ is given by 


B” = El,jj + El.kk, 


where /(s) denotes Young’s modulus, and J,(s) and J.(s) are the centroidal moments 
of inertia about the j and k axes respectively. The dyadic @(s) is defined in terms of 


A(s) d(s) + 6 
Cis cos” 6jj — sin 6 cos 6(jk + kj) + sin’ 6kk; 
and 


L(s) = | m(é) A(E)E dé, 


where / is the total length of the bar. 
The beam is assumed to be elastically supported at the axis of rotation. Consequently, 


w = w’ — ¢-B'-w” = 0, (1.2) 
where primes and Roman numerals denote differentiation with respect to s and 


e=ejjtekk, «>0, «| >0. 


In the i, j, kK coordinate system, € is the ratio of the bending moment about the k 
axis to the angle of inclination of the beam with respect to the (i, k)-plane. A similar 
interpretation holds for e, . The end s 1 is assumed to be free; therefore, 

B lew!’ . B 1 w’’)’ saad 0 (1.3) 


there. It is the purpose of this study to examine some of the properties of the eigenvalue 
problem defined by the differential equation (1.1) and boundary conditions (1.2) and 
(1.3). 

2. Positive definiteness of the Rayleigh quotient. The Rayleigh quotient, R(w 


D(w)/H(w), which corresponds to the eigenvalue problem defined by Eqs. (1.1), (1.2) 
and (1.3), has been discussed in [1]. It was found there that 
D(w) [ w’-B ew” + OLw’-w’ mAQw:-@O-w) ds + w’(0)-e -w’(0), 2:3) 
and 

H(w [ mAw-w ds. 2.2) 


A straightforward calculation shows that each of the terms can be written as a sum of 
squares, but it is not obvious that D(w) is itself positive definite since the last term is 
preceded by a minus sign. It is not reasonable, on physical grounds, to expect that this 
functional could ever be negative, for this would imply the possibility of purely imaginary 
frequencies. Nevertheless, it is necessary, from the mathematical point of view, to 
establish positive definiteness in order to apply the usual Rayleigh-Ritz procedures and 


comparison theorems 
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The proof depends on theorems concerning M-definite and N-definite eigenvalue 
problems given by Kamke [4] and Kestens [5]. These results, although actually stated 
for scalar equations, can be extended to vector equations of the type considered here. 
The major change requires the use of the integral representation of the vector solution 
in terms of a Green’s tensor rather than a scalar Green’s function. The resulting algebraic 
structure is the same as the standard situation and the proofs can be carried through by 
appropriate changes in notation. 

Let us assume that D(w) is not positive definite and there is a value of Q’ for which 
the corresponding eigenvalue ’ is less than zero. Since \” is a continuous function of 2 
and \* > 0 for 2’ = 0, there must be a value of 2’ > 0 for which \* = 0. 


Let w be that value of 2’ for which \” = 0. Thus w is an eigenvalue of the problem 
(B'-w’’)”’ — w(Lw’)’ — mAw @-w = 0, O< e585 (2.3) 

( , 1 ad _ 
|w = w’ — ¢-B'-w” = 0, s = 0, (2.4) 


[Bo -w’’ - (B 1 ew’)! = 0, @ = l. 


The Rayleigh quotient R corresponding to this problem is 


[ w’’-B'-w’’ ds + w'(0)-e '-w’(0) 
R = —_______—. (2.5) 
[ (mAw:-@-w — Lw’-w’) ds 


“0 


Since the numerator is positive definite and the denominator is of undetermined sign, 
the problem is of the type designated by Kamke [4: IT, III] as A/-definite. He has shown 
that maximum-minimum theorems of the Courant type and comparison theorems hold 
for this problem. An upper bound can be found for the positive eigenvalues w” through 


these theorems. Let 
M = max (EI, , El), OCee tt 
then direct computation shows that 
w”’-B-w"’ < Mw’’-w”. 
In addition, if 


~ 


uw = max mA, y = min mA, 0<s<bl, 


then 
Lis) < ull? — s°), 


l’urthermore if 
w= W,J + W.-K, 
then 
w:O-w = Wi, 


and 
i al 


[ mAw-@-wds > v | Wi ds. 


“0 70 
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oa 


Also, if a = max (e;", €)'), 


then 


Let w*’ be the nth positive eigenvalue corresponding to the Rayleigh quotient 
al 
| Mw’’-w’’ ds + aw’(0)-w’(0) 

: SREY SERRE (2.6) 





oe a= 
R . 


io ie 
. [ W? ds — iu | (2 — &)w'-w’ ds 
ry) 0 


0 


We then obtain from Kamke’s comparison theorems 


The differential equations and natural boundary conditions which arise from the 
Rayleigh quotient (2.6) can be found, by elementary variational techniques, to be 


MWY — yw*[(? — s)Wi)’ — w*9W, = 0, 


'W,(0) = aW4(0) — MW7’(0) = 0, (2.7) 
Wi) = Wi") = 0, 
and 
[MWiY — tuw**[(P — 8) Wi) = 0, 
+W.(0) = aWi(0) — MW2’(0) = 0, (2.8) 


Wi) = Wi") = 0. 
Let us first consider the simpler of the two eigenvalue problems, (2.8). Multiplying 
the differential equation by W, , integrating from 0 to /, and applying partial integration 
together with the boundary conditions yield 
al i 
—huw*? | (0 — 8)W2 ds = afWXO)l + | MW ads. 


70 ¥0 


Consequently, there are no positive eigenvalues w** for problem (2.8); and the assumed 
solution must be found, if it exists, from problem (2.7). 
Equations (2.7) can be analyzed by observing that they represent the transverse 
vibrations of a uniform bar rotating about an axis at s = 0 with a constant angular 
velocity w*(u)'”*. The bar is elastically restrained at s = 0 and free at s = l. The circular 
frequency of transverse vibration corresponds to w*(v)'’*. According to the Southwell [6] 
inequality, the smallest eigenvalue w**y must satisfy 

w**y > K + w*’z, 
where K > 0 (with one possible exception to be considered later). Consequently, 

w**y > w**p. (2.9) 

However, 0 < v < u, contradicting (2.9) unless w** < 0. Thus there are no eigenvalues 


w** > 0 and hence no eigenvalues w > 0. Therefore in our original problem there is 
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no value of 2’ > 0 for which \* = 0, and all of the eigenvalues of the problem defined 
by Eqs. (1.1), (1.2) and (1.3) are positive. 
With this result we can now show that for any admissible vector w, as defined in [1], 


F(w) = i 


and hence complete the proof that D(w) is positive definite for any admissible vector. 
If (2.10) were false, there would be a vector u such that 


p. 255, 
I 


1 
mAw:®@-w ds — / Lw’-w’ ds < 0, (2.10) 
i) 


l 
[ u’’-B™'-u’’ ds + u’(0)-e '-u’(0) 
p= — — —-——— >0. 
i mAu:-@-u ds — | Lu’-u’ ds 
0 0 


A result of Kamke’s [4, III, p. 257] states that there must then be an eigenvalue w’ of 
Eqs. (2.3), (2.4), and (2.5) satisfying 


0O<w< p. 


However, we have just shown that such an eigenvalue cannot exist and hence F(w) < 0, 
proving that D(w) > 0 for all admissible functions. 

The only exception to the preceding argument occurs when the constant K in the 
Southwell inequality vanishes and u = v. This takes place when the bar is simply sup- 
ported, w(0) w’’(0) = 0, and the maximum and minimum values of mA are equal, 


i.e., mA constant. In this case, 
w=sJ 


is non-trivial and yields an eigenvalue \* = O for any value of Q*. Noting the relation 
between w and the displacement u, we see that this solution corresponds to the beam’s 
remaining straight, always lying in the plane of rotation, and making a constant angle 
with the i-direction. However, the non-negative character of D(w) is preserved if we 
regard this situation as the limiting case which occurs when ¢;', ¢;' approach zero. 

3. Difference equations and continuity conditions. The preceding section, together 
with the results contained in [1], shows that the lowest eigenvalue of the twisted, rotating 
beam is described by the minimum principle 


\” = miny, D(w)/H(w), 


where the class W, of admissible vectors w consists of all vectors w having continuous 
third derivatives, piecewise continuous fourth derivatives, and satisfying w(0) = 0. 
In many instances, useful comparison results could be obtained if these continuity 
requirements were lightened. The proof of such a minimum principle, however, would 
normally require all the restrictions given above. Fortunately, Courant [7] has shown, 
by a much more delicate analysis, that these restrictions for the admissible functions 
can be reduced without disturbing any of the requisite properties of the minimizing 
function. His proof goes much deeper than that usually given since it includes an existence 
theorem as well as an analysis of the convergence of a related set of difference equations. 
Although his paper is concerned with the standard, second-order, Sturm-Liouville 
problem, he states that the results will carry over for higher order, self-adjoint problems. 
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The extension to the present vector problem is straightforward, requiring only a few 
changes in small details, and none in the basic method of proof. The algebra and notation 
are, however, considerably more awkward. 

We shall therefore give only a brief discussion of the following theorem. Let the 
class W of admissible vectors contain those vectors w which are continuous, have con- 


tinuous first derivatives, piecewise continuous second derivatives, and satisfy w = 0 
at s = 0. Then among the vectors in W, there exists a vector, w"’, which renders D(w) 
a minimum subject to H(w) = 1. This vector is the first eigenvector of the eigenvalue 
problem (1.1), (1.2) and (1.3); the first eigenvalue is \‘’° = D(w)"?). If 
H(w,u | mAw-u ds, 

then, more generally, there exists a vector w'’’, belonging to W. which minimizes D(w) 
subject to H(w) = 1, H(w), w'’’) = 0,2 1, --- ,& — 1. This vector is the kth eigen- 
vector of the system and A’ = D(w'’) is the kth eigenvalue. 

The integrals D(w) and H(w) are replaced by finite sums obtained by dividing the 
interval 0 < s < / into n equal parts of length c, . Let w = w,j + w.k and r = d@/ds. 
Furthermore, let wo, -:* , Wi», *** , Wi, represent the values of w, at the points of 
division and similarly for w, . Then w, ,, w, (0) and w,,, = w,(/). The ordinary difference 


operators will be denoted by A, A’, --- , i.e., 


l 
Aw,,, = — [Wi441 — W,, 
C. 
I ] 
A*w,,, = 3 [u — 2w,,+u | 
Ca 
Let 
P,(w) = A’w,,, — W,,,7, — 27,Aw.,, — w2.,Ar, , 


Q,(w) = Aw, — UW, t+ 27,Au,., + w;.,A7, , 


R,(w) = Aw, , — Wo >T 

N,(w) = Aw, , — WU, >7, 

G,(w) = wij,, cos’ 6, — 2w,,,W2,, sin 6, cos 6, + w3., sin’ 8, 
= [M,(w)]’, 

T(w) = « ‘IR (w) |” + «¢, ‘IN (w)}’, 

Jw) =w,+w 


L,(s) =. i m.A ke. 


D,(w) = c, >> {ET,,,P°(w) + Ele.,Q2w 
+ OL,[Ri(w) + Ni(w)] — m,.,A,4:9°G,.:(w)} + 7(w), 


and 
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H,(w) = c, >, m,A,J,(w). 


Our original minimum problem can now be replaced by an ordinary algebraic eigenvalue 
problem by asking for the vector w, which makes D,(w)/H,(w) a minimum. This yields, 
for example, a lowest eigenvalue \{"’’ and a corresponding eigenvector w;”. The eigen- 
vector w.”’ is, of course, defined only at discrete points along 0 < s < 1. However, we 
may extend w;" to a vector w"’(s) defined over the entire interval by means of inter- 
polation. Unlike the second order case, parabolic rather than linear interpolation is used. 

It can then be shown that, as n becomes large, a subsequence of the \{”" converges 
to a limit \‘” and a corresponding subsequence of w "’(s) converges to a limiting vector 
w(s). Furthermore, w(s) has all the requisite continuity properties for the solution of 
the original differential equation and is an eigenvector with \‘”” the corresponding 
eigenvalue. The proof for the higher eigenvalues follows similar lines. 

As noted before, this modification of Courant’s proof yields an existence theorem, 
shows that the continuity requirements on the class of admissible vectors can be weakened, 
and finally proves convergence (in the sense of subsequences) for the difference equations 
which result from the minimization of D,(w)/H,(w). These equations are the algebraic 
Euler equations corresponding to 

\<?* = min D,(w)/H,(w), 


and are found to be 


QAlET, ,7,Q,(we)] + Es 4: At, 1Q,..(we”) — OQ A[L,R,(wy”’)] 


27 Lipari Tra Noa (Wo) — M,4;Ara12 cos 0,.,M,.,(we) — AN? m,.:4,4,0;"., = 0, 


th the natural boundary conditions 
Ptw,") = Ow,” ) = 0, 
A[EI, ..P,(w.”)] = ALEI,,,.Q,(w,”’)] = 0, 
Aw,"o — oI, oP.(w, ) = 0, 
Aws") — ek 1. .Qo(w,”’) 0. 


As a direct consequence of these results, we obtain Courant’s maximum-minimum 
principle. Let v'"’, 7 = 1, --- ,n — 1, be a set of independent vectors and let W be the 
class of vectors w which have a continuous first derivative, a piecewise continuous 


second derivative, and satisfy 


H(w,v'’) = 0, 1=1,--- ,n-—1. 


M(v'?) = glbyD(w)/H(w), 
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where glb» stands for the greatest lower bound with respect to all vectors w in W, then 


\”” = max M(y""’). 


4. Some applications of the minimum principle. Two simple applications of the 
comparison theorems implied by the Courant maximum-minimum principle will be 
mentioned here. Until now we have considered the case in which the fixed end of the 
beam occurs at the axis of rotation s = 0. Let us now assume that the beam is mounted 
on a finite hub of radius d so that the flexible portion is in the range d < s < Il. Then 


the differential equation (1.1) holds over the range d < s < l rather than 0 < s < / 
and the boundary conditions (1.2) are applied at s = d. 
Consider two beams identical save for the fact that one has a hub radius d, , whereas 


s 


the other has a hub radius d, (d, < d,). In the first case s has the range d, < s < l, 
whereas d. < s < / in the second. We shall now establish the inequality 


2 


»™* a) < ax™* @,), (4.1) 


a result which is physically reasonable. 

We shall define D,(w) in the same fashion as D(w) in Eq. (2.1) except that the range 
of integration now runs from d, to / and the boundary term is evaluated at s = d, 
We define H,(w), D.(w) and H,.(w) correspondingly. Furthermore, let W(d.) contain 
those vectors w which, in the interval d. <‘s < l, have a continuous first derivative, 


a piecewise continuous second derivative, and satisfy w(d.) = 0. We shall normalize 
these vectors by requiring that H.(w) = 1. The maximum-minimum principle thus 


states that if 
Mv’) = glbwiu,) D.(w), 


under the condition H.,(w, v'’) = 0,7 = 1,°--- , n — 1, then 


A" (d,) = max M,(v'"’). 


Vie 


Similarly, let W(d;) contain all vectors w in d, < s < I which have continuous first 
derivatives, piecewise continuous second derivatives, and satisfy w(d,) = 0. If 


Mi(v"’) = glbwra,) D,(w)/H,(w) 
under the condition H,(w, v'’’) = 0,7 = 1, --- ,n — 1, then 
“(d,) = max Mv""’). 
vi 
Let W* be the set of vectors w* defined over d, < s < 1 which are found by continuing 


each vector belonging to W(d,) into d, < s < d, under the following restrictions. When 
s = d, , w* = 0; w* has a continuous first derivative and a piecewise continuous second 


derivative in d, < s < 1; and D,(w*) — D.(w*) < p, where p is an arbitrarily small 
number. Finally, if for a given set of vectors v'"’,7 = 1, --- ,n — 1, a vector w in W(d,) 
satisfies H,(w, v'*’) = 0, then the corresponding continuation w* satisfies 

| mAw*-v'"’ ds = 0, ~=1,°-- ,n-1; 
1.e., H,(w*, v*’) = 0. Thus, W* is included in W(d,). 
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With these definitions, we see that for any set of vectors v for which H,(w, v"’) = 0, 


n 
1 D,(w) < all D(w*) 
g dw di) H,(w) ant q we H,(w*) 
In addition, I7,(w*) > H,(w*) = 1. Therefore, if H.(w, v*’) = 0,7 = 1, o¥e n-—l, 
' D,(w*) 
glby- H (w*) < glby.D.(w*) + p = glbw.s,)D.(w) + p. 


This result may also be written as 

Mv’) < M,(v"’) + p. 
\laximizing over v'’ and noting that p is arbitrary, we have (4.1). Inequality (4.1) 
agrees with the quantitative results found for the uniform beam by Boyce [8] for zero 
angle of inclination to the plane of rotation and by Schilhansl [9] for a general angle of 
inclination. 

As a second example, we shall see how some information concerning the effect of 
twist can be obtained simply from the maximum-minimum principle. Consider a non- 
rotating, twisted beam, clamped at the end s = 0. The governing equations are Eq. 

1.1) with Q = 0, Eq. (1.2) with e = 0, and Eq. (1.3). We shall also assume that EJ, > EI, 
for 0 < s < l. The corresponding Rayleigh quotient FR is 


na st 
R= | w’’-B tw’ ds /' | mAw-w ds. 


rom the inequality FI, > EI, , we see by direct computation that 
EI.w"-w" <w"-B'-w" < Elw"’-w"’. 
This may also be written in terms of the unit vectors J and K as 
EI.(Wi” + Ws”) < w’-B' ew” < El,(W{” + Wy”). 


Consequently, if A{""", 7 = 1, 2, represents the nth eigenvalue corresponding to the 


i 


Rayleigh quotient R,; defined by 


R, [ EI (Wi? + Wy”) as / | mA(W; + W3) ds, 
70 7) 
we have from the comparison theorem 


2 


a aa” ae". 


The sealar differential equations which correspond to the Rayleigh quotients R, are 


. 
) 


(EI,Wit’)’’ — mAXNW, = 0, =a (4.2) 
W.(0) = Wi0) = Wi) = We") = 0, 


and similarly for R, with the subscript 1 replaced by 2. 

We note that Eqs. (4.2) represent the ordinary equations of transverse vibration of 
an untwisted beam of flexural rigidity HJ, and mass per unit length mA. It is important 
to note, however, that the beam is capable of transverse vibration in two directions, 
i’, or W. , and hence each eigenvalue is a double one (one corresponding to W, = 0 
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and one to W, = 0). Consequently, if we think of \;"* or \;"’ as the nth eigenvalue of 
an untwisted beam vibrating in only one transverse direction, the appropriate inequality 


becomes 


2 SA SN L, 3 & eer (4.3 


Iquations (4.3) thus give bounds for the frequencies of a twisted beam in terms of those 
of an untwisted beam. 

5. The 1-p resonance problem. The techniques used in Sect. 2 can be applied to 
yield more insight into the practical problem of 1-p resonance. This question, which 
is of interest to the aircraft propeller designer, is to determine whether there is a frequency 
of rotation which coincides with the lowest natural frequency of the twisted blade. If 
such a frequency occurs within the operating range of the propeller, undesirable resonance 
phenomena can take place. It is hoped that the following remarks will shed some light 
on the mechanism responsible for 1-p resonance. 

Lo and Renbarger [10] have set up the equations for a uniform, rotating bar whose 
plane of bending is inclined at an angle y to the plane of rotation. The governing differ- 
ential equation, in non-dimensional form, is 


wy — sa [(1 — a )w’ |’ — (8 — a’ sin y)w 0, <2 < i. 

Here a is the non-dimensional rotational frequency and 8 a similarly defined frequency 
of vibration. Boyce [11] has shown for the clamped bar that if a < 12.36, there is no 
angle of inclination y for which 1-p resonance will take place. On the other hand, if 
a’ > 15.2, there will always be a value of y for which one can have 1-p resonance. Similar 
results hold for other end conditions. In other words, if the beam rotates at a speed 
above a fixed limit, inclination to the plane of rotation will produce 1-p resonance. 

Now let us consider whether it is possible to have 1-p resonance in an untwisted beam 
vibrating perpendicularly to the plane of rotation. In particular, we ask whether we can 
find variable EJ and mA distributions which will permit this type of resonance. 

Equations (1.1), (1.2), and (1.3) reduce in this case to 


(EIw’’)’’ — Q[L(s)w’)’ = mAd*w, 0<s< l; 
u w’ — ekIw’’ = O, = 6: (5.1) 
EIw"’ = (EIw’’)’ = 0, s = l. 

For 1-p resonance \* = 2’, and the differential equation becomes 


(EIw’’)” = Q[L(s)w’)’ + mAYw 


with the boundary conditions of Eqs. (5.1). This defines an eigenvalue problem for ©,” 


and the corresponding Rayleigh quotient is 


| Elw'” ds + €'[w’(0)}° 
= — , 


al a 


| mAw ds — | 


#0 /0 


L(s)w’’ ds 


If there is a positive eigenvalue 2°, 1-p resonance will take place. 
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M max HJ, y = min mA, wm = max ml. 


Again, the extended comparison theorems state that the positive eigenvalues Q* are 
bounded above by the positive eigenvalues w which correspond to the Rayleigh principle 
al 


M | w’”? ds + € '[w’(0)}’ 


st st 


y | w ds — yu | (? — s*)w” ds 


Jo /0 


The resulting eigenvalue problem is 


Mw'* + day wl[(U? — s)w’)’ w vw, 
w(0) e 'w’(0) — Mw’’(0) = w"() = w’’(D = O. 


However, this system is identical with Eq. (2.7). We have previously found that there 
are no positive eigenvalues in this case unless the beam is simply supported and the 
mass distribution is uniform. It is easily seen from Eqs. (5.1) that there is always an 
eigenvalue Q° in this case and I-p resonance takes place. Furthermore, this result 


remains true for the simply-supported, twisted beam with uniform mass distribution 


as is seen from the solution 
Ww sK. 


Qur analysis has thus shown that, except for the simply-supported bar of uniform 
mass distribution, there is no LI or mA distribution which will produce 1-p resonance. 
In other words, if the section is not inclined to the plane of rotation, this phenomenon 


cannot take place 
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BOOK REVIEWS 
(Continued from p. 334) 


The molecular theory of solutions. By I. Prigogine, with the collaboration of A. Bellemans 
and V. Mathot. Interscience Publishers, Inc., New York, and North-Holland Pub- 
lishing Co., Amsterdam, 1957. xx + 448 pp. $13.25. 


It has long been re¢ 
even the most elaborate versions fail to allow for changes in lattice parameters with composition, and 
more important, neglect the effect of composition upon the free volumes of the constituent molecules 
Furthermore, no lattice model will yield an equation of state for the mixture. On the other hand, em- 
pirical relations such as the Benedict-Webb-Rubin equation of state, although lacking a proper theo- 
retical basis, have been extremely successful in predicting the properties of certain mixtures from the 


cognized that the lattice theories of liquid mixtures are inherently inadequate 


properties of the pure components. 

Professor Prigogine and his co-workers have combined the best features of both approaches into 
a treatment in which the mixture obeys the same reduced equation of state as the pure components 
but with composition-dependent reduced variables. The theoretical justification of this is given in 
terms of the cell model for liquids, the difficult calculation of the free volume being avoided through 
the hypothesis that the reduced free volume in the mixture is equal to that of the pure components 
at the same reduced pressure and temperature. Although calculations made by this method seem some- 
what more cumbersome than those involved in the simpler lattice theories which are most commonly 
used to discuss solution thermodynamics, there is no doubt that the results obtained are in considerably 
better agreement with experiment. 

In addition to the detailed exposition and application of his own work on mixtures, Prof. Prigogine 
has given an excellent review of alternative approaches, including one dimensional solutions and the 
concept of conformal mixtures introduced by Longuet-Higgins. Quantum effects, especially in isotopic 
mixtures, are also treated. 

This is probably the most original and stimulating discussion of the statistical mechanics of solu- 
tions to appear within the last several years. It should be useful both to those doing research in the 
area and to those who wish to obtain a general acquaintance with the recent literature. 


STEPHEN PRAGER 


Ordinary difference-differential equations. By Edmund Pinney. University of California 
Press, Berkeley and Los Angeles, 1958. xii 262 pp. $5.00. 


Mr. Pinney has written a detailed account of practical methods of obtaining quantitative infor- 
mation on the solutions of difference-differential equations. He has directed his book towards the reader 
interested in applications, and the care with which the methods are described, and numerous illustra- 
tions make the book a valuable reference. Although nonlinear equations are discussed in the latter part 
of the book, these equations are nearly linear, and the book is primarily a study of a highly specialized 
linear analysis applicable to linear systems with constant coefficients. 

Integral transform and integral representation methods are described and illustrated in Chapter I. 
The Euler-Laplace transform is used in Chapter II to obtain infinite series solutions of a general system 
of linear difference-differential equations with constant coefficients. The solution of the characteristic 
equation is studied in Chapter III. Nyquist’s stability criterion is properly called Cauchy’s index 
theorem. A great number of specialized first and second order linear equations with constant coefficients 
are studied in detail in Chapters IV through VIII. Chapter IX is described as the “first practical theory” 
of nonlinear difference-differential equations. A method is given for obtaining approximate solutions 
of a nearly linear difference-differential equation. It is a method of successive approximation that begins 
with a solution of the nearby linear system. The method is illustrated in Chapters X and XI. Existence 
and uniqueness theorems for a class of linear integro-differential systems are given in an appendix. 


J. P. LASALLE 
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THE BLASIUS EQUATION WITH THREE-POINT BOUNDARY CONDITIONS* 


BY 
L. G. NAPOLITANO 


Polytechnic Institute of Brooklyn 


Abstract. The Blasius equation subject to three-point boundary conditions, de- 
scribing the interaction between two parallel streams, is solved by way of a series in 
terms of ascending powers of the ratio \ = (u, — u2)/u, , where the u,’s are the outer 
streams’ velocities. 

The first three terms of the series are analytically expressed in terms of the repeated 
integrals of the complementary error function (2” erfe ») and of the repeated integrals 
of the square of the successive integrals of the complementary error function (j"2" erfe 7). 
These functions often appear in problems leading to extended heat-conduction type of 
equations. A recurrence formula for j"7" erfe 7 is established and formulae relating the 
functions 7” erfe (—7) and j”j" erfe (--n) to available tabulated values of the functions 
2” erfc (n) are derived. 

The first three approximations to the Blasius function and to its first two derivatives 
are also presented in tabulated form with four significant figures. Test on the convergence 
of the series has been made by comparison with some exact solutions obtained by high 
speed computing machine. The comparison, extended to the physically essential quan- 
tities, shows that: 

(1) The Blasius function itself is slightly less accurate than its second and first 
derivatives. 

Two terms of the series for \ up to 0.5 and three terms for \ up to 0.7 yield 
extremely accurate results. The errors in the first two derivatives of the Blasius 
functions are always contained within less than one per cent. 

1. Introduction. The solution of the Blasius equation with three-point boundary 
conditions has per se a considerable academic interest. The availability of closed form 
solution has, however, become a practical necessity in view of the recent findings which 
have shown the essential and unique role played by this equation in isobaric mixing 
flows. It can indeed be said that all the types of plane two-dimensional interactions 
between two streams are governed by the Blasius equation with three-point boundary 
conditions. 

The reducibility of the basic Prandtl equations to the Blasius equation with pertinent 
boundary conditions was first shown by Gortler [1]** for the laminar incompressible 
mixing of uniform streams and, subsequently [2], for the turbulent case also. The present 
author showed [3] that the same happens for the compressible laminar case. A recent 
investigation by the present author [4] has also brought forth some evidence of an 
empirical correlation existing between turbulent and laminar compressible mixing. 
It was found that, under the assumption of a unitary turbulent Prandtl number, the 
velocity profiles are considerably independent of Mach numbers and density ratio and 

*Received December 19, 1957. This work was sponsored in part by the NACA under Contract 
NAw 6480. Some of the results have appeared in Ref. [9] which was issued in a limited number of copies 


not for circulation. 
**Numbers in brackets refer to the bibliography at the end of the paper. 
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can therefore be deduced from the solution of the Blasius equation. This statement has 
an even larger implication insofar as it applies whenever the density ratio can be given 
a parabolic dependence on the velocity ratio. Thus the field of application of the Blasius 
function is widened to include a large variety of interactions between streams of different 
gases [5]. Solutions of the Blasius equation are furthermore needed in problems of laminar 
and turbulent mixing of non-uniform constant vorticity streams. These problems are 
solved by means of a series solution in terms of the “vorticity numbers’’: the zeroth 
order terms is the Blasius function and the coefficient of the equations for the higher 
order terms are all functions of the Blasius function and its derivatives [6]. 

The present paper is mainly concerned with the solution of the Blasius equation and 
not with its derivation for which reference is made to the pertinent literature. The solu- 
tion is obtained as a series in terms of the parameter \ = (wu; — U2)/u, . The first three 
terms (up to \*) are given in explicit closed form and in tabulated forms. Owing to the 
complicated nature of the terms of the series, its convergence could not be formally 
established. The results of the present method are, however, compared with some exact 
solutions obtained by high speed computing machine calculations. 

In the course of the mathematical treatment there often appeared successive repeated 
integrals of the complementary error function (symbolically indicated by 7 erfe 7) 
and successive repeated integrals of the square of the functions 2” erfe n. As this feature 
is common to a large variety of physical problems which can be reduced to extended heat 
conduction type of equations, a summary study of those functions is presented in Appen- 
dixes A and B. Formulae necessary to compute the values of 7” erfe (—7) in terms of 


the already tabulated values (up to n = 11) of 2 erfe 9 are developed. A recurrence 
formula is established for the functions j"2" erfe 7 defined as the successive integrals 
of the functions (2” erfe n)*. Finally relationships giving the functions j"7" erfe 7 in 


terms of the functions 7” erfe y are derived which afford a rapid evaluation of the func- 
tions themselves. 
This work is part of a program of investigation on mixing phenomena carried out at 
the Polytechnic Institute of Brooklyn under the supervision of Prof. Antonio Ferri. 
2. Solution of the Blasius equation. The equation to be solved is 


f’"’ + off’ =0 (1) 
and it is subject to the following three-point boundary conditions 


lim f’ l lim f’=1-—A f(O) 0 <= K < 8. (2) 


§ =) 


Quantities related to the mixing of two streams are expressible in terms of the function 


f(¢) and its derivatives as follows 


y = 2m, x) “*f(d), 


: —“_ « 
f= sy ) +k, 
My 


“u= uf’(o), 


1/2 
Uy P op . 
v= (“) (ge — kf’ — ff]. 
* : 


3) 


1959] BLASIUS EQUATION WITH THREE-POINT BOUNDARY CONDITIONS 399 

In these equations x and y are space coordinates whose origin is taken to be at the 
point where the interaction begins, u and v are the corresponding velocity components, 
y is the stream function defined by u = y, ;v = —y, and v is the kinematic viscosity coeffi- 
cient. The quantity k is an arbitrary constant whose presence follows from an interesting 
property of the Prandtl boundary layer equations. These equations are invariant under 


the transformation 
yi = y + 8(2), 
“1, = 2, 
u(x, y) = ula, , y), 


ds 


v(x, y) = v,(z7, ,y:) — ula, , Wi) =-- 
dx 


Asymptotic boundary conditions on the 2-component of the velocity remain also 
unchanged while boundary conditions on the y-component are changed. The mathe- 
matical implication of this fact lies in the freedom of choosing arbitrarily the third 
boundary condition for the Blasius equation. The physical implication is the resulting 
indeterminacy of the wake orientation insofar as the transformation back into the physical 
plane cannot be performed unless k is known. Additional physical considerations, such 
as the one suggested by von KaArman that a free wake be acted upon by a zero resultant 
force in the y-direction, will uniquely determine this constant k and thus will fix the 
orientation of the wake. The solution herein presented relates to a wake whose streamline 
through the origin satisfies the equation y = —2k(vx/u,)'” 

The following series solution of Eq. (1) in ascending powers of \ is sought 


f= Dar‘, . (5) 


ISquation (5) is substituted into Eq. (1) and the coefficients of the successive powers 
of X are set equal to zero. The zeroth order approximation must satisfy the following 


equation 


fo’’ + 2fofs’ = 0 (6) 
with 
lim fg = 1, 
4 ‘=a (7) 
{.(0) = 0. 
The pertinent solution is fy = ¢ and it corresponds, physically, to zero mixing. The 


equation for the first approximation is, by taking the zeroth order solution into account 
fi’ + 2¢fr' = 0 (8) 


with 


lim ff = 0, lim f; = —1, 
[-+@ f--« (9) 
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The solution of Eq. (8) is 


a$ rt 


fi = (x) i | dt | e” db — Lt. (10) 


Finally, the 7th approximation (¢ > 2) must satisfy the equation 


fi" + 2th? = f° RO, (11) 
where the #;(¢) are functions, at most, of the (¢ — 1)th solution and are given by 
RAQ) = —2 ¥ facaftea/fl!- (12) 
A ) 
Equation (11) is subject to the boundary conditions 
lim f/(0) = 0, Fie) = 6. a>? (13) 
and admits the solution 
f-=@)''C,.f -@” | dt | e db | R,(a) da (14) 
with 
C, F = (a) | é - db | R (a) da. ( 15) 


Successive approximations to the solution f can thus be readily determined to any order. 

The task is considerably simplified if Eq. (14) is expressed in terms of the comple- 
mentary error function and of related functions. In this problem, indeed, as well as in 
several other problems leading to extended heat-conduction type of equation, the solution 
can be expressed rather simply in terms of the following functions 


" erfe n = | | cee | (erfe n,) dne +++ dn, , 
erfe 7 = | | cee | ” erfe n, g dn. *** dn 
T ani | 2” erfe t)-(2" erfe t) dl. 


Of these functions only the first ones, usually referred to as repeated integrals of the 
complementary error function, have been studied. Hartree [7] has shown some of their 
properties and applications and, more recently, Kaye [8] has tabulated them, up to the 
eleventh repeated integral, for positive value of the argument. 

The functions 2” erfe (—7) are considered in Appendix A wherein their expressions 
in terms of the functions 2” erfe 7 and their asymptotic behavior are presented. A sum- 
mary study of the functions j"2” erfe n is given in Appendix B. Therein the existence of 
a recurrence formula is proved and expressions relating the functions j"2" erfe (+n) 
and T’,,,(-7) to the repeated integrals of the error function are given. 

If simplified notation such as 


t” erfe n = 0", 


jt erfen =j'r’, 
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are adopted, the first approximations to the function f can be given simple analytical 
expressions as follows: 


First approximation— 


f, = 3[t — (x) "I, 
ff = —}hi = —}erfc f, 
ft! = (x)? e*,” 


Second approximation— 


5 xs mb itinee dite At ee OE ed 
| (167) [3 rj iP | + r E (x)' 2 ty, 


0 om} 
2 2 1-0-2 


$  4(n)'? 
fs T em) '?e —44+ 7). 


TABLE 1 


Blasius equation with three-point boundary conditions—V alues of fj 


r h(¢ f(-$) fof) fol —$) fo) fs —$) 
0.00 .0000 0000 .0000 .0000 .0000 0000 
0.01 — .0050 0050 0004 — .0006 0002 — .0002 
0.02 — .0099 0101 0009 — .0009 0004 — 0004 
0.04 — .0195 0204 0017 — .0032 .0009 — .0009 
0.06 .0290 .0310 .0025 — .0034 .0013 — .0013 
0.08 0382 0418 0032 — .0041 0017 — .OO18 
0.10 — .0472 0528 .0038 — .0052 .0021 — .0023 
0.12 — .0559 0640 0044 — .0065 .0025 — .0028 
0.14 — .0645 0755 .0050 — .0077 .0029 — .0933 
0.16 0728 .0872 .0054 — .0090 .0033 — .0038 
0.18 - .0809 .0991 .0058 — .0104 0037 — .0043 
0.20 0888 .1112 .0063 O118 0041 — .0048 
0.30 1250 .1750 .0078 — .0196 0058 — .QO75 
0.40 1560 2440 0075 — .0295 0074 — .0104 
0.50 — . 1823 3177 0069 — .0381 0085 — .0136 
0.60 2041 3959 0056 — .0503 .0098 — .0171 
0.70 2220 1779 .0037 — .0610 .0107 — .0209 
0.80 — .2365 5635 .0016 — .0705 01138 — .0250 
0.90 - .2480 .6520 — .0005 — .0778 0018 — .0295 
1.00 2570 7430 — .0024 — .0867 0120 — .0342 
1.20 — .2691 .9309 — .0060 — .1017 .0122 — .0437 
1.40 2758 1.1243 — .0085 — .1126 0121 — .0529 
1.60 2792 1.3208 — .0102 — .1199 0119 — .0609 
1.80 2809 1.5191 — .0112 — .1238 0117 — .0674 
2.00 — .2816 1.7184 — .0117 — .1268 0116 — .0720 
2.20 — .2819 1.9181 — .0119 — .1280 0115 — .0750 
2.40 2820 2.1189 — 0120 — .1286 .0115 — .0765 
2.60 2821 2.3179 — .0121 — .1288 .0115 — .0771 
2.80 — .2821 2.5179 — .0121 — .1289 0115 — .0773 
3.00 2821 2.7179 — .0121 — .1289 .0115 — .0773 
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Third approximation— 
f= PF, +F(-©)fi ; 
. fit. f=} 3/1 | aa 
po = —- |= 4+ 2 — 4+ 0 +5\- ie 2" 
' 2 4 2Qr 2 \4r 16 J a 
——— : a "te 
+ —(2( 7 — (9)"] — = T Gran) 
7 (2¢ 1 - qa + Fen? 
| 1 —— = eo 
cata -— Gy ~*+eF + - [i(d)]° dt, 
2 Lt T 2 4 J; 
> ra l e —_ 1 (a) 5 4 e 
 — GY = i aS S735 bal , = c a > ' 
Se [i+ ; * Fe 
, : — a 
2 1/2 -2¢ 2° 
— (2 — 1) —@ “tt -—5I*tT 
TABLE 2 
Blasius equation with three-point boundary conditions—V alues of f; 
fi( fi(—-$) fo( ¢) fx —$) {X¢) {x -5) 
0.0 — .5000 — .5000 0454 0454 .0222 .0222 
01 — ,4943 ~— 5056 0440 0468 .0220 .0224 
02 — ,4887 — .5113 0426 0482 .0218 0225 
04 — 4774 — .5225 .0398 0510 .0215 0229 
06 — 4662 — .5338 0370 0538 0211 .0233 
08 — .4550 — 5450 0342 0566 0208 .0236 
10 — .4438 — .5562 0314 0594 0204 0240 
12 — .4326 5674 0287 0621 .0201 0244 
.14 — .4215 - 5785 0260 0648 .0197 0247 
16 — .4105 — .5895 0234 0674 £0193 0251 
i8 — .3995 — .6005 0218 0700 .0190 .0255 
20 — .3886 — .6113 0183 0725 0186 .0259 
.30 - 3357 — .6643 0080 0839 .0166 0280 
410 — .2858 7142 — .0033 0929 .0146 0304 
50 — .2397 ~ .7602 - 0109 0990 0123 0332 
60 1981 8019 - 0163 .1018 0098 0364 
.70 ao S613 — .8389 0195 1014 0075 0399 
80 1289 — .8710 0210 0980 0052 0433 
.90 — 1015 - 8984 0209 0920 .0034 .0450 
1.00 ~ 0786 — 9213 - 0197 .0841 0018 0475 
.20 — .0448 — .9551 - 0154 0648 .0000 0476 
40 — ,0238 — .9761 0105 0451 — .QO10 0439 
60 — .0118 —0 .9882 0064 0285 — .0012 .0369 
.80 — .0054 —0Q .9945 — 0035 .0123 — .0008 0278 
2.00 — .0023 —0.9977 — .0017 0086 — .0004 0010 
2.20 — 0009 —0.9991 0008 0041 — .0002 0107 
2.40 — .0003 —() .9996 0003 0018 — .0001 .0045 
2.60 — .0001 —Q) .9999 — .gool 0007 — .0001 0018 
2.80 0000 - 1.0000 0000 0003 .0000 0002 
3.00 .0000 - 1.0000 0000 .0001 .0000 0000 


——__— 
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TABLE 3 
Blasius equation with three-point boundary conditions—V alues of {*’ 
$ fi'(s) fi'(-$) f2'(¢) 2( —f) fs'(¢) f3'(—$) 
00 .5642 . 5642 — .1410 — .1410 — .0176 — .O176 
01 5641 .5641 1410 — .1410 — .0176 — .0176 
02 5640 .5640 1409 — .1409 — .0177 — .O177 
04 5633 .5633 1404 — .1404 — .0177 — 0177 
06 5622 5622 . 1395 — .1395 — .0178 — 0178 
O08 5606 5606 1384 — .1383 — .0179 — .0179 
10 5586 .5586 1370 — .1367 — .0180 — .0180 
12 5561 .5561 - .1352 — .1348 — .0180 — .0182 
14 5532 5532 — .1332 — .1326 — .0182 — OI184 
16 5499 .5499 — . 1309 — . 1300 — .O0186 — 0188 
18 5462 .5462 1283 — .1271 — .0188 — 0193 
20 5421 5421 1255 — .1239 — .0190 — .0200 
30 5156 5156 1083 — .1031 ~ .0201 — .0225 
40 1808 .4808 OS874 — .0760 — .0216 — .0255 
50 1394 4394 0650 — .0448 — .0235 — .0300 
.60 3936 3936 — .0430 — .0121 — .0265 — .0345 
70 3456 3456 0230 0195 — .0250 — .0350 
SO 2975 2975 0061 0478 - 0218 — .0320 
90 .2510 . 2510 0070 0708 — .0171 — .0205 
1.00 2075 . 2075 0165 O875 .O118 — .0100 
1.20 1337 .1337 0246 .1011 0072 0100 
1.40 0795 0795 0233 0927 0031 0260 
1.60 0436 0436 0176 0722 OO1LS 0435 
1.80 0221 .0221 0114 0491 0028 0480 
2.00 0103 0103 0065 0297 0018 0440 
2.20 0045 0045 0033 0160 0007 .0350 
2.40 0018 0018 OO15 0078 0004 0220 
2.60 0006 0006 0006 0035 0002 0105 
2.80 0002 .0002 0002 0014 0001 0050 
3.00 0001 0001 0001 .0005 0001 0008 


The first three approximations to the functions f(¢), f’(¢) and f’’(¢) have been com- 
puted to six significant figures. Values of f;(¢), f/(¢) and f/’(¢), (¢ = 1, 2, 3) are tabulated 
in Tables 1 and 3 with four significant figures. 

3. Accuracy of the solution. The convergence of the series given by Eq. (5) could 
not be formally established, owing to the complexity of the relative terms. Practical 
indications about the rapidity of the convergence, however, are derived by comparing 
the results with those of some exact solutions obtained with the D12 Differential Analyzer 
presently in operation at the Centro di Calcolo Elettronico of the University of Naples.’ 

The following indicative values for \ were considered: \ = 0.2678; \ = 0.4796; 
\ = 0.5541; = 0.6915. Comparison was extended to the following physically meaningful 
quantities: 

i) f’(0), (proportional to the component u of the velocity along the £-axis); 


'These calculations are part of a larger program of high speed machine solution of turbulent mixing 
flows sponsored by the United States Air Force through the Air Force Office of Scientific Research, 
Air Research and Development Command, under Contract AF 18(600)-693, Project No. 17500. The 
cooperation of Prof. Giorgio Savastano, Associate Director of the C. C. E., is gratefully acknowledged. 
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ii) f’’(0), (proportional to the shear stress along the x-axis); 
iii) lim... [¢f’ — f], (proportional, for k 


at the edges of the wake). 
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0, to the y-component of the velocity 


Values obtained from the exact solutions and from the first, second and third approxi- 


mations to the function f(¢) are listed in Table 4. 


h = 0.26783 


f’(0) 

f’’(0) 

lim [¢/’ — f] 
[+o 


lim [¢/’ — f] 


= 0.47961 
f 


‘'(0) 
f’’(0) 
lim [¢f’ — f 
lim [¢f’ — f] 


f’(0) 


{’'(O 
lim [¢f’ — f] 
lim [¢f’ — f] 


f+o 


; 
\ = 0.69147 


/’(0) 
f’"(0 
lim [¢f’ — f] 
c+o 
lim [¢f’ — f] 


Comparison between exact and approximate solutions 


Exact 


0.8698 
0.1406 
0.0762 


0.0865 


0.7421 
0.2660 
0.1568 


0.2154 


0.6870 
0.3158 
0.1937 


0.3003 


TABLE 4 


I Approx. 


0.8661 
0.1511 
0.0755 


0.0755 


0.7602 
0.2706 
0.1353 


0.1353 
0.7229 
0.3126 
0.1563 
0.1563 
0.6543 
0.3901 
0.1951 


0.1951 


The following comments are proper: 


II Approx. 


0.8693 
0.1410 
0.0764 


0.0847 


0.7706 
0.2382 
0.1381 


0.1649 
0.7368 


0.2693 
. 1600 


~ 


0.1959 


0.6760 
0.3227 


rare | 


0.2009 


0.2567 


Approx. 


.8697 
. 1407 


0762 


0862 


. 7406 


2663 


. 1580 


2090 


6833 
3169 
1971 


2822 


(1) The best agreement is obtained for f”, followed, in order, by f’ and f. 


— 


| —¢ | large. 


(2) The accuracy decreases with \ and, for a given X, is greater for ¢ large than for 


(3) The y-component of the velocity at the lower edge of the wake always exhibits 


the maximum percental error. 


(4) The first two approximations are more than satisfactory up to values of 
The errors in f’(0) and f’’(0) 
(5) Three terms of the series are 


are less than one per cent. 
needed for values of \ greater than 0.5. With three 


= ().5. 


terms the errors are nearly always contained within less than one per cent up 
to \ equal to 0.7. The only exception lies again in the value of the y-component 
of the velocity at the lower edge of the wake. The approximate value is 6% 


smaller than the exact one. 
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APPENDIX A 


The error function and its repeated integrals. The error function is defined by 


2 a 
erf 7 = 773 e dt (Al) 
(7) " Jo 
and the complement of the error function, erfe 7, by 
2 “ t? s <)) 
erfe 7 = 1 — erf yn = ain | e at. (A2) 
(7) ; 


The nth repeated integral of erfe 7 is symbolically defined as 


@ 


t" erfe 7 = [ (i""* erfe t) dt (n > 1) (A3) 


with 
i’ erfe n = erfe n. (A4) 


The functions given by Eq. (A3) are tabulated, up to n = 11, in Ref. [8] for positive 
values of the variable 7. 
It is of interest to extend the definition of 7” erfe » to negative values of 7 as follows 


ny 


Pale(—~ « [ (i"" erfe ¢) dt, (A5) 


and to see whether it is possible to express them as functions of the known tabulated 


values of 2” erfe 7. 
In consistence with the notation given in Eq. (A3), the first derivative of the error 


function is defined by 
d —} 2 —7? { 
——(erfc 7) = 7 erfe n = —we". (A6) 
dn (7r) 
It is, apparently, an even function, so that 


i’ erfe (—n) = t"' erfe n, (A7) 
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whereas 
0 


i’ erfe (—n) = 2 — 7’ erfe 7». (A8) 
By repeated application of the following recurrence formula [7], valid for n > 1 
2ni” erfe n = i" ” erfe n — 2mi"' erfe 7 (A9) 


it is then easily verified that the following general identity holds 


— ee ae ‘ ht Qs i ed Eee ee 
eric (— 7) (—1)" 3 eric 9 + 2 - E ae Fe n'i"~" erfe (0). 
- Bs 


h=0 


A10) 


In deriving Eq. (A10) it has been taken into account that Eq. (A9) yields, for » = 0 
2ni” erfc (0) = i” ~ erfe (0) 
or 


, l 
2 erie (UY) = = . 
2"(4n)! 


Equation (A10) gives the required relation between the repeated integrals of the error 
function for negative and positive values of the independent variables. Thus, for instance, 
the first three integrals in the negative range of y are simply expressed by 


i erfe (—n) = 2n + 7 erfe n, 
i’ erfe (—n) = —? erfen + n° + 27° erfe (0), (A11) 
i* erfe (—n) = 7° erfe n + 2m’ erfe (0) + 7°/3. 


Apparently these repeated integrals will not converge as » — — ©. Their asymptotic 


behavior, for | —7 | large, is 
— oe Ma ee SPT \ wae ST 
i" erfe (—n) ~ (—1)""* —a =i + 2 — - "~" erfe (0). (A12) 
erf¢ n in)? (Qn)" ° ? 5 hi erfe ( \] 
APPENDIX B 
Repeated integrals of the functions (7” erfe »)°. Let the successive integrals of the 


functions (z" erfe 7)” be symbolically indicated by 


a4" erfc 7 = | ie a” erfe t) dt m > 0) (BI 


with 
ji” erfc » = (2” erfe 7)’. (B2) 


It is desired to express these functions in terms of the repeated integrals of the error 


function. 

The relationship is immediate for the two particular cases: n = —1 (m any positive 
integers) and m = 1 (n > OQ). 

Indeed when m = 1 and n = —1, it is by definition 


jt erfe n = | | a rad | dt, (B3) 
Ja T 
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>\? 


so that [see Eq. (A3)] 
rf a 1/2 
jt erfe n + (: erfe [n(2)°’"]. (B4) 
\ 


Repeated integrations easily yield the required relation between j"7' erfe 4 and 
the repeated integrals of erfe y as 
“mm 1 (2) a *m—1 \1/2 - 
yj". erfe 7 = 7st erfe [n(2)'’*]. (B5) 
(xr) 
A corresponding expression, valid for m = 1 and any n > 0 can be obtained by 
repeated integrations by parts. As it is easy to verify, the following identity will result 


2in — h + 1)]n!2""" 
a nar [2(n h- In!2 Pe need ae 
ju” erfe 7 2 —] On + Din —h pit erfe 7:2 erfe 7 


(B6) 
19' 2+n 
pm htl oe yoy | yn 1/2 eee, rgy!/2 
ni erfe n) | + (-—1)"" (nr) (Qn + py ere {(2)°“~n]. 
In deriving Eq. (B6), the following identity 
len d nm ‘a 2 n n—I - 
} tw’ erfen = — in (7" erfe n) = 22” erfc 7-2" erfe 7 (B7) 
ay 

which constitutes an obvious extension of the definition (Bl) to the case m = —1, 


has been taken into account. 

In the most general case use must be made of a recurrence formula. This formula 
can be derived by successive and repeated integration by parts of Eq. (B1). By taking 
Iq. (AY) into consideration one obtains 

1 *m -n—1 


(2n + m)j"i" erfe » = 3j" “7” erfe n — nj” @ erfe yn — j"t" erfe 7 (B8) 


valid form > l andn > 0. 

The recurrence formula (B8) together with Eqs. (B5) and (B6) afford a rapid com- 
putation of the functions j"2" erfe n. Their extensions to negative values of the argument 
are readily accomplished by means of Eq. (A10). 

The first few functions j”2” erfe n, are explicited and their values at 7 = 0 are given. 
A simplified notation such as 


2” erfe n = 1 (B9) 


“mon 


jt" erfe n = 7" 
is adopted. Thus 


9\1/2 
jv i 1°2 1 = n(i°)? = (2) 1°[n(2)'”7], 
Th 


sa in i sical : B10 

pi? = 0)? — Anji” — (2)-*i[9(2)””], ates 
ji = Bit — bald? — Bi, 
1:20 


ji = 3) — int -— dt, 
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and 
6 , an 2 
n°(0) = (x)*? [2 (2) °"] 47-2 (Q) l-—-- 
rs 
(B11) 
**7/)\ -1/279)! 27m) _ l 
371(0) = (7) 1(2) _ 1] 47° 2(Q) -— 7 
To conclude, the integral j 
PT nn(Q) = | i” erfe t-7" erfe t dl n> m) 
will be evaluated. 
Repeated integrations by parts and consideration of Eq. (B6) yield 
~ —] a ' , - , m+» * . TE 
Tani) = a 1(2  erfe n ir > (-1 |2 eric 7:2 ' erfe 7] (1312) 
when n — m = 2h + 1, and 
Tani?) = 1) "92 “ erfe » + ps (—1)""|i"*” erfe n-t" "*' erfe n] B13) 
when n — m = 2h. 
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—NOTES— 


EFFECTS OF GRAVITATIONAL OR ELECTROMAGNETIC 
FIELDS ON FLUID MOTION* 


By CHIA-SHUN YIH (University of Michigan) 


Summary. It is shown in this paper that the effect of gravity on a stratified fluid is 
to inhibit steady motion in the direction of gravity; that, for small values of the magnetic 
viscosity, the effect of a main magnetic field is to make steady weak motions of a fluid 
independent of the distance along the lines of force; and that, again for small values of 
the magnetic viscosity, the effect of a uniform electric field is to make steady weak 
motions of a fluid rotationally symmetric with respect to an axis in the direction of the 
field. These results, together with a similar one of Proudman (1916) [1j for a fluid with 
general rotation, enable one to state that the effects of rotation, gravity, and electro- 
magnetic fields (for small magnetic viscosity) are to endow the fluid with a certain 
anisotropic rigidity by “stiffening” it along the vorticity lines, the isopycnic surfaces 
or lines, or the lines of foree, as the case may be. In the case of a weak steady motion or 
relative motion, this “stiffening” has the effect of reducing the a priori number of physical 
or at least mathematical dimensions of the motion by one. 

1. Introduction. Since Proudman [1] showed that the effect of rotation on weak 
steady relative motion of an inviscid fluid is to make it two-dimensional, in the sense 
that the motion is independent of the distance along the axis of rotation, experiments 
performed by Taylor [2], Long [8], and others have largely supported his assertion. 
Recently, Gariél’s experiments with a stratified fluid and Lehnert’s with mercury have 
indicated similar effects of gravity and of a magnetic field, respectively. The purpose 
of this paper is to describe these effects, to predict the effect of an electric current on 
weak steady motions of a fluid, to explain all these effects mathematically and, as far 
as possible, to bring out their similarity to one another. The fluid is assumed to be 
inviscid and incompressible throughout and, in the case of an electromagnetic field, to 
have negligible magnetic viscosity. 

2. Effect of rotation. For convenience of comparison and for completeness, Proud- 
man’s demonstration [1] of the effect of rotation will be briefly presented. If the general 
rotation of the fluid is Q, the axis of rotation can be chosen to be the z-axis. Two other 
axes, both rotating with angular velocity ©, can be chosen to provide a rotating cartesian 
frame of reference, with coordinates (x, y, z). The velocity components relative to this 
rotating frame in the directions of x, y and z will be denoted by u, v, and w respectively. 


The equations of motion are then (Morgan [4]): 


Du _ _oP (1) 
Dt — ax’ 
Dv oP 
i  e.. 2 9 
Di 2Qu ay” (2) 
Dw oP 
- fe 9 = —-— 7 
Dt oa Oz’ (3) 


*Received August 2, 1957. 
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in which 


D ey 0 ¢ ; 0 
= _ a S pa w 
Dt at Ox oy Oz 
and 
) 9 2 
r - O(yY +2) +¢ } 
p 2 


with p, p, and @ as the pressure, density, and the body-force potential, respectively 
The equation of continuity is 

Ou, ow , dw a 
" — a —— as (, dD) 


Ox Oy Oz 


lor steady and weak relative motion, Eqs. (1) to (3) become 
oP , 
: 0, >) 
Ox 


oP 
= 20Ow 7 
oy 


oP 


Oz 


—20On 


By cross-differentiation of hqs. (7) and (8), it can be readily shown that 


Ov , Ow 
: a 0. )) 
oY O02 
It follows from hq 5) that 
Ou 
: Q). 10) 
O72 


Furthermore, by differentiation of qs. (7) and (8) with respect to 2 it ean readily 
b I 


seen with the aid of Ikq. (6) that 


fan Ow : 
: 1] 


Equations (10) and (11) indicate that the motion is independent of x. Hence, the steady 
and weak relative motion is two-dimensional. In a loose way, one might say that the 
general rotation has endowed the fluid with a certain rigidity along the vorticity lines. 

It can be readily shown that for symmetric motion the radial velocity component 
so that the streamlines are spirals wound around circular cylinders). This 
1 be directly compared with the vanishing of the vertical velocity com- 
weak motion of a stratified fluid to be treated in the following 


vanishes 
situation car 
ponent in the steady 
section, for in the case of rotation one may imagine a virtual gravity acting in the radial 
direction, according to Mach and Einstein 

3. Effect of gravity. One now turns to the effect of gravity on steady and weak 


motions of a stratified fluid. Let the z-axis be vertical. The equations of motion are 
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Du Op ‘ 
Dt = ax’ (12) 
Dv Op psiiale 
> iar. oe (1% 
P Dt ay’ 13) 
Dw Op ’ 
Phe a, (14) 


in which D/Dt, as usual, stands for the substantial differentiation. The equation of 


incompressibility is 


Dp 


= ( 4 
Dit 0 15) 


which permits the equation of continuity to be given by Eq. (5). 
The density can be expressed as the sum of a prevailing density p, , which depends 


only on z, and a perturbation density p’: 


p = p,(z) + p’. (16) 
For steady weak motion, Eqs. (12) and (138) become 
Vp Op as 
= @, ~~ «© (17) 
OX OY 


and hence, from qs. (14) and (16) one obtains 


, , 


) 0 
<! = 0, <a 
OX oY 


= (). (18) 


IXquations (17) and (18) are correct to the first order. But from Eqs. (15) and (16) it 


follows that 


y , y , d 4 , 
u = + y ss +w < +w e 0 (19) 
Ox OY dz Oz 
30 that, from Eqs. (18), 
w= 0 20) 


to the second order. Thus gravity inhibits steady velocity in the vertical direction, and 
if the disturbance is not strong enough, steady motion can occur in horizontal planes 
only. This situation is directly comparable to that for symmetric relative motion with 
general rotation, provided that that motion is steady and weak, as mentioned in the 
stiffened.” 


ee 


last section. In the present case, it is the isopycnic surfaces or lines that are 
Another result comparable to Proudman’s can also be obtained. If two-dimensional 
motion is considered to start with, so that v is zero, Eqs. (5) and (20) yield 


} 
— = @. (21) 
OX 


Thus the motion will be one-dimensional, in the sense that the state of affairs does not 
change with z—in much the same way that Proudman’s motion is independent of 2. 
In both cases x is measured in a direction normal to that of gravitation or virtual gravi- 
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tation created by rotation. The essential one-dimensional character of steady weak 
motions of a stratified fluid has been experimentally demonstrated by Gariél [5]. 

4. Effect of a magnetic field. One now considers a uniform magnetic field H, acting 
in the z-direction, and imagines a weak steady motion to take place in this field. The 
fluid is assumed to be homogeneous. The perturbation magnetic field will be denoted 
by h, so that the total magnetic field is 


H = H, +h. (22) 
The vector equation of motion is (all electromagnetic variables being measured in emu), 
with g as the body force per unit mass, 


Dv 


’ Di = —grad p + pg + uj X H (23) 


in which v is the velocity vector, u is the magnetic permeability and j the current density. 


Since 
curl H = 4nj, (24) 
Eq. (23) may be written 
= grad p’ (2) 1ix (am) 25) 
‘) — = —grad p grat —s +T Glv - \ao 
? Di ee ees ae 4a 


in which 

p’ = p+ po (26) 
with @ again denoting the potential of the body force. The intensity E of the electric 
field is related to the magnetic field by 
bie OH - 
curl E = —p —. (27) 
ot 


Furthermore, the current density, electric field, magnetic field, and velocity field are 


related by the equation 


j = o(E + uv X Hi) (28) 
in which ¢ is the electric conductivity. The fluid being homogeneous, the equation of 
continuity is again given by Eq. (5). 

Since the motion under consideration is weak and steady, Eq. (25) becomes, after 
linearization: 
Op’ 
I 0, (29) 
OX 
op’ _ A, (dhs 7 dhe), (30) 
Oy Sor \ Ox Oy 
Op’ H, (oh, a) 31) 
Oz 4r \ Ox oz] 
Irom Eqs. (29) to (31) one obtains by cross differentiation 
re) ‘ 
curlh = 0 (32) 


Cd 


JS 


e 


— 
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or, by virtue of Eq. (24), 

re) 

—j=0 33 

Ox ( ) 
so that the current density is independent of x. Although Eq. (33) is not needed to 
prove the independence of the fluid flow on z, it contains a result which is worthy of note. 


The linearized form of Eq. (28) is 
j. = oF, , (34) 
j, = o(E, + vow), (35) 
j. = o(E, — pH). (36) 


Since the motion under consideration is steady, Eq. (27) becomes 
curl E = 0. (37) 


From Eqs. (24) and (34) to (37) one obtains, by the aid of Eq. (5), 


a ae. 2 
ies HV» (38) 
in which 
n = (dps) (39) 


is the magnetic viscosity. If the magnetic viscosity is negligible, Eq. (38) becomes 


<= 0 (40) 
so that the motion does not change with x. Thus the magnetic field seems to have en- 
dowed the fluid with a certain rigidity, as if the material along the lines of forces had 
been stiffened to give the fluid a fibrous structure. If the motion is axisymmetric (e.g., 
if a sphere is towed through the fluid), from the equation of continuity one can easily 
show that the radial velocity is zero. The motion is then one dimensional in any meridianal 
plane or, if no revolving motion exists, even truly one-dimensional. It may be mentioned 
that for small values of the magnetic viscosity the lines of force always move with the 
fluid, or, to use Alfven’s expression, are ‘frozen’? into the material (Cowling [6], pp. 
5-6). In the case of steady weak motions, the word “frozen” used by Alfven acquires a 
much stronger meaning, as explained above. 

5. Effect of an electric current. Consider now a uniform electric current jo flowing 
in the z-direction, with (r, 6, z) as the cylindrical coordinates. From Eq. (24) one obtains 


H, = Axjor (41) 


the other two components of the main magnetic field being zero. Eliminating j and E 
from Eqs. (24), (27), and (28), one obtains © 
0H P 
ie curl (v X H) + 1VH (42) 
c 
which, for steady motion and negligible magnetic viscosity, becomes 


curl (v X HD) = 0. (43) 
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To the first order, one has, with (u,v, w) for (v,, v9, v.), 
v XH = (—Awu, 0, Hu) 
and 
lo 0 rs] 1 @ { 
curl (v X H) = 4- = (Hw), — - (Hu) + — (How) |, — = (Hew) ?. (44) 
r 00 . or em Oz . ) r 00 ) 


From Eqs. (41), (43) and (44) it follows that 


Ou i 
—— = U, (45) 
06 
O(ru) , dalrw) 
Se ae ee ae (46) 
or Oz 
Ow on 
oP 0 17) 
00 
But the equation of continuity is 
(ru) Ou O(rw) 
: 7 + mt 0). 
or og Oz 
Thus Eq 16) yi Ids 
Ov 
0 18) 


which, together with Eqs. (45) and (47), states that the motion is independent of 8. 
The meaning of this statement can be fully realized physically by imagining an object 
to be towed slowly and steadily along a circular path about an axis parallel to the main 
electric current. According to the result just obtained, the fluid will be pushed along 
circular paths or along the lines of force of the main magnetic field induced by the main 
electric current. What has been said in the previous section about the “stiffening” of 
material along lines of force can be said in the present case also. 

6. Concluding remarks. From the foregoing it may be concluded that, under the 
assumptions stated, the effects of rotation, gravity, or an electromagnetic field are to 
stiffen the fluid along the vorticity lines, isopycnic surfaces or lines, or lines of force, 
respectively. In the case of weak steady motions or relative motions, this ‘stiffening’ 
has the effect of reducing the a priori number of physical or at least the mathematical 
dimensions of motion by one, in the various senses stated. In steady weak motions of a 
stratified fluid the velocity in the direction of gravity vanishes, and simple cases exist 
in which motion normal to vorticity lines or lines of forces are completely inhibited. For 
the motion of a conductive fluid due to the disturbance of an object slowly and steadily 
towed through it along a circular path, a torus-shaped mass of fluid will revolve around 
the axis of symmetry if a main uniform electric current is flowing longitudinally. 

As already mentioned, the concept of a virtual gravitational field created by rotation 
helps to unify the understanding of many phenomena of fluid motion involving gravity 


and rotation. 
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A GENERALIZATION OF LATTA’S METHOD FOR THE SOLUTION OF 
INTEGRAL EQUATIONS* 


By MARVIN SHINBROT (Lockheed Aircraft Corporation) 


1. Introduction. The integral equation 


b 


o(z) | ke —Of)dt, a<x<b (1.1) 


with a and > finite arises in many applications. As examples, we mention the problem 
of optimization of filters havirg a finite memory [1] and the problem of diffraction 
through a slit |2, 3]. In the first of these problems, the kernel represents the autocorrela- 
tion of the inputs to the filter; in the second, k(x) = H4'’(| x |), the first Hankel function. 

Its importance notwithstanding, no general solution of (1.1) is known; with but two 
exceptions, only special methods exist for application to special kernels. The first excep- 
tion is the method developed not long ago by Latta [4, 5]. The second—which is now 
subsumed under Latta’s method—applies to the conceptually trivial case when the 
kernel has a rational Fourier transform [1]. 

Latta’s method can be applied in the circumstance that k(x) satisfies a linear differ- 
ential equation (of any order) with linear coefficients. This restriction is still very far 
from being moderate, of course, and since no general solution of (1.1) appears to be 
forthcoming, it appears worthwhile to see whether there are other integral equations 
having the form of (1.1) which can be reduced to one of Latta’s type. 

In a recent paper [6], Pearson considered (1.1) with 

k(x) = p(x) log | x | + q(x), 
where p and qg were polynomials. Although Pearson’s method has no apparent connection 
with Latta’s, the fact that log | 7 , is a Latta kernel leads one to consider kernels of the 
form 
k(x) = p(x)j(x) + g(x), (1.2) 
where j satisfies a differential equation while, as in [6], p and q are polynomials. It is 
to consideration of this case that the present paper is devoted. 
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One remark about (1.2). Although the restriction of p and q to polynomials seems 


rather severe, the situation is not as bad as it looks. The range of the argument (x — 2) 
of k in (1.1) is the finite interval (a — b, b — a). Consequently, just so long as k(x) is 


given by (1.2) with p and q continuous, the kernel can be successfully approximated 
by (1.2) with p and q polynomials. Thus, given, say, a singular kernel k, if the singularity 
alone is of Latta’s type (e.g., log | x | or| 2 |” or, for that matter, H'’(| x |)), the 
equation can be solved approximately by writing k in the form (1.2) where 7 contains 
the singularity and by approximating p and q by polynomials. 

The paper begins with a brief discussion of the reduction of an equation with a 
kernel (1.2) to an equation with a Latta kernel. In the rest of the paper, three examples 
are worked out: the first two have logarithmic singularities; the last has the singularity 
isl. 
2. The idea of the method. Consider now (1.1) with k(x) as in (1.2). It will be 
convenient to follow Latta and use operational notation—thus, we set 


b 
[ k(x — t)f(t) dt = Af 


ad 
| x — t)f(t) dt = cy, 


Ja 


where the kernel of If is the 7 of (1.2). 


Now, the polynomial g of (1.2) is immaterial to the analysis; if g(x) = >> q,2", we 
have 
| gx — Of) dt = p> a | (x — t)"f(t) dt 
e a 
= 2 gt", 
where the constants g* depend on the moments f° t’f(t) dt, v = 1, --- , n of f. Thus, the 


term involving q(x) can be incorporated into the function g(x) of (1.1); after the resulting 

equation has been solved for arbitrary values of the g* , the solution can be substituted 

back into (1.1) to determine them. We may write, therefore, without loss of generality, 
k(x) p(x)j(2x), 

in place of (1.2). 

Set p(x) = >> p,x”. We then have 


ob 


\f > pn, > (- y( )e / | ix — OU f(H dt 


a 


= >'p, >, (—1)"(")2" 'Tx’f. (2.1) 
n y V 


Now, Latta’s condition that j(z) must satisfy a differential equation allows one to relate 
T'zf (and, as a consequence, Iz’f) to I'f by a linear differential equation. Using this fact, 
one can derive from (2.1) a differential equation relating ['f and Af. However, Af is 
known by (1.1) to be equal to ¢(z). Solving this differential equation between If and 


Af will thus result in an equation 


Wx) = Tf (2.2) 
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which must be satisfied if (1.1) is. The new equation (2.2) can then be solved by Latta’s 
method since its kernel, j, satisfies his assumptions. 

This will now be illustrated by examples. 

3. First example. For our first example, we set 

k(x) = x log | z |. (3.1) 

This kernel itself seems to be of Latta’s type, since it satisfies a differential equation 
(xk’’ — 1 = O) of the appropriate kind. However, k(x) is non-singular; consequently, 
we may expect that no solution will exist unless the function ¢ of (1.1) is properly chosen. 
Since this fact violates a second assumption made in [4] and [5]—that there is a unique 
solution—Latta’s method cannot be applied directly. We have chosen (3.1) as our first 
example in spite of the apparently complicating factor of possible non-existence of a 
solution because of the simplicity of the corresponding relation between Af and If. 

To begin, consider (3.1). With the notation of (1.2), we have j(x) = log | x |, so that 


(x — )j(e -)—1=0 
d/dx). As discussed in [5], this implies 
I’xf = cI’f — mm, (3.2) 
where yuo is the zero order moment of f. Now, obviously, 
Af = zIf — Izf. 
Therefore, 
A’f = zI’f + Tf — Iaf 
= If+ mo, 
using (3.2). Thus, if f satisfies (1.1) with k(x) as in (3.1), it must also satisfy (2.2), where 
¥=¢' — m- (3.3) 


IXquation (2.2) can be solved by Latta’s method if g(x) is a polynomial or an expo- 
and [5]) or, for general ¢, 





nential polynomial (¢ is restricted to this class of functions in [4 
: aot ; we 
hy Carleman’s formula ([7]; see also [6]). As an example, we choose g(x) = x,b = —a = 1. 


By (2.2) and (3.3), we then have to solve 


rf=1— pm. (3.4) 
Since 1 — po is a constant, the solution, according to Carleman’s formula, say, is 
Ko — ] 2 2 9 &) 
= —(l1—27). (3.5) 
x log 2 
Integrating, we can find wo from the equation 
1 
Mo — 1 2\-1/2 
= -— (1 — 2°)’ dx 


x log 2 J_; 


Mo — 1 


log 2 
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Indeed, 
l 
“ 1 — log 2’ 
and 
(1 — 2*)” 
f m(1 — log 2) 


That this really is a solution can be verified by substitution into the integral equation. 


4. Second example. Now, let 


k(x) = (1 + az) log | z |. (4.1) 
Equation (1.1) with k(x) as in (4.1) and g(x) = 1 was solved for small a by Pearson 
in [6]. 
Here again we have j(x) = log | z |, so that (3.2) continues to hold. Also, 


Af = (1 + az)I'f — of zf, 
so that 
A’f = (1 + az) I’f + alf — al’ xf 
I’f + alf + apo , 


using (3.2). Thus, in this case, f can only satisfy (1.1) if it also satisfies (2.2), where 


ome }« 


vy’ tay = go’ — ap. (4.2) 
It remains, then, to solve (2.2), where y is some solution of (4.2). 
We now consider Pearson’s example, setting g(x) = 1, b = —a = 1, and solving 


(1.1) for small a. Equation (4.2) gives 
(x) = —po + Ae 
(A — po) — Aar + O(a’), 


Il 


where A is a constant. 
Define f,(z) by the equation 


rf, = l, (4.3) 
so that 
7 a ae 
sin r log 2 
(see (3.4, 5)). Of course, 
al 
I’zf, = -| fo(a) dx 
: l 
~ log 2’ 


by (3.2), so that 
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x 
Pho = Tog 2’ 
since zf, is an odd function. 
We wish to consider 
rif=y 
= (A — wo) — Aar + O(a’) (4.4) 


= (A — po) Ifo — (Aa log 2)T xf, + O(a’). 


Since the equation ['f = y has a unique solution [7] whose singularities can only be of 


the form (1 — 2z°)~'” at worst [5], it follows that 
f = (A — po) fo — (Aq log 2)zf, + Ofa?(1 — 2”)7'”?] 
‘a —A Aa (4.5) 
= (1 — 2’) | = Sigs ora’) | 
m log 2 T 


The constant A can now be evaluated by integrating (4.5). We obtain 


- or a Ho — A (2 
te [ fiz = "> + OC), 


so that 
A = pl — log 2) + O(a’) 
and 
f=(1-— ay rl + me — log 2)x + or) | (4.6) 


To evaluate yo , recall that we must have 
l1=A o 
j (4.7) 
= (1 + ax)If — af sf. 
Now, (4.6) is equivalent to 


f = —(uo log 2)fo — (wom log 2)(1 — log 2)zfo + Ofa*(1 -- x*)~'”*). 
Hence, 
If = —(uo log 2)I'fo — (uoa log 2)(1 — log 2)l'xf. + O(a’) 
= —py log 2 — woa(l — log 2)r + O(a’) 
by (4.3) and (4.4). Also, 
—(uo log 2)T'xfo + O(a) 
= —pr + O(a). 


II 


I'af 


Therefore, (4.7) results in 
1 = (1 + ax)[—po log 2 — poa(l — log 2)z] + woar + O(a’) 


—po log 2 + O(a’), 
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so that 


Mo = - 5 + Oe’), 


log 2 


giving finally, 


2-179 | _ a (_t ) fa? 
f i 2) | mr log 2 r \log 2 te + 01a" |, 


which agrees with Pearson’s result. 
5. Final example. For our final example, we choose 
k(x) = (1 + az’) | 2 |" ' +1, 0<rv<1; 


this will illustrate a kernel with a non-logarithmic j and a non-zero q [see (1.2)]. The 
integral equation is 


ab b 
g(x) = | [1 + a(z — #)*)|x2—t|"" fijdt + [ f( dt. (5.1) 
The last term in this equation is a constant, up , and so if we set 
% =9— bo, (5.2) 
we have 
ab 


¢,(t) = | k(x — of(d dt 


= Af, 
where 
k(x) = (1 + az’) [2 |". 
Now, if j(z) = |x|", 
ai’ +(l—»j =0, 
and so if If = {2 j(x — é)f(t) dt, 


I’cf = cI’f + (1 — Myf. (5.3) 
Also, 
A,f = (1 + aa’)I'f — 2axT2f + aS x’f. (5.4) 
We obtain from (5.3) that 
xf = 2°0''f + 2(2 — v)xl’f + (1 — »\(2 — Ty; (5.5) 


therefore, differentiating (5.4) and using (5.3) and (5.5), we obtain 
AYf = Tf + ar(1 + v)Pf. 
Thus, in order to solve (5.1), we must first consider the equation 


rf = y, (5.6) 
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where y is a solution of 
Vv’ + axl +r)y~y = ¢', , 


Iquation (5.6) can be solved by Latta’s method if y (and, therefore, g) is an exponential 
polynomial. The details of the computations, though complicated, are not impossible, 


but it would serve no purpose to complete them here. 
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ON THE FIRST STABILITY INTERVAL OF THE HILL EQUATION* 
By C. R. PUTNAM (Purdue University) 


Let \ denote a real parameter and let f = f(t) be a real-valued, continuous periodic 
function of period 1. It is known (Liapounoff) that the Hill equation 

zr’ + (A + f(d))x = 0 [’ = d/dt, f(t + 1) = fd] (1) 

is stable for \ = 0, so that every solution of the equation x’ + f(t)x = 0 is bounded, 


whenever 


f>0, f#0 and [sass (2) 


see, e.g., [1], [6]. Moreover the constant 4 of (2) is the best possible in the sense that 
(2) is not sufficient for the stability of «’’ + fx = 0 if the 4 is replaced by 4 + « (e,a 
positive constant) [3]. If A, and A, denote respectively the left and right end-points of 
the first stability interval of (1) then the first two conditions of (2) imply A, < 0 while 
all conditions together imply \, > O (and so \ = 0 is interior to the first interval of 
stability of (1)). Actually the inequality 4, > 0 is implied by the single condition 


/ f' dt <4, where f'(t) = max (0, f(d] (3) 


(see [6]); moreover, the estimate 


4, >4- [ {f‘dt= a(1 _ if f* at) (4) 
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easily follows. The object of this note is to obtain a “‘best possible” refinement of (4), 


namely, tf (3) holds, then 
1 al 
i> (1 —¢ | f’ at), (5) 


where the equality holds only if f = 0. 

It is known that when A = X, , Eq. (1) possesses a “half-periodic”’ solution x = y(t) 
satisfying y(t + 1) = —y/(t); see, e.g., [5]. In case f = 0, this function is sin zt which 
satisfies 2’ + x°x = 0 (i.e., (1) for f = 0 and \ = X, = 2’); thus the equality of (5) 
holds in this case. It is clear then that the italicized assertion above can become false 
if the x of (5) is replaced by x + e, even if f # 0. That the assertion can become false 
if the + of (5) is replaced by } — e€ is obvious from an earlier remark concerning the 
stability criterion furnished by (2). 

In order to prove the italicized assertion, let y(¢) denote the half-periodic solution 
of (1) for \ = X, considered above, so that y”’ + [A, + f(OJy = 0, and let a, a + 1 denote 
two zeros of y. A multiplication by y of both sides of the last equation, followed by an 
integration and an application of the inequality f < f*, leads to 


nati nar a+ 


1 
A, | y dt > | y” dt — | fy’ dt. (6) 


Ja 


Since 2y(t) = fi y’ ds — f7*' y’ ds, it follows that 2 | y(t) | < f2*' | y’ | ds and hence, 


by the Schwarz inequality, that 4y*(t) < [2*' y” dt. Consequently 


[~ fy dt < ; | 


nati na+l 


f* dt | y”” dt; 


a 


hence, by (6), 
aati al sa+l 
. ] . ” 
A; | y dt> (1 =3 | f at) | y”” dt. (7) 


In view of (3), which implies \, > 0, and the Wirtinger inequality 2° [2*' y’ dt < f2*' y” dt, 
where the equality holds only if y = const. sin x(t — a) (see, e.g., [2], p. 184), relation (7) 
implies (5) and the proof of the italicized assertion is now complete. 

Remark. The referee has pointed out an interesting parallelism between the in- 
equality (5) and an inequality in a recent paper by J. Peetre [4], see pp. 16-17, for the 
least eigenvalue of an eigenvalue problem for partial differential equations with a Rie- 


mann metric. 
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A STABILITY CRITERION FOR QUASI-HARMONIC VIBRATIONS* 
By AUREL WINTNER (The John Hopkins University) 


1. Let w = w(t) be a continuous function which is periodic, choose its period to. be z, 
and call 


xr’ +w (iz = 0 (1) 


a stable oscillator if all solutions x(¢) (i.e., if two, linearly independent solutions) stay 
bounded as | ¢ | — @. It is understood that w(t) is assumed to be real-valued, which 
means that w(t) = 0. Here the sign of equality is allowed (for certain t-values, which 
ean be isolated or form intervals), but it will be assumed that w(t) does not vanish 
identically, i.e., that the case x = 0 of (1) is excluded. 

According to Liapounoff [1], (1) must be stable if 


a® 


i w(t) dt < 4/x. (2) 


But it is also known [2] that if 4/7 is replaced by any numerical constant greater than 
t/x, then the assumption which thus results from (2) will allow the instability of (1) 
(for certain frequency functions w(t), to be chosen in an appropriate manner). In other 
words, 4/z is the best absolute constant in (2). 

This is not, of course, at variance with the fact that if w(t) is the constant 1, then 
(1) is stable but (2) is violated (since x > 4/7). But it can be expected that the absolute 
constant 4/z will be improved if it is assumed that w(t) keeps away from 0 by a fixed 
amount. In this regard, a result of P. 8S. Gusarova [3] states that if 


w(t) = 1, (3) 


then 


. 


[ w(t) dt S 12/x, (4) 
“0 
a relaxation of (2), will guarantee the stability of (1). 

2. This appears to clear up the situation in a satisfactory manner. Actually, the 
stability criterion (3)—(4) is unsatisfactory from the point of view of most applications. 
lor, from such a point of view, the stability of (1) ought to be ensured by imposing a 
limitation on a characteristic which measures the relative degree of the non-constancy 
of w(t), rather than on a characteristic measuring the amount by which w(t) keeps away 
from 0. But the criterion (3)-(4) is of this latter kind. 

Needless to say, the criterion of the desired type will also have to involve a limitation 
of the mean-value of w(t) over a period; see (2) or (4). In the desideratum, the latter 
limitation, supplied by the bounds 4/7, 12/z in the respective criteria (2), (3)—(4) of 
Sec. 1, will have to be expressed in terms of the characteristic measuring the non-con- 
stancy of w(t). 

3. The simplest choice of such a characteristic is the quotient 


0 = (M — m)/(M + m), (5) 
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where 
m= min w (2), M = max w‘(0), (6) 
Ostss Ostss 
hence 
Osm=M<eo and M>O0O (7) 


(in fact, if 17 = 0, then (1) reduces to the unstable, but trivial, equation x’ = 0, excluded 
7) that the quotient (5) exists 


since the beginning of Sec. 1). It is clear from (6) and 


and is subject to the inequalities 
7 ae Me 2? (8) 


and that (1) is the closer to the state of a harmonic oscillator (8 = 0, w(t) = const. > 0) 


the closer is 3 to 0. For this reason, the following criterion, the deduction of which is 
the purpose of this note, is along the lines of the desideratum formulated in See. 2: 

If the coefficient function of an oscillator (1) does not vanish identically and is of period x, 
then all solutions x(t) of (1) stay bounded, as | t | — ©, whenever the mean-value of the 
squared “instantaneous frequency” w(t) satisfies the inequality 


oe oe Se 
| w(t) dt < VU — wv)” (9) 
T 


70 
where 3 denotes the characteristic of w(t) defined by (5) and (6), and Q is a certain numerical 
constant, 
Q = 0.8043 --- (10) 


It will remain undecided whether (10) is the best absolute constant allowed in this 
stability criterion. But it is seen from (10) that the value of the product 7Q certainly 
exceeds the value of 4/z. Hence, if # is reasonably small, then the criterion (9) is more 
favorable than Liapounoff’s criterion (2). 

Note that, the quotient (5) being a relative (“percentile”) characteristic of w(t), 
nothing like an absolute limitation (3) is involved in (9). 

4. The sufficiency of (9) for the stability of (1) will be concluded from the following 
criterion, proved by Borg [4] as a counterpart of Liapounoff’s criterion (quoted at the 
beginning of Sec. 1): If f(¢) is a real-valued, continuous function of period 7’, and if 
f(t) is non-negative throughout but does not vanish identically, then all solutions x(t) of 


2’ + f(t)x = 0 (11) 
will stay bounded (as | ¢ | — o@) whenever the coefficient of (11) is restricted by the 
condition 

aT 
| fd) dt < 64P'/37", (12) 
where 

uae 

P = [ (1 + sin’ ¢)"'”’ de. (13) 
In order to identify (11) with (1), let {() = w(t) and T = x. Then (12) can be written 


in the form 
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1 r 
= [ w(t) dt < 64P*/3r". (14) 
Tv “0 


Consequently, it is sufficient to show that (9) is implied by (14). 
5. First, if f,(¢) and f,(t), where u < ¢ S », are arbitrary continuous functions satisfy- 
ing the inequalities 


0O<m, Sf() S$ M,;, where i = 1,2, 


then, by “the (partial) converse of Schwarz’s inequality” [5], the quotient 


t [ filt) auf fel) u/| A(Of.(0) a) 


cannot exceed the square of 
(M,M,/m,m,)' + (m,m,/M,M,)'”. 
If this is applied to 
f,(t) w(t), f.(j = 1 and u=0,0=-g, 


it is seen from (6) that 


ior | w(t) dt < [(m/M)'? + (M/m)' “| f w(t) a| 


Consequently, (14) is certainly satisfied if 


1 a* : > ad 2 
| w(t) dt < 16F” /= 


r J, 3°7[(m/M)"? + (M/m)"?) te 


In view of the end of Sec. 4, this means that the proof will be complete if it is shown 
that (9) is implied by (15). 
6. It is readily verified from (5) that 
(m/M)'? + (M/m)'? = 2/(1 — 8’)'”. 


Hence (15) can be written in the form (9) if Q is defined by 


Q = 2-3°'*P”/z’. (16) 
Finally, the substitution a = sin ¢ transforms (13) into 
»l 
r | (1 — a‘)? da. (17) 


But (17) is the celebrated constant of Stirling [6], the complete lemmiscatic integral 
(of first kind). Since 


P = 1.3110:--, (18) 


(10) follows from (16). 
7. As an illustration, consider Mathieu’s equation 


x’’ + (a + b cos 2t)r = O, (19) 
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where a, b are constants. This is an equation of the form (1), with a positive w(t) of 


period z and of mean-value a, if 0 < b < a. Since (6) reduces tom = a — band M = 
a + b in this case, it is seen from (5) that 


1 — 8 = 4c/(1 +c)’, wherec = b/a. 


Hence, condition (9) becomes a < Qe'”*/(1 + c). 
Since c = b/a, the latter inequality can be written in the form 


a/b < 4Q’/(a + b)*, where 0 < b <a. (20) 


Consequently, Mathieu’s equation (19) is stable (and, what is more, ‘“disconjugate’’) 
whenever its parameters a, b satisfy the conditions (20), where Q is the value (10), 


defined by (16) and (17). 
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A HEAT CONDUCTION PROBLEM 
INVOLVING A SPECIFIED MOVING BOUNDARY* 


By R. E. GIBSON (Jmperial College, London) 


1. Introduction. Physical problems governed by an equation of heat conduction 
type occasionally arise which require a condition to be satisfied on a boundary moving 
at a specified rate. In the problems with which we are concerned the movement of the 
boundary is generated by the accretion there of new material and not by deformation 
of the medium. Typical examples are: the dissipation of excess pore water pressure in 
a clay layer upon which fresh material is being deposited [1], and the type of thermal 
problem considered by Benfield [2] relating to the theory of the formation of the earth 
on the dust cloud hypothesis. The moving boundary is also a feature of thermal problems 
where a change in state of the conducting medium occurs [3]. Here, the motion of the 
boundary separating the phases is not given and must be determined as part of the 
solution. We can hope, therefore, that our problems, where the motion of the boundary 
is given, will prove to be less intransigent. 

While working on one-dimensional problems of this type, we found that exact solu- 
tions could be obtained in two cases, namely, when the rate of movement of the boundary 
is proportional to t'’* and when it is constant. Although in few actual problems may 
these conditions be realized, we believe our results to be of interest since it appears 
unlikely that closed solutions can be obtained for an arbitrary motion of the boundary. 

2. Problem with spherical symmetry. We consider first a case with spherical 


*Received October 22, 1957 
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symmetry. Fort < 0in0 < r < o the temperature is constant, v(r, t) = v.. Att = 0 
the formation of a sphere begins, within which an homogeneous rate of heat production 
A(t) per unit volume occurs. Additional material, at a temperature v, , arriving on the 
surface of the sphere also participates in the generation of heat at a rate A(t), the current 
radius of the sphere being R(t). 

The appropriate equation for the temperature is 


1 Ov 10f.ov A(t) ; 
= <= “ - "r.C«< 3 < ’ 
k ot r (; 2) + K OsrskO (1) 


where the constants x and K are respectively the diffusivity and thermal conductivity 
of the medium. The problem is to determine v(r, t) for t > 0 subject to the conditions 


v[R(t), t] = v0 , (2) 
Ov 
ay (0, 4] = 0. (3) 


3. Separation of variables. If in (1) we take as new variables u = r(v — vo), § = 
r/R(t) we obtain 


ou chon op Oe , OU 
kK ae t K oF = R ry ERR a’ (4) 
where the prime denotes differentiation with regard to t. The boundary of the sphere is 
now “fixed” at — = 1. We now investigate the conditions for which 
u = f(OG(E) + Eh(d) (5) 
is a solution of (4). This requires, firstly, that 
r A(t) = R“(h’/R — hR’). (6) 
The variables f and G are found to be separable only if R = ct'”*, and then (7) 
GP. i ee 
G E dé” t dg pi Ss eal 8) 


where \ is an arbitrary constant. From the second of (8) 
f= Det, (9) 


where D is an arbitrary constant. Now, in terms of the new variables the boundary 


conditions (2) and (3) are 


f(0GQ) + h() = 0, (10) 
f()G(0) = 0. (11) 
From (9) and (11) we have 
h(t) = —DG(iyt 


and therefore, from (5) 


u = Dt{G() — &G(\)). (12) 
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The constant D is found from (6) to be given by 


K A(t) = —DG(1)t*7(A — 3). (13) 


The rate of heat production is therefore restricted to the form 
A(t) = Aot’, s>-l. (14) 
From (12)—(14) we then have 
- the of fi F oe | 
" © + 1) [* ~ GQ) 


The function G is now found from (8), one of the arbitrary constants arising being 
determined from (11). Finally, we arrive at the following expression for the temperature, 


: 5 3s é | 

KA ets Pils r ; es J l 2 “) 

K (s+ 1) 5 8 e¢ ~ 4\0et 
Fiis+ 5.54 ] 


which expression may be reduced to a finite combination of elementary functions if 
either s or (s + 3) is an integer [4]. In particular, for a constant rate of heat generation 
(g = 0) 


4. Integral equation. In the following section we derive a solution to (1), valid 
for an arbitrary rate of heat generation, when the motion of the boundary is 
R(t) = nt. 
We take 
KA(t) = AF’(b) 
and consider a solution of (1) of the type 
v=v,+ Fit 


pen ae 2 17) 
aie oe | p g(p) [exp - (7 p txt — exp ir -- p) 4xt] dp. : 


Apart from a constant factor the third term represents the temperature at (7, ¢) due 
to an initial temperature g(r) in the region r > O when heat generation is zero every- 
where. Since the sphere vanishes at ¢ = 0 no initial condition on the temperature need 
be satisfied, and we are therefore at liberty to regard g as an arbitrary function which 
must be chosen so that 

v[R(t), t)] = v% , t>0 


and this requirement is met with if g satisfies 


R? = a(p ph p- 
IR(t)F(t)t'” exp = | sin] ‘Xp - dp. 18) 
an P Akt 7 ‘ag tae" 
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In the special case* when R = ut that part of the kernel involving the hyperbolic function 
may be absorbed in the unknown function. The following changes of variable 


p= o, t= 1/4ep 
then reduce (18) to the form 
[ (He dg = 62” iad P(ZE) ex — dS: (19) 
where 
git” ue”? 
H(¢) 2¢ ’ sinh & —_ 


Proceeding formally from (19) we find that 
2 ue 
q(p) PP cosech =f ” Seete r() exp (7 ‘a as p ‘p) d dp 
8x°’" :. Omi - 4kp/ — 16x" 
which, together with (17), leads to an expression for the nh gmt satisfying the 
governing equation (1) and the boundary conditions (2) and (3). In particular, when 


heat generation takes place at a constant rate 


F = Aj,t/K 
the temperature is given by 
Ao aie we ‘A 
1 % +> [ ( ) exp — 
k ur \nrt 4xt (20) 
re sinh Of 4@? coth #2 — 1b exp — & 
| p sinh 7 ius coth & = 1 exp — 4 ao | 


5. One dimensional problems. Problems one-dimensional in some other coordinate 
system may be treated along lines similar to those discussed in the preceding sections. 


For example, the temperature within an infinite cylinder of radius R = ct’’, when 
heat is generated according to (14), is found to be 


ee: [1 _ Fie + 2, 1, r°/4et) | fe ’ “))] 
oT Ke +) Tertid Tau) PP aNnt 


while for the conditions considered in Sec. 4 the temperature is given by 


wherein 


F(t) a | A(r) dr. 


*The solution of Sec. 3 can also be attained directly from (18) after suitable changes of variable 


and the use of the convolution theorem. 
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6. Numerical evaluation. In the spherical problem the quantities of physical 
importance are the temperature at the centre and the gradient at the boundary. When 
heat is generated at a constant rate and the radius increases according to (7), these 
quantities are found directly from (16). We note, firstly, that the temperature difference 
[v(0, t) — v] is a constant fraction 1/(1 + 6xe~*) of the same difference if no leakage of 
heat across the boundary is allowed, and, secondly, as time goes on the temperature 
gradient at the boundary increases without limit. 

On the other hand when R = ut this gradient approaches an asymptotic value 
—xA,/uK which may be seen by considering the behaviour of 


uA Ov 2 8 T/4 ; 3 —y? ry /2 
75 e€ ye” cosech yT'”-dy 


R,t]=7,-1- 
KA, or | 

T(= pt/x) > @. 
The infinite integral has been evaluated numerically along the lines suggested by Goodwin 


[5], and the relation between the gradient and the time factor T' is shown in Fig. 1. 


“pKOv([Rt] 
k A.dr 








r 8 9 10 pet /k 


Fic. 1, Variation of boundary temperature gradient with time. The dashed line represents 


the asymptotic value. 
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ON A GENERALIZATION OF A RESULT OF WINTNER* 
By RICHARD BELLMAN (The Rand Corporation) 


In a recent note [1] Wintner proved the following interesting result. 
Theorem 1. Consider the two equations 
(a) u’’ + f(ju = 0, (b) wv’ + g(t» = 0. (1) 


If there exist two linearly independent solutions of (la), u,; and uz, , such that 
[ (iu PtiuP)lf-¢ldt< « (2) 


then every solution of (1b) can be written in the form 
v = cu, + cu, + of | u, | + |u|). (3) 


This is an extension of known stability results, see [2], to which it reduces if we assume 
that all solutions of (la) are bounded ast — o~. 

Let us now show that we can obtain a generalization of this result following the 
method used in our book [2], to establish the Hukuwara stability theorem, of which 
this will be an extension. 


Theorem 2. Consider the vector-matrix systems 


dx _ ie dy _ 
(a) — A(t)x, (b) tn By. (4) 
Let X(t) be the solution of 
aX _ AX, XQ =! 5) 
_—* ° nella ‘ 
If 
[|| Be) — A} |] || XO | || XW |] dt < @ (6) 
then every solution of (4b) may be written 
y = Xc + o( || X ||) (7) 


ast— o~, 


The norms of matrices and vectors are taken to be respectively }>;,, | x,; | and 
Ds | a | 
Proof. Write 


dy 


= A(t)y + [B() — A(dly. (8) 
dl 


Then, if (0) b, we have 


y = X(b 4+ [ X()X°>(s)[B(s) — A(@®]y ds. (9) 


*Received February 6, 1958. 
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432 
Hence 
{| y < || X(d) bl] + | | X(?) X'(s) |! || Bis) — A(s) | || y {| ds. (10) 
Thus, if we set 
u(t) = || X- ‘(0 B(t) — A(é) |! || y( II, an) 
v(t) = || B(t) — A(d) X(t) X'(é II, 
we obtain the scalar inequality 
uscuvtov | u ds. (12) 


This yields, as a consequence of the fundamental inequality, [2], or directly, the estimate 


| uds < ¢, | v(s) exp | v ar| ds. (13) 


By assumption {“ v ds < o. Hence the integral 


a x 


| X~"(s)[B(s) — A(s)]y ds (14) 


“0 


converges. This means that we can write (9) in the form 
y= X()b+ X(t) | X7\(9[B@) — Aly ds — X() | X-“@[BG) — Ally ds (15) 


which yields the stated result. 
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BOUNDS ON THE ERROR IN THE 
UNIT STEP RESPONSE OF A NETWORK* 


By PAUL CHIRLIAN (New York University) 


Summary. Bounds have been placed upon the error, in the unit step response of a 
network, caused by a departure of the transfer function from its desired value. The 
theorems developed place bounds on the error when it is caused by: 

1) A departure of the amplitude function from zero, for frequencies above a stated 

cut-off frequency 

2) Any deviation of the amplitude or phase function from its desired value. 
These bounds are more readily evaluated than the actual errors, and thus prove useful 
in the design of networks. 

*Received November 29, 1957. This paper is based on a portion of a thesis which has been accepted 
by the faculty of the Graduate Division, College of Engineering, New York University, in partial 
fulfillment of the requirements for the degree of Doctor of Engineering Science. 
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Introduction. In many analysis or synthesis procedures, the transfer function of 
« network only approximates the desired value. It is often necessary to know how this 
approximation affects the transient response. Existing methods require tedious analytic 
or graphic analysis to determine the deviation of the transient response. More readily 
evaluated bounds, which eliminate the need for this analysis would be useful. Such 
bounds have been developed in the theorems presented in the body of this paper. 

The first considers the error, in the unit step response, caused by an assumption 
that the transfer function is zero above a given cut-off frequency. 

The remaining theorems deal with the errors, in the unit step response, which result 
from deviation of either the amplitude or phase functions from their desired values. 

Throughout this paper, the transfer function of the network in question will be 
written as T'(w)é where the amplitude function 7T(@) is never negative and T(0) > 0. 
In addition the magnitude of the per cent error in the unit step response, referred to its 
final value will be written as | %6 |. That is, if A(¢) is the actual unit step response and 
1,(/) is the approximate unit step response, then 

.; . | A(t) — A(t) | max X 100 
iad 7(0) 

Bounds on the error in the unit step response caused by a band-limiting approximation. 
It is often convenient to assume that the transfer function of a network is zero for all 
values of frequency above a given cut-off frequency w, . This assumption, calied a band 
limiting approximation, introduces an error into the calculated value of the unit step 
response. It is desirable to bound this error so that the approximation can be justified. 
This bound can be obtained by evaluating an integral of the form 


r” Tw) 1 
|] T(0) w 
This evaluation may be simplified by replacing 7(#)/7T(0) by a new, more readily 
integrated function G(w), where G(w) > T(w)/7T(O) for all w > w. and J2, [(G(w)/wdw 
exists. The best bound on the error will be obtained when G(w) = 7'(w)/7T(0) and this 
procedure should always be considered as a first choice. 
Theorem | will bound the error caused by the band limiting approximation. 


dw. 


Theorem 1. If T(w)/T(0) < G(w) for w > w, , then 


100% f° Glw) 
| — dw. 


T Jae W 


% § < 


c 
Proof. For any network the unit step response is given by 


ap ray 


A(t) = 3T(0) + | 
WT Jo 


(w) 


sin [wf + O(w)] dw. 


Thus in this case the error is given by 


100% | ¢” Tw) 


|; = si ») ) 
> 6 . [ oT(0) sin [wt + Aw)] dw 


The bound is then obtained by utilizing the fact that T(w)/7T(0) < G(w) and | sin [wt + 
#(w)| | < 1. Theorem 1 will now be applied to a particular function which falls off as 


|/w", where n > 0. This is of use in many commonly encountered networks. If G(w) = 


e(w./w)" the following is obtained. 
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Corollary la. If T Ww T(0) _* El W w)” for Ww > Ww and if] =~ 0 then 


an ee 100€% 
2 , 
nw 
Some typical values given by this bound are 
If « 0.1 and n 1 then | %6| < 3.18% 
If ¢ = 0.707 (8 db down) and n = 3, then | %6| < 7.5%. 


Bounds on the error in the unit step response caused by an approximation of the 
amplitude or phase function. When an arbitrary transfer function is not realized 
exactly, but is only approximated, an error appears in the unit step response. 

The following theorems place bounds upon this error. The first of these deals with 
an approximation of the amplitude function. In this case the transfer function of the 


desired network is 7'(w)e’’ “>. Often this transfer function is approximated by the transfer 
function [7(w) + 7'.(w)] exp [j@(w)]. The amplitude error function 7',(w) is a real function 
of w which may be positive and negative. Then if {% | 7.(w) |/w dw exists, | %é | can 


be bounded. 


Theorem 2°. If the amplitude error function T .(w) exists, then 


100% [ T .(w) | 1 
T Jo T(0) | a 


dw. 


Proof. This theorem is proven by formally writing the expression for the error and 
then utilizing the fact that | sin [wt + 6(w) 5 

If a function M/(w) can be found such that | 7’.(w)/7(0) | < M(w) for all w and the 
J% [A1(@)/w| dw is more readily evaluated than the integral of Theorem 2, then 


Ww 


lA 


term 
it may prove convenient to replace | 7.(w)/7T(0) | by A/(w) in this theorem; the bound 
developed in this case will be weaker than the one originally presented. However, in 
some cases, it may be expedient to use the simpler integral to obtain a rough estimate 
of the error. 

The next two theorems will place a bound upon the error in the unit step response 
produced by an approximation of the phase function. The transfer function of the desired 
network is 7T(w) exp [j@(w)]. Many times this transfer function is approximated by 
T(w) exp {j[0(w) + o(w)|}. The phase error function ¢(w) is a real function of w which 


may be positive and negative. If the integrals 


x ry 1 ~@ say \ 2 
Tw) | b(w) T(w) [(w) | 
= = da and | , dw 
7(0) w , Jo TO) ow 


exist then | %6 | may be bounded. 
Theorem 3. If the phase error function $(w) exists, then 
100% f® T(w) 1 . 2 | bw) | 
07 nN od tf | {9 2 ' / 
6 — 2 sin sin d(w) | ? dw 
, tr Jo T(O)w | 2 id ~~ ) i 
4 n© IY fae \ 72 ) 
100% | T'(w) 1 } [o(w) | mi { 


: p(w) dw. 
T Jo IO) w | 2 oe f - 


'The proof of this theorem is similar to one given by A. H. Zemanian, An approximate method of 


evaluating integral transforms, J. Appl. Phys. 25, 262-266 (Feb. 1954) 
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Proof. The first bound is obtained by formally evaluating the error and then utilizing 
the fact that | sin [wf + @(w)] | < 1 and | cos [wt + @(w)] | < 1. The second equation 
is obtained from the first by realizing that | sin x | < | z |. The first of these equations 
will provide a smaller bound than the second. However, the second equation is more 
readily evaluated and if ¢(w) is small, in the range of integration, then the results obtained 
from this second equation can be quite satisfactory. 

If the normalized amplitude function T(w)/T(Q) is itself bounded, the results of 
Theorem 3 are readily extended to obtain a bound which is more easily evaluated. 


Corollary 3a. If T(w)/T(0) < M then 


00M% r°1s... — | 
%Msel< m Sb f [ . $9 sin? [#2 | + | sin ¢(w) > dw, 
Tv Jo W | 2 J 
and 
7 _100M% 1 )- 
S315 = | : {lo} + | dw) de. 
T Jo Ww 2 


When networks whose amplitude functions are such that Theorems | or la indicate 
that a negligible error is introduced by neglecting the response at frequencies greater 
than w, then the upper limit of integration can be changed from © to w, in the equation 
of Theorems 3 and 3a. This technique should especially be used when the effects of phase 
correction is to be determined since extremely large phase errors usually occur at fre- 
quencies above w, . These phase errors usually have little bearing on the transient 
response. 

Conclusion. The theorems which have been developed here place bounds on the 
error, in the unit step response, caused by a deviation of the transfer function from its 
desired value. 

In transient analysis procedures, it is often convenient to neglect the response of a 
network above a given cut-off frequency. Theorems 1 and la bound the error produced 
by such procedure. In addition, knowledge of results obtained in Theorems 1 and la 
can simplify the design of a network, since these results indicate a frequency above 
which the response is no longer important. Theorems 2, 3, and 3a bound the error caused 
by a deviation of either the amplitude or phase function from its desired value. Deviation 
of this type will occur in synthesis procedures, where the desired transfer function 
cannot be exactly realized with a finite number of realizable elements. 

When a network is designed to produce a desired unit step response, the error of the 
transfer function, in itself, has no direct significance. However, this is often the known 
error. Theorems 2, 3, or 3a can then be employed to relate the error in either the ampli- 
tude or phase function, to the error in the unit step response. 

The use of these theorems, therefore, can provide valuable information to the net- 
work designer. 

Acknowledgments. The author is indebted to Professor James H. Mulligan, Jr. 
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ON AN INITIAL VALUE PROBLEM CONCERNING 
TAYLOR INSTABILITY OF INCOMPRESSIBLE FLUIDS* 


By G. F. CARRIER AND C. T. CHANG (Harvard University) 


Abstract. By taking the effect of surface tension and viscosity into consideration, 
Bellman and Pennington [2] have generalized the original treatment of the problem 
of Taylor instability [1]. They claim, however, that a problem in which the motion 
starts from rest cannot be treated with their linearized formulation. In this paper, the 
initial value problem is treated and, although the algebra becomes more complicated, 
the linearized analysis suffices. In particular, the cut-off wave number as found in [2] 


is not modified. 


1. Introduction. The problem of Taylor instability [1], asks about the instability of 
a small disturbance imposed on the interface separating two infinitely extended fluids 
when the system is undergoing a constant acceleration directed from the lighter towards 
the heavier fluid. 

The original treatment was limited to incompressible nonviscous fluids. Surface 
tension effects at the interface were ignored. Referring to a simple sinusoidal initial 
disturbance of amplitude a and wave number k, the result showed that the time history 
of the surface disturbance with respect to the accelerated frame of reference is given by 


n(x, t) = a cosh (nt) cos kz. (1.1) 
The growth factor ‘‘n” is related monotonically to the wave number ‘“‘k’’ by the 
expression 


2 Po — Pr 
———_—— g®k. (1.2) 
P2 + Pi 
In the above, p; , p2 are the densities of the lighter and the heavier fluid respectively, 
while 
gQ=-gton (1.3) 
with g as the usual gravitational constant, and g, , the imposed acceleration. 
Bellman and Pennington [2], later generalized Taylor’s original result by taking the 
effect of surface tension and viscosity into consideration. They assumed a standing 


ai 


wave type solution and obtained the following expression for the growth factor ‘“‘n” 


saa ae T 
. <a — i, (1.4) 


=. 3 a 
P2 + Pi . P2 + pi 


where T is the surface tension of the interface. Equation (1.4) indicates the existence 
of a cut-off wave number k, , i.e. a wave number such that no instability arises when 


k > k. , where 
g* |'? i 
k, = | -. = ee) o* | : (1.5) 


a 


*Received November 14, 1957. 
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Their result also shows that viscosity tends to retard the growth of the interface. 
However, they note a dilemma, for according to their solution it is impossible to start 
a motion from rest. They attributed this seeming paradox to the linearization performed 
in their analysis. 

The object of this paper is to show that this paradox is easily resolved. As a con- 
sequence, the time history of the interface is given by a more complicated expression 
than that of Eq. (1.1); in particular, it is no longer possible to define a simple growth 
factor. However, this algebraic complication, neither compromises the existence of the 
cut-off wave number nor changes its value as given by Eq. (1.5). 

2. Formulation of the problem. Tor simplicity, we shall confine our attention to 
an air-liquid system, i.e. we may approximate the density ratio by p:/p2 = 0. The more 
general problem can be treated in precisely the same way. The flow field is taken as the 
half plane bounded above by the free surface y = n(x, t). The medium is assumed to be 
incompressible but viscous; surface tension effects are to be included. The whole body 
of the fluid is assumed to undergo a constant vertical acceleration g, directed from air 
towards the liquid. Referring to a stream function y, , the linearized momentum equation 
with respect to the accelerated frame of reference can be written as 


00 » > > 
——{[Vy-r~VV yy] =0 (2.1) 
Oy ot 
and y itself must satisfy 
Q i ° . 
= Vy=rV' Vy. (2.2) 
Cc 


The coordinate axes are chosen in such a way that the z-axis is taken in the unper- 
turbed plane of the free surface while the y-axis is taken normal to it. 

With some amount of algebraic manipulation it can be shown that y must satisfy 
the following boundary conditions (they imply ¢,, = ¢,, = 0) at the free surface: 


T 
v\ l rz zu _ xt - g* ry = amu) 
¥ M w/t v ! ‘ v p v (2.3) 
= a + Vows — Q), at y = 0. 
The usual kinematic requirement at the free surface is 
% = — Vey - (2.4) 
With no loss in generality, we can restrict the analysis to a single Fourier component 
in «; thus we may write 
tke 
(x, y, ) = oly, Ye 
Yr, y oly, Ye’, (2.5) 
n(x, t) = o(te’. 
We may also define 


fly, 8) = [ e ‘oly, t) dt 
J0 (2.6) 


o 


u(s) = / e "‘o(t) dt 
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and obtain in place of Eqs. (2.2) and (2.3) the following equations 


a (1+ oe) Nae — ef - EL * 
— — k : ;— kK = — 5 — K DY, 0 . 
| , J a vk 14 , | v dy ' wy ( 


for —~ <y <0, 


bo 
“I 


and 


> > Y hd o 2 8) 
s{k°’f + fu] + 8°f + (9° es Kt, 
p 


Tos 
= s¢(0, 0) — (9° + K* oO), at y = 0. 
k\ p 


rom these expressions one observes that the problem is determined theoretically 
once the initial surface disturbance and the initial velocity field throughout the domain 
are prescribed. 

3. Illustrative example. To illustrate our discussion, let us take 


/ 


¥(r, y, O) = 0 (3.1) 


and 
n(x, O a cos kz, (3.2) 
i.e. the motion is started from rest with an initial disturbance imposed at the free surface 
in the form of a sinusoidal wave having amplitude a and wave number k. 
Using the foregoing method, it can be verified that 7, is given by 


f 4 \ 
N: —asvk°B > A a exp a. = l)k*t] erf¢ pangs (vk*t)' “I? cos kx (3.3) 
where a; are the roots of the quartic 
Wax a a: fi —— i >| 9 
\7) = 7 FT 21 ir+U+86 (3.4) 


with 


riven by 


y 


the quantity A; is 


and the parameter 8 is given by 
ar ‘Ie (* | ae 
sie v” pk (! ~ AE } |: (3.6) 


Here, k, is the cut-off wave number corresponding to p,/p. = 0, i.e. 
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Equation (3.4) indicates that for each j the Re (a’) < 1 when | k/k, | > 1, but this 
inequality does not hold when | k/k, | < 1. Thus we are in agreement with the principal 
result obtained in [2] but note that the model of that paper has broader validity than 


was claimed. 
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WAVES PRODUCED BY A PULSATING SOURCE TRAVELLING 
BENEATH A FREE SURFACE* 


Quarterly of Applied Mathematics XV, 249-255 (1957) 


By PAUL KAPLAN (Slevens Institule of Technology) 


There appears to be an error in the analysis performed in the subject paper by 
Dr. Tan [1], which has led to erroneous conclusions in regard to the character of the 
resulting wave pattern and the asymptotic conditions at infinity that are used to arrive 
at a unique solution. Other treatments of this same problem have been carried out by 
Secker [2], Kaplan [3], and Wu [4], and they have arrived at different conclusions from 
those of Tan. Their results show two harmonic wave trains on the downstream side 
of the source when r > 1}, corresponding to the roots A, and K, , in agreement with 
the results of Tan. However, when 0 < + < 4, Tan indicated that four harmonic wave 
trains were on the downstream side, with no disturbance found to propagate to infinity 
upstream. This led to his conclusion that the imposition of an asymptotic upstream 
condition of “‘vanishing disturbance at infinity” is sufficient to render his solution 
unique. The results of Becker, Kaplan and Wu were all in agreement, and showed that 
when 0 < 7 < } there were three harmonic wave trains on the downstream side, and 
one harmonic wave train on the upstream side. It can be shown that the wave in the 
upstream side corresponds to the root K, and that the error in Tan’s analysis was not 


recognizing the fact that, in the notation of [1], 


lim JImK, = 0- O<7r<} 
not 0+, as he indicated. 

If this error were corrected, and the resulting analysis carried out in the same manner 
as done in the paper, the results would have been in agreement with that of the other 
authors. The imposition of the requirement of no disturbance upstream would not be 
valid, but the use of the technique of Rayleigh’s artifice in order to secure uniqueness 
would not be affected, in spite of the artificial nature of this factor. A more realistic 


*Received November 18, 1957. 
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approach, which has more physical appeal, is to formulate the problem as an unsteady 
initial value problem with the motion assumed to start from a state of rest, and the steady 
state oscillatory solution is obtained when the time t — and ‘‘transients” have dis- 
appeared. The behavior of the waves at infinity, i.e. the radiation conditions, are deter- 
mined directly from that type of analysis by only imposing boundedness conditions at 
infinity, without any a prior? assumptions about the character of the waves. This tech- 
nique has been utilized by Kaplan and Wu in their treatments of this problem 
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BOOK REVIEWS 


(Continued from p. 396) 


Statistical analysis of stationary time series. By Ulf Grenander and Murray Rosenblatt. 
John Wiley & Sons, Inc., New York, and Almqvist & Wiksell, Stockholm, 1957. 
300 pp. $11.00. 


The basic probability model considered in this monograph consists of a discrete parameter sto- 


chastie process y x, + m,, where x; is a stationary stochastic process and m,; is a linear combination 
of known regression vectors. Often in the book 2, is specialized to a linear process, which can be described 
roughly as an infinite linear combination of indepent identically distributed random variables. Since z; 
is stationary, the convariance sequence 7, = E(x; , 2:4,) depends only on an n. The hermitian matrix 

) is non-negative definite, so by a well-known theorem in Fourier series there exists a nondecreasing 
function F(A) which has ro, rm. , *** , a8 its Fourier-Stieltjes coefficients. This function is called the 
spectral distribution function (sdf) of 2. . The central problem of the book is to study the distributions 


f certain statistical estimators of the sdf and of its derivative, the spectral density (sd). These esti- 
mators are quadratic functions of a finite partial realization of the process y,; . The estimation of the 
regression coefficients and tests of hypotheses concerning them are also considered. The continuous 
parameter case is also discussed in several places. Some of the material in the later chapters is here 
published for the first time. 

The first three chapters of the book contain background material. The first chapter presents the 
basic spectral theory of stationary processes. The second chapter contains «a discussion of the least 
squares theory of prediction, interpolation, and estimation of regression constants, when the spectrum 
is known. The language of linear manifolds and function spaces is used to introduce the fundamental 

leas. In the next chapter some examples of the older approaches to time series are given. The models 
considered are finite parameter ones such as the model in the first paragraph above but with independ- 
ent, identically distributed stochastic elements x; , and also the autoregressive case in which a finite 
moving average of the elements y; is set equal to a succession of independent, identically distributed 
random variables. A feature of this chapter consists of a number of carefully polished proofs of some of 


} 


the leading classical results including Slutzky’s theorem. 
The authors start taking up the new approaches in the fourth chapter, in which statistical estimators 
. continuous sd are introduced and their asymptotie biases and variances are derived. Emphasis is 
placed on quadratic estimators which are weighted averages of the periodogram (‘“‘spectrograph esti- 
mates’’). For linear processes they are shown to be as good in a certain sense for estimating the sd at a 
particular point as any of the more complicated types of quadratic estimators. There is a rigorous dis- 
cussion of the asymptotic bias and variance of a certain estimator of the sd in the time-continuous case. 
This is one of the few major demonstrations presented for that case in the book. 

In the next chapter some sketches of the role of spectral analysis in the study of elecéronic devices, 
turbulence, and water waves are presented. In the sixth chapter, distribution theory is again taken up. 
The discussion there centers upon a linear process and upon an estimator of the sdf which consists 
merely of the indefinite integral of the periodogram. The major result can be described roughly as 
follows: Under certain regularity conditions, the asymptotic distribution function, say, H(a), of the 
maximum distance (properly scaled) between the estimator F'*(\) and the sdf F/(\) for 0 SX S = is 
equal to the probability that a certain normal process 2(\) stays inside (—a,a) for 0 < \ S z. The co- 
variance function of z() unfortunately depends upon the sdf, so it is difficult to use the result to set up 
asymptotic confidence bands except in special cases. Extensions to more general estimators and to the 
case where a regression is present are considered. The chapter concludes with numerical illustrations 
in which the estimation techniques are applied to some artificially generated numerical series. 

Various questions arising in regression analysis are considered in the seventh chapter. There is a 
comparison of the asymptotic efficiency of the least square estimator and the minimum variance esti- 
mator. This is followed by a careful study of least square estimation. The exposition is somewhat intri- 
cate and depends upon a device for handling the regression vectors all at once which the authors call 
the regression spectrum. The last chapter deals with miscellaneous problems. Perhaps the most interest- 
ing consists of the problem of prediction when the spectrum must be estimated. The authors present 


a consciously non-rigorous derivation of the mean square error of prediction in this case. The book con- 
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cludes with some problems for the reader and an appendix which presents some of the classical results 
connected with the names of Hardy and Littlewood on functions analytic in the unit circle. 

The stated prerequisite for reading this book is ‘‘a knowledge of statistics and basic probability 
theory equivalent to that contained in H. Cramer, ‘‘Mathematical Methods of Statisties.”’ ’’ With this 
understanding, the exposition is highly self-contained. However, the student who has reached the level 
of mathematical competence at which he can just comfortably read Cramer’s book, but nothing more 
difficult, will find some rough going here. To feel really at home, he will need to have had a good deal 
of experience with the methods of rigorous mathematical analysis, and in particular he should have a 
reasonably close acquaintance with the modern theory of trigonometric series. The proofs which the 
authors give for previously known results are often elegant but tend to be unmotivated and highly 
condensed. 

The study of time series has always presented many challenging and urgent problems. The authors 
claim, with considerable justification, that the highly specified models used in the past have often been 
unsatisfactory. They embark upon a program designed to investigate whether the spectral methods 
which have been so useful in the engineering sciences and in physics can be developed in such a way as 
to give a better theory for general purposes, The investigation is undoubtedly of great importance, but 
as with every other known approach to time series analysis, there are some discouraging facts to be 
faced. For example, knowledge of the spectrum is an end in itself in many engineering and physical 
problems, and this is not always true in other contexts. The statistical estimation of the spectrum is 
somewhat tricky technically, particularly if not too many data are at hand. (For instance, the relation 
between the sample size and the maximum lag of the products used in the estimator must. be carefully 
acjusted.) Accurate point estimates of the sd are elusive. The spectrum of a stationary process does 


not uniq iely determine the process without further assumptions so even if one could estimate the 


spectrum of a given series exactly, not all questions about its structure would be cleared up The random 


elements of the processes encountered in the social sciences may be generally non-stationary, and the 


adaptation of spectral methods to non-stationary series is still somewhat problematical. The analytical 
derivation of exact sampling distributions in the spectral approach may be very difficult. (The statistical 
results in the present monograph are almost all of an asymptotic nature. ) 

These problems are all openly recognized by the authors, and the y are careful to point out the tenta- 
tive character of parts of their monograph. They state in the preface that one of their main purposes 
is ‘‘to stimulate research in time series analysis which will lead to practically useful and theoretically 
sound methods.”’ The reviewer feels that the publication of this book is an important milepost in the 
development of the mathematical theory of time series. The book is rigorous and accurate, and at the 
same time it is full of interesting results and techniques which will probably insure its continuing value 
as a reference whichever way the research in the field turns 
There are the usual minor misprints which characterize a first edition. Perhaps the only accident 


which might confuse the reader takes place on the bottom line of page 86 and page 87, where the z 


rR 


should all be y's 
J. H. Curtiss 


Vector spaces and matrices. By Robert M. Thrall and Leonard Tornheim. John Wiley 
& Sons, Inc., New York, and Chapman & Hall, Ltd., London, 1957. xiii + 318 pp. 
$6.75. 

This book presents the theory of vector spaces and matrices simultaneously at two levels, one 
concrete and one axiomatic, the concrete approach sometimes preceding, sometimes following, the 
axiomatic one. The authors have succeeded in providing a text that aims to develop in the student 
conceptual grasp rather than mere manipulative skill, as well as to aid him in the applications of the 
subject to various branches of pure and applied mathematics. The exposition, throughout clear and 
concise, is more detailed in the chapters forming roughly the first half of the book (often more so than 
is usual), and more advanced in the later ones. Each chapter contains a good number of selected exer- 
cises which frequently supplement preceding theory. Among the features generally not found in texts 
on this subject is a development of the general theory of simple algebraic extensions of a field (chapter 
8) and an introduction to linear inequalities and game theory (chapter 11). 

Table of contents: 1. Vector spaces. 2. Linear transformations and matrices. 3. Systems of linear 
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equations. 4. Determinants. 5. Equivalence relations and canonical forms. 6. Functions of vectors. 7. 
Orthogonal and unitary equivalence. 8. Structure of polynomial rings. 9. Equivalence of matrices over 
a ring. 10. Similarity of matrices. 11. Linear inequalities. There are appendices on mathematical in- 
duction, relations and mappings, and a bibliography. 

H. A. ANTOSIEWICz 


Water waves—the mathematical theory with applications. By J. J. Stoker. Interscience 
Publishers, New York and London, 1957. xxi + 567 pp. $12.00. 


The theory of water waves as a branch of mathematical hydrodynamics may be said to begin with 
Laplace’s work on tides, and during the nineteenth century it attracted the attention of many of the 
most eminent physicists. Their work is described by Lamb who devotes roarly one-third of his book 
to it. After 1900 the subject became unfashionable, interest in wave motion shifted to electromagnetism 
and atomic physics, while interest in the application of hydrodynamics shifted to aerodynamics; only 
a few mathematicians kept working on the numerous problems left unsolved in this field. Some out- 
standing theoretical work was done during this period which is only now receiving its merited recog- 
nition. New activity began during the recent war, and is still continuing. Much new and interesting 
mathematics was discovered, although much still remains to be done. The progress in theoretical knowl- 
edge has been accompanied by advances in hydraulics. A leading part in these developments has been 
plaved by the Institute of Mathematical Sciences of New York University where the well-known author 
of the present. work is Professor of Mathematics. 

This book is the first connected account of the subject since Lamb’s Hydrodynamics, but so rapid 
has been the expansion of the subject that a book of nearly 600 pages cannot now cover all topics, and 
the author has been faced with the difficult problem of selection. The contents of the book are as follows: 


Part I 1. Basic hydrodynamics. 
2. The two basic approximate theories (small-amplitude and long-wave ) 
Part II Subdivision A. Waves simple harmonic in time. 
3. Simple harmonic oscillations in water of constant depth. 
4. Waves maintained by simple harmonic surface pressure in water of uniform depth. Forced 
oscillations. 
5. Waves on sloping beaches and past obstacles. 
Subdivision B. Motions starting from rest. Transients. 
6. Unsteady motions. 
Subdivision C. Waves on a running stream. Ship waves. 
7. Two-dimensional waves on a running stream in water of uniform depth. 
8. Waves caused by a moving pressure point. Kelvin’s theory of the wave pattern created 
by a moving ship r 
9. The motion of a ship, as a floating rigid body, in a seaway. 
Part III 
10. Long waves in shallow water. 
11. Mathematical hydraulics. 
Part IV 12. Problems in which force surface conditions are satisfied exactly. The breaking of a dam. 


Levi-Civita’s theory. 


As is explained in the introduction, the problems included in the book are chosen mostly from 
among those on which the author or his colleagues at New York University have worked, and this 
gives the book a personal character. Completeness is not aimed at, even in the problems treated here, 
and the literature cited is confessedly incomplete, even as regards recent work. The mathematical and 
physical equipment needed by readers differs in different parts of the book, and some of the most inter- 
esting sections (e.g. on ship wave patterns) may be read with little previous knowledge of hydrody- 
namics. There is great variety in the problems. ‘‘The physical problems range from discussion of wave 
motion over sloping beaches to flood waves in rivers, the motion of ships in a sea-way, free oscillations 
of enclosed bodies of water such as lakes and harbors, and the propagation of frontal discontinuities in 
the atmosphere, to mention just a few. The mathematical tools employed comprise just about the 
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whole of the tools developed in the classical linear mathematical physics concerned with partial dif- 
ferential equations, as well as a good part of what has been learned about the non-linear problems of 
mathematical physics. Thus potential theory and the theory of the linear wave equation, together 
with such tools as conformal mapping and complex variable methods in general, the Laplace and 
Fourier transform techniques, methods employing a Green’s function, integral equations, etc. are used. 
The non-linear problems are of both elliptic and hyperbolic type.”’ But, in spite of the diversity of the 
material, a thread of continuity runs through the book. Friction is neglected except in the chapter on 
mathematical hydraulics 

In the space allowed for this review it is not possible to discuss in detail any section ol the book. 
and short comments on a few topics must suffice. One-half of the book treats small-amplitude waves, 
for which the non-linear free-surface condition is replaced by a linearized approximation. In the di- 
vision dealing with time-periodic waves, the section on sloping beaches does not claim to describe ade- 
quately the behavior of real waves on beaches, the interest lies in the mathematical boundary-value 
problem, and in the behavior of solutions near the corner. The literature on this problem alone is now 
so big that only a few of the new methods could find a place in the book (There is a section on break- 
ing waves in Chapter 10). Among other problems treated in this part of the book is Sommerfeld’s problem 
of diffraction by a half-plane for which a new and very pretty solution is given. Many other methods 
and applications of time-periodic linearized waves receive brief mention. The next subdivision, devoted 
to transients, contains among other things a discussion of boundary conditions at infinity (radiation 
conditions) and of the principle of stationary phase. The next subdivision, on waves in a running stream, 


gives a treatment of unsteady waves, leading to the appropriate radiation condition, and goes on to 


describe ship wave patterns due to pressure points moving along straight lines and circles, with photo- 

graphic illustrations. The last chapter in this subdivision treats ships in a seaway in a systematic manner 

The infinitely thin ship (Michell ship) and its limitations are more clearly explained here than any- 
: | 


where else in the literature, and there is reason to believe that this will be the opening of a new phase of 


development in this neglected interesting and technically important subject. A similar treatment for 


planing ships and other ships of thin section is also given. (For a recent extension of this work see A.S 


Stoker, Comm. Pure Appl. Math, 10, 1957, 399-490.) The third part of the book deals 


Peters and J. J 
in Chapter 10 an interesting variety of problems is treated; in Chapter 1! 


with long-wave theor, 
friction is regarded as dominant, and the theory is applied to flood routing on the Mississippi with 


success, using parameters determined from earlier floods and progressively cerrected. Such calculations 
can be made quickly enough to give useful predictions about the later stages of a flood. Comparison 
with other methods of calculation would have been interesting 

I perio he 


The last part of the book discusses exact solutions, and gives a new proof of the existence o 
progressive waves on deep water, first established by Levi-Civita. 
believes that this book has completely achieved its aim. A picture emerges of the 


The reviewer 
many problems have been solved, but some of the most interesting major 


present state ol the subje ct: 


problems remain unsolved, including the following: 
(1) In the linearized theory of time-periodic waves (Laplace’s equation with the simplest: mixed- 
type boundary condition) the uniqueness theory is still missing 


I 
2) The theory of unsteady non-linear wave motion is still very incomplete, with corresponding 


difficulty in the interpretation of experiments 


i 
(3) The developn ent ol experimental techniques and the detailed comparison of theory and ex- 


periment still leave much to be desired 
I’) URSELL 


Turbolenza di Parete. By Carlo Ferrari. Levrotto & Bella, Torino, 1957. 116 pp 


It is now generally recognized that the inner and outer parts of the turbulent boundary layer have 
quite different turbulent structure, and as Professor Ferrari shows in his lecture notes, acceptance of 
this single premise is almost sufficient for rejuvenating the mixing-length theory. Indeed, the inner 
sublayer is dominated by the wall shear to such an extent that dimensional argument and Taylor series 
yield quite detailed information, and the assumption of complete self-preservation in the outer laver 


(critically examined by Townsend, ‘‘Structure of Turbulent Shear-flow,’’ Cambr. Univ. Press, 1956) 
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then suffices for deducing the logarithmic distributions of mean velocity and skin friction. On the further 
assumption of a turbulent shear diffusivity constant across the outer layer (cf. again Townsend), the 
mean velocity distribution can be completed, and this is carried out in detail for channel flow and the 
boundary layer without, and finally with, pressure gradient. 

Of course, little understanding of turbulent structure is gained, but on the other hand, it is possible 
to make the theory appear almost independent of assumptions derived from experiments. In fact, no 
experiments are mentioned, beyond the Laufer-Townsend energy balance of the boundary layer without 
pressure gradient. An advantage of this approach is that the calculations can be extended in a natural 
way also to heat transfer and compressibility, in the manner pioneered by Reichardt and Ferrari, and 


this extension is presented for the case without pressure gradient. 


R. E. Meyer 


Quantum mechanics. By ¥. Mandl. Academic Press, Inc., New York, and Butterworths 
Scientific Publications, London, 1957. x + 267 pp. $6.50. 
This is a second edition of a justly popular book. A final chapter on the Dirac equation, requiring 
1 rudimentary knowledge of the special theory of relativity, has been added to the first edition, pub- 


lished in 1954 
Apart from this, 
tion, remain a most valuable part of the text. 


there are no essential changes. The exercises, together with hints for their solu- 


RicHarp BELLMAN 


The numerical solution of two-point boundary problems in ordinary differential equations. 
By L. Fox. Oxford at the Clarendon Press, 1957. xi + 371. $9.60. 


Chis book presents the most complete treatment of the two-point boundary problem by numerical 
methods which has so far appeared. The nearest rival is Collatz’s “‘Numerische Behandlung’’ which 
devotes about 100 pages to boundary problems for ordinary differential equations. The points of view 
of Collatz and Fox differ considerably, the former leaning more toward analytical, the latter more 
toward purely numerical methods. 

The scope of Fox’s book is shown by an examination of the contents. Chapter 1 gives a brief in- 
troduction to the general problem, and emphasizes the distinction between two-point and _ initial- 
point problems, a distinction which is far more serious for numerical methods than for purely analytical 
methods. Chapter 2 presents finite differences and difference equations, supplying the basie compu- 
tational background needed in later chapters. In Chapter 3, the author studies the solution of algebraic 
equations, both linear and non-linear, by direct methods, by iteration, and by relaxation. The method 
of relaxation is then applied to simultaneous finite-difference equations. 

Up to this point the material has been preparatory. In Chapter 4 the author takes up second order 
differential equations with two-point boundary conditions. Various types of boundary conditions are 
treated by a variety of methods. Chapter 5 is entitled ‘First-order equations,’’ and the reader may 
wonder what goes on, since the solution of a first-order equation cannot satisfy two independent con- 
ditions. A somewhat over-simplified explanation is the following. Obviously (4/4: — Yn-.)/2h is in general 
a better approximation to the first derivative y,’ than is (Yn41 — yn)/h, where h is the step interval. 
But when the first expression above replaces the derivative in the differential equation we are led to a 
second order differential equation. This differential equation has an extraneous oscillatory solution 
which must. be suppressed if possible. Nevertheless it may be easier to solve the second order difference 
equation by relaxation than to solve the original differential equation by step-by-step methods. Tech- 
niques for suppressing the oscillation are discussed. 

In Chapter 6, equations of higher order are examined. If central differences are used (as is often 
desirable because of greater accuracy for a given number of terms) the same problem of extraneous 
oscillatory solutions bobs up whenever the order of the equation is odd. This phenomenon is related 
to questions of stability for certain methods of solving differential equations step-by-step. 

Chapter 7 is devoted to eigenvalue problems associated with ordinary differential equations. When 
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a linear differential equation is replaced by a difference equation the eigenvalue problem is approx- 
imated by the problem of latent roots and vectors of a matrix. Hence Chapter 7 contains a rather full 
treatment of latent roots and vectors. In Chapter 8, Dr. Fox discusses various ways to use initial-value 
methods in order to solve two-point problems and eigenvalue problems, pointing out some advantages 
and some dangers of this procedure. In Chapter 9 he considers accuracy and precision of boundary value 
methods. In linear problems or ‘‘linearized’’ problems this involves essentially an examination of the 
latent roots of an appropriate matrix. Also involved are error estimates for the various formulas of 
finite differences used in the solution. The book concludes with miscellaneous methods and comments. 

The text abounds in interesting worked examples which not only illustrate the procedures being 
explained but also reveal the author’s skill and long practical experience in computational methods. 
As stated in the preface, the author has concentrated on methods ‘which are suitable for paper-and- 
pencil and desk-machine computation. Most of the methods are adaptable to automatic high-speed 
machines.” 

It is possible that the real strength of the book, namely its wealth of computational ideas, may be 
a disappointment to some readers who want to find a simple cut-and-dried procedure which can be 
immediately programmed for a high-speed digital computer. In fact the reviewer would like to have 
seen at least some space given to solutions on high speed machinery, even though he agrees with the 
author that a good grasp of paper-and-pencil methods should precede electronic methods 

All in all Dr. Fox has made a welcome and significant contribution to the literature of numerical 


methods, 
W. E. MILNE 


Progress in semiconductors. By Alan F. Gibson, R. E. Burgess, and P. Aigrain. John 
Wiley & Sons, Inc., New York, 1956. Volume I, viii + 220 pp. $8.00. Voluine IT, 
vii + 280 pp. $10.50. 


This consists at present of a two volume summary (Vol. I, 1956 and Vol. II, 1957) or review of 
work in semiconductor research and applications. It is intended that additional volumes be published 
annually. Each volume contains review articles written by specialists in the field. One of the outstand- 
ing features of these two volumes is the very large reference list accompanying each article. In view of 
the thousands of papers on semiconductor work published thus far, such a venture as this seems worth- 
while even if only for the reference lists. The articles are, however, excellent whether for specialist 
or non-specialist. These annual volumes are particularly recommended for anyone who wants to be 


aware of the main features of semiconductor work without reading hundreds of papers. 
R. TRUELL 


An introduction to probability theory and its applications. Volume I. By William Feller, 
John Wiley & Sons, Inc., New York, and Chapman & Hall, Ltd., London, 1957. 
xv + 461 pp. $10.75 


The first edition of this book, which is well known as a standard reference and popular textbook 
in probability theory, has been modified somewhat for this second edition. Although the style and the 
general outline of topics remain nearly the same, the detailed wording has been changed in many places. 
There is more emphasis on suggested topics for first reading and changes in wording appear to be moti- 
vated by a desire to improve the book’s usefulness as a text. Despite the fact that many sections have 
been made more concise, the second edition is somewhat larger than the first as a result of the addition 
of two new chapters on ‘Fluctuations in coin tossing and random walk’’ and ‘‘Compound distribu- 
tions, Branching processes.’’ This edition also puts greater emphasis on waiting times and the connec- 
tion between recurrent events and Markov chains. The high quality of the first edition is certainly main- 
tained if not, in fact, considerably improved. The book’s many admirers, however, continue to wait 


patiently for Volume II 
G. NEWELL 
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Mathematics and wave mechanics. By R. H. Atkin. John Wiley & Sons, Inc., New York, 
1957. xv + 348 pp. $6.00. 


This book consists of eight introductory chapters (198 pages) dealing with mathematical prepa- 
ration for the study of quantum mechanics. The remaining six chapters deal with quantum mechanics 
and some of the more important applications. The form of the book should make it particularly valuable 
to seniors in physics, applied mathematics or chemistry, who want to begin the study of quantum 
mechanics. It will also be valuable to graduate students who lack mathematical background for the 


study of quantum mechanics. 


R. TrRuew. 


Principles and techniques of applied mathematics. By Bernard Friedman. John Wiley & 
Sons, Inc., New York, and Chapman & Hall, Ltd., London, 1956. ix + 315 pp. $8.00. 


This book has the laudable intention “to show how the powerful methods developed by the ab- 
stract [mathematical] studies can be used to systematize the methods and techniques for solving prob- 
lems in applied mathematics.’’ More specifically, the abstract theory of linear spaces and linear operators 
on these spaces is used extensively in an attempt to systematize the use of Green’s function and eigen- 
function expansions in the solution of ordinary and partial differential equations. 

The first two chapters contain an excellent presentation of the theory of abstract linear spaces and 
operators which is required in the remainder of the work, and in this sense the book is self-contained. 
After introducing the Dirac 6-function in Chapter Three, the author introduces Green’s function for a 
second order linear ordinary differential equation and expresses solutions in terms of the Green’s func- 
tion. The eigenvalue problem for ordinary differential equations is considered in Chapter Four, and 
Chapter Five diScusses cursorily Laplace’s equation, the diffusion equation, and the wave equation. 

The reviewer feels that the strong emphasis on the theory of linear operators in the development 
of the theory of the last three chapters is obtained largely at the expense of rigour. For example, the 
proof of the completeness of the eigenfunctions u, = sin nz (n = 1, 2, ---) promised (on page 197) to 
be given later could not be found in the book by the reviewer. The formal approach, particularly in 
the last chapter, seems to obscure the underlying physical nature of the solutions to the equations 
considered. The reviewer also feels that the scope of the book is too limited for it to warrant such a 


broad title 
R. T. Suretp 


Proceedings of the Fifth Midwestern Conference on fluid mechanics. Edited by A. M. 
Kuethe. University of Michigan Press, Ann Arbor, 1957. viii + 388 pp. $8.00. 
The Volume collects 26 papers, among which a discussion of low-density Couette flow by H. T. 

Yang and L. Lees stands out as clearly deserving careful reading by any students of the slip-flow re- 

gime. T. Y. Li caleulates the hypersonic boundary layer corresponding to a prescribed pressure distri- 

bution, and W. N. MacDermott describes early operational experiences with an arc-heated blow-down 
hypersonic tunnel of about 0.03 sec. running time. H. H. Korst and W. Tripp predict the base pressure 
at the blunt trailing edge of a supersonic aerofoil at incidence; experiments confirm again the validity 
of the Korst model and suggest that an appreciable change in base pressure can arise from external 
pressure waves incident on the base region. A simplified model for the head-on interaction of shocks with 
rigid grids is developed by H. O. Barthel, and the characteristics systems of the equations governing 
interaction of a magnetic field with one-dimensional unsteady flow of a conducting gas are classified by 

S. I. Pai. R. W. Truitt proposes that the flow near the nose of a body in hypersonic flight be approxi- 

mated by that generated when the body moves slowly towards a fixed concentric shell. 

Some effects of chemical reactions on inviscid hydrodynamic stability are discussed by 8. F. Shen, 
and Gortler draws attention to H. Oser’s thesis on the vortex sheets arising in the low-frequency per- 
turbations of rotating inviscid fluids, and to the analogous effect in fluids with stable density stratifi- 
cation. M. Lessen gives the compressible perturbation equations for flow with vortex-symmetry, and 
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A. de Neufville gives solutions for vortices decaying due to viscous shear diffusion. Some functions 
useful for the caleulation of two-dimensional, incompressible laminar heat transfer from bodies with 
variable pressure and wall-temperature are tabulated by C. R. Guha and C. 8. Yih, J. M. Robertson 

r. Eckert and T. F. Irvine present some 


rbulent flow in rough pipes, while EK. R. ‘ 
W. Smith reports 


reviews data on turbul 
measurements on entrance length, transition and friction in non-circular pipes. D. 
A. San lborn proposes fn 


floating-element measurements of turbulent skin friction on a flat plate, V. 
and W. K. Mueller 


of mean profiles differing from Coles’ (which reviewer prefers), 


two-parameter iamlly 
discusses why slow pulsations in pipe flow should reduce the average turbulent heat transfer slightly 


below th it for ste 1dy mean flow 
P. Kaplan generalises Stokers work on the surface waves generated by pressure distributions to 
the cas iniformly moving oscillatory disturbances and shows that permanent advancing waves are 
ead of the disturbance when the product of its frequency and speed is less than g/4 
R. Seikel calculate surface wave effects on the performance Oo 

P. Breslin presents a simplified calculation of the 


Stelson and J. P. Murtha 


formed a far ah 
A. G. Strandhagen and G 
hydrofoil of infinite span by the lifting line model, J. 
suction exerted o1 very slender body by a prope ller at its stern, and F. EK. 
ss of some bodies moving in water. An ellipse-fitting approximation lor caleu- 


by the ‘“added-mass’”’ model is proposed by 


f a submerged 


measure the virtual n 
‘nsional impact of bodies on water 

‘ ee er > ; =e al 
reports experiments contiirming Ravleigh’s prediction of the instability of 


W. H. Guilinger and E. Saibel discuss the programming of the approximate 


of adiabatic slider bearings 


MI. Z. v. Krazywoblocki and C. G. Whittenbury 


lations 
Cy I ula | n ! 


As 


one-dl! A little philosophical accompaniment is pro- 
vid 


\IEYER 


Differential equations: geometric theory. By S. Lefschetz. Interscience Publishers, Inc., 


New York and London, 1957. x + 364 pp. $9.50 
Chis book is strongly to be recommended to anyone entering the domain of differential equations, 


either as a pure or applied mathematician 


The author has pres nted in a very readable and flowing 
tals of the theory of the existence and stability of equilibrium and peri dic solutions 


style the fundamer 
ar ordinary differential equations 
The title does not dc the contents justice In addition to presenting the geometric theory, in the 
ritical points, and periodic solutions, the author gives a very lucid and detailed 
Poinearé-Lyapunov. Starting with the 


of nonline 


sense of trajectories, 
ideas of these 


stability theory of 
Malkin, 


discussio: he al yt 
authors. he continues to the more recent work of Perron, Bellman, Levinson, Persidiski, 


Massera and Antosiewicz 


iddition to a quite detailed discussion of systems of n-th order, and then of the 


Re mark ibly, in 
second order, he presents a survey of various approximate techniques devised by Van der Pol and 
xtensively developed by Krylov-Bogoliubov, and of the more recent stroboscopic method of Minorsky 
These met ais 2 llustrated by examples 


RicHarv BeELLMAN 


Integraltefeln zur Quantenchemie. By H. Preuss. Volume II. Springer-Verlag, Berlin, 


Géttingen, Heidelberg, 1957. 143 pp. $8.60. 

This volume consists of two parts, in the first of which the integrals and the method of calculation 
are described along with the details of the one electron wave function used in the calculation of mole- 
cular energies. The second part of this volume is concerned with the numerical tabulation of the results 
of the calculation of the interaction integrals for the K and L Shells as a function of the screening param- 
eters used in the normalized Slater functions. 

These tables constitute the second volume of such integrals for quantummechanical calculations. 
This work is part of a program of the Max Planck Institute for Physics in Géttingen to prepare such 


tables for the more important types of calculation. 
RoHN TRUELL 
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